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

GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  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)
             acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑aa                                                                                                                                    <  1:1243163/151‑1 (MQ=255)
               aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGCTGCGGTcccggca                                                                                                                                     >  1:999257/1‑144 (MQ=255)
                     ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGcc                                                                                                                             <  1:671760/150‑1 (MQ=255)
                     ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCCGG‑AAAAACGccc                                                                                                                            >  2:1105009/1‑151 (MQ=255)
                      tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCTGCGAACCGGC‑CAAAACGcc                                                                                                                             >  2:847057/1‑149 (MQ=255)
                       aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTCCGGTCCGGC‑CAAAACGccc                                                                                                                            >  2:776305/1‑149 (MQ=255)
                        aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGC‑AAAAAAGccc                                                                                                                            >  2:226824/1‑148 (MQ=255)
                             cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTg                                                                                                                    <  1:77756/151‑1 (MQ=255)
                               gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGGTCCGGC‑AAAAACGCCCGAGAAGTg                                                                                                                    >  2:2640/1‑149 (MQ=255)
                                 gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATc                                                                                                                 <  1:556482/150‑1 (MQ=255)
                                  aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATcgc                                                                                                               <  1:423616/151‑1 (MQ=255)
                                  aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGCGGCGGGAGGGGGGGGGGTACTGAGAGCCGGCGAAAAACCCCCggaagagggag                                                                                                                 >  1:728255/1‑140 (MQ=255)
                                       ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCgccgcc                                                                                                           <  1:125075/150‑1 (MQ=255)
                                       ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCgccgcc                                                                                                           <  1:829709/150‑1 (MQ=255)
                                       ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGccgccg                                                                                                          <  1:371230/151‑1 (MQ=255)
                                           tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGACCTGCGGTCCGGC‑AAAAAAGCCCGAGGAGGGGGTGccgccgggg                                                                                                       >  1:239063/1‑147 (MQ=255)
                                              ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGTTCCGC‑AAAAAAGCCCCGGGAGGGATCGCCGCCGGGGAgg                                                                                                    >  2:1222590/1‑150 (MQ=255)
                                                gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACa                                                                                                 <  2:912459/151‑1 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGGAGGACaa                                                                                                >  1:364506/1‑150 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGCGGGCGGAGATGCTACGGCCCGGC‑GAAAACCCCCCAGCAGCGAACGCCGCCGGGAAgg                                                                                                    >  2:239120/1‑146 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGATCCCGC‑GAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGGCaa                                                                                                >  2:788563/1‑150 (MQ=255)
                                                  ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGAATTGCCGAACTAAACGGCGAAGGCGAAGAGGCGGGGCGGGGGGGGATACTGCGATACGGG‑CAAAAATCCCGAGAAGGGTGCGCCGCCGCggaggaaa                                                                                                 >  2:871886/1‑147 (MQ=255)
                                                    aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGATCTGCGGTCCGGC‑AAAAAAGCCCGGGAAGGGGACGCCGCCGGGAAgggaaaa                                                                                               >  1:1098649/1‑144 (MQ=255)
                                                       tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGGCCTGTGATCCGGC‑CAAAACGCCCGGGGAGGGAACGCCGCCGGGTAAGAAAAAcatg                                                                                           >  1:634811/1‑147 (MQ=255)
                                                       tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGC‑AAAAAAGCCCCGGAAGGGATCGCCGCCGCGAAGGAGaaaattg                                                                                           >  1:838543/1‑146 (MQ=255)
                                                       tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGGGATCCCGG‑CAAAACGCCCGAGCAGGGATAGCCGCCCTGAAGGACAAAATg                                                                                            >  2:740503/1‑149 (MQ=255)
                                                        gcggGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                             <  1:713202/147‑1 (MQ=255)
                                                        gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                             <  1:660065/147‑1 (MQ=255)
                                                        gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTTCCGCGATCCGGC‑CAAAAACCCCGCGGAGTGAACGCCGCCGTGGAGGACCAACAgg                                                                                           >  2:1236584/1‑149 (MQ=255)
                                                          gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGTATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGACCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGAACAAACTg                                                                                            <  2:600982/146‑1 (MQ=255)
                                                          gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGG‑AAAAAAGCCCGAGAAGTGATCGCCGCCGTGGAGGAAAAACTGGaaa                                                                                        >  1:673281/1‑150 (MQ=255)
                                                          gcgAGGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                       <  2:152441/151‑1 (MQ=255)
                                                           cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGATCCCGC‑AAAAAAGACCGAGGAGAGGTAGCCGCCCGGAAGGGCAAACTTGaaaa                                                                                       >  2:271219/1‑149 (MQ=255)
                                                           cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGC‑AAAAACGCCCCCGAAGGGATCGCCGCCGGGAAAGACAAACTCGaaa                                                                                        >  2:1161958/1‑149 (MQ=255)
                                                           cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                        >  2:372391/1‑149 (MQ=255)
                                                            gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                        <  2:107181/148‑1 (MQ=255)
                                                               gTTGCGAAGGTCCAGATTGGCCCGGAGAGGCGCCGGGTCATTGTCGAACTAACCGGCGGAGGCGAAGTGGCGGGCGGGGTGGTGCTCCTGCGATTTGGC‑AAAAACGCCCGCGAAGTCATTGCCGCCGTGAATGACAAACTGGAAACGCTg                                                                                   <  1:862833/150‑1 (MQ=255)
                                                                  gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                               <  2:773227/151‑1 (MQ=255)
                                                                     aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                            <  2:846158/151‑1 (MQ=255)
                                                                     aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                            <  1:365800/151‑1 (MQ=255)
                                                                        gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCCCGC‑CAAAAAGCCCCAGAAGAGATTGCCGCCGGGGAGGGCAAACTGGAAACGCCGgaaa                                                                               >  1:413191/1‑145 (MQ=255)
                                                                        gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGC‑AAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGGAACGCTGAAAAGTAGTc                                                                         >  2:343882/1‑151 (MQ=255)
                                                                           cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGTGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGC‑AAAAAAGCCCCAGAAGGGATCGCCGCCGTGGAGGACAAAATGGAAAAGCTGAAAAGGAGTATGc                                                                      >  1:819869/1‑151 (MQ=255)
                                                                           cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                       <  2:1191379/150‑1 (MQ=255)
                                                                           cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGc                                                                      <  2:20470/151‑1 (MQ=255)
                                                                            aGATTGGCCCGGAGATGCGCCGGGACATTGCCGACCTAAACGGCAAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGGTCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                        <  1:197836/148‑1 (MQ=255)
                                                                             gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                                    <  2:263019/151‑1 (MQ=255)
                                                                               ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGC‑CAAAAAGCCCGCGAAGGGCTCGGCGCCGGGGAGGACAAAATAGAAAAGCCGGAAAGAAGTCTggc                                                                     >  1:452502/1‑146 (MQ=255)
                                                                                 ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                                <  1:1165581/151‑1 (MQ=255)
                                                                                   cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAAGAGTCTGCCGGAAgg                                                               >  2:822277/1‑150 (MQ=255)
                                                                                      ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                            <  1:543529/150‑1 (MQ=255)
                                                                                      ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCCCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCg                                                             >  2:1132670/1‑149 (MQ=255)
                                                                                       gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCCCGC‑AAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGGAggggggg                                                          >  1:934714/1‑146 (MQ=255)
                                                                                         gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                        <  2:1079058/151‑1 (MQ=255)
                                                                                         gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                         <  1:622883/150‑1 (MQ=255)
                                                                                            gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                      <  1:484110/150‑1 (MQ=255)
                                                                                            gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                      <  2:460656/150‑1 (MQ=255)
                                                                                            gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                     <  1:320252/151‑1 (MQ=255)
                                                                                                 ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGC‑AAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGGAAAGGAGTCTGCCGGGAGGCGGGGAGATAGtt                                                  >  2:940067/1‑149 (MQ=255)
                                                                                                   ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                              <  1:174722/151‑1 (MQ=255)
                                                                                                   ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGC‑AAAAAAACCCCCGAAGTGATCGCCGCCGTGGAGGGCAAACCGGGAACGCTGAAAAATAGTCTGCCGGGAGGCGGGGGGATAGGTaaa                                               >  1:543586/1‑148 (MQ=255)
                                                                                                    gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACCAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                             <  2:274392/151‑1 (MQ=255)
                                                                                                      aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCCGCGATCCCGG‑AAAAACCCCCGAGAAGGGATCGCCGCCGCGAAGGAGAAACTGGGAAAGCTGAAAAGCAGGCTGCCGGGAGGCGTggggggaggtaaaa                                              >  2:258366/1‑136 (MQ=255)
                                                                                                         gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                          <  2:6245/149‑1 (MQ=255)
                                                                                                          ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGATCCGGC‑AAAAAAGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCCGGAGGCGGGGAGATAGTTACAACATACg                                        >  2:1057436/1‑150 (MQ=255)
                                                                                                           cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                      <  1:1336480/151‑1 (MQ=255)
                                                                                                                  aaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGC‑AAAAAAGCCCGAGGAGGGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAATCTGCCGGGAGGCGTGGGGATAGTTACAACATACGATAGcagcc                               >  2:382179/1‑151 (MQ=255)
                                                                                                                     ggcgaaggcgaaGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGATAGTTACAACATACGATCGcagcca                              >  2:143851/1‑149 (MQ=255)
                                                                                                                       cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                          <  2:766418/151‑1 (MQ=255)
                                                                                                                        gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                          <  1:1154217/150‑1 (MQ=255)
                                                                                                                         aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                        <  2:299155/151‑1 (MQ=255)
                                                                                                                           ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                       <  2:634830/150‑1 (MQ=255)
                                                                                                                           ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                      <  1:980536/151‑1 (MQ=255)
                                                                                                                                      gggcgggGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc           >  1:1214037/1‑151 (MQ=255)
                                                                                                                                                gtgATCATGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTc  >  2:456224/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGC‑AAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  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: ‑