Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5790TG42.5% 128.9 / 10.2 81V289G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (5/41);  new base G (34/0);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.31e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  NC_000913/599434‑599727
                                                                                                                                                 |                                                                                                                                                     
gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                                  >  2:372479/1‑151 (MQ=255)
 acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                                 <  1:1049938/151‑1 (MQ=255)
  cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcc                                                                                                                                                  >  2:1336547/1‑149 (MQ=255)
   gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAttccg                                                                                                                                                >  1:856426/1‑146 (MQ=255)
    acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGggggtgctgc                                                                                                                                                  >  2:1211213/1‑143 (MQ=255)
    acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGc                                                                                                                                               >  2:1103095/1‑150 (MQ=255)
       ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAt                                                                                                                                            <  1:1273171/150‑1 (MQ=255)
          aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCgg                                                                                                                                        <  1:1300374/151‑1 (MQ=255)
             acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaa                                                                                                                                    <  1:1243163/151‑1 (MQ=255)
               aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGCTGCGGTcccggca                                                                                                                                     >  1:999257/1‑144 (MQ=255)
                     ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGcc                                                                                                                             <  1:671760/150‑1 (MQ=255)
                     ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCCGGAAAAACGccc                                                                                                                            >  2:1105009/1‑151 (MQ=255)
                      tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCTGCGAACCGGCCAAAACGcc                                                                                                                             >  2:847057/1‑149 (MQ=255)
                       aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTCCGGTCCGGCCAAAACGccc                                                                                                                            >  2:776305/1‑149 (MQ=255)
                        aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGccc                                                                                                                            >  2:226824/1‑148 (MQ=255)
                             cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTg                                                                                                                    <  1:77756/151‑1 (MQ=255)
                               gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGGTCCGGCAAAAACGCCCGAGAAGTg                                                                                                                    >  2:2640/1‑149 (MQ=255)
                                 gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                 <  1:556482/150‑1 (MQ=255)
                                  aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                                                                               <  1:423616/151‑1 (MQ=255)
                                  aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGCGGCGGGAGGGGGGGGGGTACTGAGAGCCGGCGAAAAACCCCCggaagagggag                                                                                                                 >  1:728255/1‑140 (MQ=255)
                                       ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                           <  1:125075/150‑1 (MQ=255)
                                       ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                           <  1:829709/150‑1 (MQ=255)
                                       ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                          <  1:371230/151‑1 (MQ=255)
                                           tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGACCTGCGGTCCGGCAAAAAAGCCCGAGGAGGGGGTGccgccgggg                                                                                                       >  1:239063/1‑147 (MQ=255)
                                              ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGTTCCGCAAAAAAGCCCCGGGAGGGATCGCCGCCGGGGAgg                                                                                                    >  2:1222590/1‑150 (MQ=255)
                                                gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACa                                                                                                 <  2:912459/151‑1 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGGAGGACaa                                                                                                >  1:364506/1‑150 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGCGGGCGGAGATGCTACGGCCCGGCGAAAACCCCCCAGCAGCGAACGCCGCCGGGAAgg                                                                                                    >  2:239120/1‑146 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGATCCCGCGAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGGCaa                                                                                                >  2:788563/1‑150 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGAATTGCCGAACTAAACGGCGAAGGCGAAGAGGCGGGGCGGGGGGGGATACTGCGATACGGGCAAAAATCCCGAGAAGGGTGCGCCGCCGCggaggaaa                                                                                                 >  2:871886/1‑147 (MQ=255)
                                                    aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGATCTGCGGTCCGGCAAAAAAGCCCGGGAAGGGGACGCCGCCGGGAAgggaaaa                                                                                               >  1:1098649/1‑144 (MQ=255)
                                                       tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGGCCTGTGATCCGGCCAAAACGCCCGGGGAGGGAACGCCGCCGGGTAAGAAAAAcatg                                                                                           >  1:634811/1‑147 (MQ=255)
                                                       tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCCGGAAGGGATCGCCGCCGCGAAGGAGaaaattg                                                                                           >  1:838543/1‑146 (MQ=255)
                                                       tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGGGATCCCGGCAAAACGCCCGAGCAGGGATAGCCGCCCTGAAGGACAAAATg                                                                                            >  2:740503/1‑149 (MQ=255)
                                                        gcggGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                             <  1:713202/147‑1 (MQ=255)
                                                        gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                             <  1:660065/147‑1 (MQ=255)
                                                        gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTTCCGCGATCCGGCCAAAAACCCCGCGGAGTGAACGCCGCCGTGGAGGACCAACAgg                                                                                           >  2:1236584/1‑149 (MQ=255)
                                                          gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGTATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGACCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGAACAAACTg                                                                                            <  2:600982/146‑1 (MQ=255)
                                                          gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGGAAAAAAGCCCGAGAAGTGATCGCCGCCGTGGAGGAAAAACTGGaaa                                                                                        >  1:673281/1‑150 (MQ=255)
                                                          gcgAGGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                       <  2:152441/151‑1 (MQ=255)
                                                           cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGATCCCGCAAAAAAGACCGAGGAGAGGTAGCCGCCCGGAAGGGCAAACTTGaaaa                                                                                       >  2:271219/1‑149 (MQ=255)
                                                           cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAACGCCCCCGAAGGGATCGCCGCCGGGAAAGACAAACTCGaaa                                                                                        >  2:1161958/1‑149 (MQ=255)
                                                           cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                        >  2:372391/1‑149 (MQ=255)
                                                            gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                        <  2:107181/148‑1 (MQ=255)
                                                               gTTGCGAAGGTCCAGATTGGCCCGGAGAGGCGCCGGGTCATTGTCGAACTAACCGGCGGAGGCGAAGTGGCGGGCGGGGTGGTGCTCCTGCGATTTGGCAAAAACGCCCGCGAAGTCATTGCCGCCGTGAATGACAAACTGGAAACGCTg                                                                                   <  1:862833/150‑1 (MQ=255)
                                                                  gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                               <  2:773227/151‑1 (MQ=255)
                                                                     aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                            <  2:846158/151‑1 (MQ=255)
                                                                     aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                            <  1:365800/151‑1 (MQ=255)
                                                                        gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCCCGCCAAAAAGCCCCAGAAGAGATTGCCGCCGGGGAGGGCAAACTGGAAACGCCGgaaa                                                                               >  1:413191/1‑145 (MQ=255)
                                                                        gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGGAACGCTGAAAAGTAGTc                                                                         >  2:343882/1‑151 (MQ=255)
                                                                           cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGTGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAAAGCCCCAGAAGGGATCGCCGCCGTGGAGGACAAAATGGAAAAGCTGAAAAGGAGTATGc                                                                      >  1:819869/1‑151 (MQ=255)
                                                                           cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                       <  2:1191379/150‑1 (MQ=255)
                                                                           cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGc                                                                      <  2:20470/151‑1 (MQ=255)
                                                                            aGATTGGCCCGGAGATGCGCCGGGACATTGCCGACCTAAACGGCAAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGGTCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                        <  1:197836/148‑1 (MQ=255)
                                                                             gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                                    <  2:263019/151‑1 (MQ=255)
                                                                               ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAAAGCCCGCGAAGGGCTCGGCGCCGGGGAGGACAAAATAGAAAAGCCGGAAAGAAGTCTggc                                                                     >  1:452502/1‑146 (MQ=255)
                                                                                 ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                                <  1:1165581/151‑1 (MQ=255)
                                                                                   cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAAGAGTCTGCCGGAAgg                                                               >  2:822277/1‑150 (MQ=255)
                                                                                      ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                            <  1:543529/150‑1 (MQ=255)
                                                                                      ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCCCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCg                                                             >  2:1132670/1‑149 (MQ=255)
                                                                                       gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCCCGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGGAggggggg                                                          >  1:934714/1‑146 (MQ=255)
                                                                                         gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                        <  2:1079058/151‑1 (MQ=255)
                                                                                         gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                         <  1:622883/150‑1 (MQ=255)
                                                                                            gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                      <  1:484110/150‑1 (MQ=255)
                                                                                            gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                      <  2:460656/150‑1 (MQ=255)
                                                                                            gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                     <  1:320252/151‑1 (MQ=255)
                                                                                                 ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGGAAAGGAGTCTGCCGGGAGGCGGGGAGATAGtt                                                  >  2:940067/1‑149 (MQ=255)
                                                                                                   ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                              <  1:174722/151‑1 (MQ=255)
                                                                                                   ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAACCCCCGAAGTGATCGCCGCCGTGGAGGGCAAACCGGGAACGCTGAAAAATAGTCTGCCGGGAGGCGGGGGGATAGGTaaa                                               >  1:543586/1‑148 (MQ=255)
                                                                                                    gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACCAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                             <  2:274392/151‑1 (MQ=255)
                                                                                                      aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCCGCGATCCCGGAAAAACCCCCGAGAAGGGATCGCCGCCGCGAAGGAGAAACTGGGAAAGCTGAAAAGCAGGCTGCCGGGAGGCGTggggggaggtaaaa                                              >  2:258366/1‑136 (MQ=255)
                                                                                                         gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                          <  2:6245/149‑1 (MQ=255)
                                                                                                          ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGATCCGGCAAAAAAGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCCGGAGGCGGGGAGATAGTTACAACATACg                                        >  2:1057436/1‑150 (MQ=255)
                                                                                                           cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                      <  1:1336480/151‑1 (MQ=255)
                                                                                                                  aaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGGAGGGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAATCTGCCGGGAGGCGTGGGGATAGTTACAACATACGATAGcagcc                               >  2:382179/1‑151 (MQ=255)
                                                                                                                     ggcgaaggcgaaGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGATAGTTACAACATACGATCGcagcca                              >  2:143851/1‑149 (MQ=255)
                                                                                                                       cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                          <  2:766418/151‑1 (MQ=255)
                                                                                                                        gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                          <  1:1154217/150‑1 (MQ=255)
                                                                                                                         aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                        <  2:299155/151‑1 (MQ=255)
                                                                                                                           ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                       <  2:634830/150‑1 (MQ=255)
                                                                                                                           ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                      <  1:980536/151‑1 (MQ=255)
                                                                                                                                      gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc           >  1:1214037/1‑151 (MQ=255)
                                                                                                                                                gtgATCATGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTc  >  2:456224/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  NC_000913/599434‑599727

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: