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

GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATC  >  NC_000913/599434‑599718
                                                                                                                                              |                                                                                                                                               
gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                         >  2:372479/1‑151 (MQ=255)
 acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                        <  1:1049938/151‑1 (MQ=255)
  cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcc                                                                                                                                         >  2:1336547/1‑149 (MQ=255)
   gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggggggcc                                                                                                                                         >  1:516662/1‑142 (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)
                                                                                                                                              |                                                                                                                                               
GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATC  >  NC_000913/599434‑599718

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: