Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG48.6% 99.2 / 20.2 73T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (4/33);  new base G (35/0);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.93e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356058‑1356351
                                                                                                                                                |                                                                                                                                                             
gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGGGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTtaaa                                                                                                                                                       >  2:1161465/1‑150 (MQ=255)
 tCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGAGGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAAt                                                                                                                                                       <  2:1673791/150‑1 (MQ=255)
  cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                                     <  1:1300782/151‑1 (MQ=255)
  cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                                     <  2:380582/151‑1 (MQ=255)
  cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                                     <  2:699497/151‑1 (MQ=255)
          cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAATCCAGAGAgg                                                                                                                                             >  2:1608029/1‑151 (MQ=255)
            cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCg                                                                                                                                           <  2:531637/151‑1 (MQ=255)
                     cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAAAGCCCGATAGGGGGCGCTCAt                                                                                                                                   >  2:1420127/1‑150 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAATGTCAAAAAgggggagtttcttc                                                                                                                                 >  1:633271/1‑138 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGGCAGAGAGGCGAAGGTTATcc                                                                                                                                 >  1:809093/1‑150 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAATACCAGATAGGGGGCGGTCATcc                                                                                                                                 >  1:369382/1‑150 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATcc                                                                                                                                 <  1:998540/150‑1 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAATGCCAGAAAAGGGACGGtcctcc                                                                                                                                 >  2:1160719/1‑150 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGTTTAATTCCAGATAGGCGGCGTTCATcc                                                                                                                                 >  2:1727290/1‑150 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGATTTGGGGGGTTAAAAGCCAGAAAGGGGGCGTTAAttc                                                                                                                                 >  2:1380158/1‑148 (MQ=255)
                        tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATccc                                                                                                                                <  1:908431/150‑1 (MQ=255)
                         aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGCCCGATAAGGGGCGGTCATccccg                                                                                                                              >  2:338215/1‑149 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGGCAGAAAAGGGACGTTTATCCCCgggg                                                                                                                           >  1:1650438/1‑149 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                           <  2:1088357/151‑1 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                           <  1:626261/151‑1 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                           <  1:1010815/151‑1 (MQ=255)
                            aTAAAGGCCGCCGTTTCAGACGTACCAGAATAGATGGATTGCATCAGGCGCTGGTTATTAAGCGCCAGTGCAAGCCCATGCCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGCTAAATGCCAGCAAGGCGAAGTTCAGCCCAAGAc                                                                                                                           <  1:870675/151‑1 (MQ=255)
                                      ccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCCGCGTGTAATGCCAGATAGGGGACCTTTATTCCACGgcgccccgctca                                                                                                                 >  1:271700/1‑143 (MQ=255)
                                         tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAAc                                                                                                            <  2:1715579/150‑1 (MQ=255)
                                         tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAAc                                                                                                            <  1:1495636/150‑1 (MQ=255)
                                         tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACg                                                                                                           <  1:1658050/151‑1 (MQ=255)
                                              gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGCCAGGTAGGCGGCGGTTAACCCCGGACGccacggttaaagccaa                                                                                                          >  2:1488186/1‑135 (MQ=255)
                                                 gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTGAAAGCCCCAAATGCGGCGCTTATCCCAGGATGCGGTGTTAAACGCAAAcgcg                                                                                                  >  1:584047/1‑151 (MQ=255)
                                                  tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAATGCCCGATAGGCGGCGGTTATCCCAAGAACCAGCGGTTaaaacaaaagcgg                                                                                                   >  1:1624627/1‑141 (MQ=255)
                                                  tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAAAGCCAGATAAGGGGCGGTTATCCCCGGGAGCCGCGGTAAAAGCAAACGCgg                                                                                                 >  1:1544021/1‑151 (MQ=255)
                                                  tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAATGCCAGAAAGGCGACGTTCATCCCCGGACGCAGCGGTAAAAGCAAACGCgg                                                                                                 >  2:1675491/1‑151 (MQ=255)
                                                    ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAAAAAGGGGGAGTTCATCCCAGGACGCAGCGGTTAAAGCAAAAGCgggg                                                                                               >  2:422410/1‑150 (MQ=255)
                                                      ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAATGCCAGATAAGGGACGGTTATCCCCGGGCGCCGCGGTAAAAGCAAACGCGGggcg                                                                                             >  1:1358930/1‑151 (MQ=255)
                                                      ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGGCAAATAGGGGACGTTCATCCCCCGGCGCAGCGATTAACCACAAGGCgggggg                                                                                             >  2:513896/1‑148 (MQ=255)
                                                      ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCCCATAGGGGACGTTCTTCCCCTGACCCCGCGGTTAAAGCCAACGCGGggcg                                                                                             >  2:1358262/1‑151 (MQ=255)
                                                       taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAATGCCAGATAGGGGGCGGTCCTCCCCGGACGCAGGGGTTAAAGCAAACGCGGggcgg                                                                                            >  2:310577/1‑150 (MQ=255)
                                                          taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGCCAGAGAAGCGACGTTCATCCCCGGGCGCAGCCGTTAAAGCAAACGCcgggtggca                                                                                          >  1:1341130/1‑144 (MQ=255)
                                                               ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAAAGTCAGAAAAGGGAGGTTTATCCCACGACACAGCGTTCAAAACAAAGGCGGGGCGTCAAGCaaa                                                                                   >  2:551985/1‑151 (MQ=255)
                                                                    gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGCGGGGTTTAAAGCCCGAAAAGGGGCGGTCCTCCCCGGACTCAACGGTTAAAGCAAACGCCGGGCGTCCCGCaaaaaaga                                                                              >  1:861101/1‑147 (MQ=255)
                                                                          ggCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAATGCCCGATAGGCGGCGTTCATCCCAGGACGCAGCGGTAAACGCAAACGCGGGTCGGCACGCCAAATgggaag                                                                           >  1:1409670/1‑144 (MQ=255)
                                                                           gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTTGTCACGCAAAATGGATAgc                                                                          <  2:409673/148‑1 (MQ=255)
                                                                           gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgc                                                                          <  1:301416/148‑1 (MQ=255)
                                                                              cTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCACAAAGGGGACGGTCATCCCCAGCCGCAGCGGTAAAAGCAAACGCGGGGCGGCCCGCcaaaaagaaagccgg                                                                       >  1:1231609/1‑137 (MQ=255)
                                                                              cTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCCCAAAAGGGACGCTCATCCCCCGACGCACCGTTTAAAGCAAACGCGGGTCGTCCCGCAAAAAGGGTAactgg                                                                       >  1:1699250/1‑148 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGTTTAATGCCAGAAAGGGGACGTTCATCCCAGGACGCCGCGGTAAAAGCAAGCGCGGGGCGGCACCCAAAAAGGGAAGctggctggg                                                                  >  2:1255180/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGTTAAATGACCGAAAGGGGGCGGTTTTCCCAGGAAGCAGCGTTAAAAACAAACAGGGGGGGGCACCCAAAAAGGAAAgcacggcggg                                                                  >  2:1495638/1‑143 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCCGAAAAGGGGCGTTCATCCCCGGACGCAGCGTTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGAAAGctggctggg                                                                  >  2:1213220/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCCGATAGGCGGCCTTTCTCCCCCGACGCCGCGTTTAAAGCCAACGCGGGGGGGCACCCAAAAATGGTAGctgggtgg                                                                   >  2:1589394/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGCGACGTTCATCCCCGGACGCCGCGGTTAAAGCAAGCGCGGGGCGTCACGCAAAAAGGGTAGctggctggg                                                                  >  2:760926/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGCGACGTTCATCCCAGGGCGCAGCGGTAAACGCAAACGCGGGGCGGCACCCCAAAAGGATAGctggctgg                                                                   >  2:1175704/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTGTTAAGCCAGCTAGGCGACCGTTATCCCACGACGGCGCGTCAAATTACACCCCGGGCCCTGCCGCAAAAATGATAgctgggggg                                                                   >  1:1421884/1‑146 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATTCCCGAAAGGCGACGTTTATCCCCGGACCCCGCGTTTAAAGCAAGCGCGGGTCGTCACGCCAAAAGGGTAAccgggtggg                                                                  >  1:1663972/1‑150 (MQ=255)
                                                                                     ttaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGCGGCGGTCATCCCCGGACGCAGCGTTAAAACCAAGCGCGGGGCGTCACCCCAAAAGGGTAACCGGCGGGgagcag                                                             >  2:483992/1‑151 (MQ=255)
                                                                                       aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggc                                                           <  2:1444754/151‑1 (MQ=255)
                                                                                               gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGGTTAATTCCAAATAGGCGGCCTTCATCCCAGGACGCAGCGGTTAACGCCAAGCCGGGGCGCCCCCCCAAATGGGTAACTCGCCTGCCGCaggccaggcccg                                                   >  2:138812/1‑149 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                     <  2:489024/148‑1 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                     <  2:793/148‑1 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                     <  2:797196/148‑1 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                     <  1:1510306/148‑1 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                     <  1:1716926/148‑1 (MQ=255)
                                                                                                ggCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTGTGGCGGTGTTAAATGCCAGATAGGCAACGTTCACCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                     <  1:294661/147‑1 (MQ=255)
                                                                                                  ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                                <  1:1651533/151‑1 (MQ=255)
                                                                                                  ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                                <  1:533246/151‑1 (MQ=255)
                                                                                                  ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                                <  1:310480/151‑1 (MQ=255)
                                                                                                  ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCGGGCTGGCAGCAGGCCAGGCCAGAAc                                                <  2:450797/151‑1 (MQ=255)
                                                                                                           ggCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGGTAAAAGGCCGGAAAGGGAAGGTTCTTCCCAGGCGGCGGGGTTTAAGCAAAGGCGGGGGGCTAAGCAAAAATGGATGGGTGGCGGGAGCAGACCCCGCCAGAACGgcgccatt                                       >  1:1310814/1‑147 (MQ=255)
                                                                                                                    cGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                              <  2:1578583/150‑1 (MQ=255)
                                                                                                                    cGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGAATTCCCCGGTc                               <  2:14496/149‑1 (MQ=255)
                                                                                                                       gATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGAAAGGCGACGTTTATCCCCGGGCGCAGCGTTAAACGCAAGCGCGGGGCGCCACGCAAAATGGATAGCTGGCGGGGAGCAGGCCAGGCCAGAAAGCCCCATTTCCCCggccgggg                           >  2:83208/1‑145 (MQ=255)
                                                                                                                             ttAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                    <  2:1525186/151‑1 (MQ=255)
                                                                                                                              tAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc                   <  1:1479974/151‑1 (MQ=255)
                                                                                                                                    cggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg              <  1:1658339/150‑1 (MQ=255)
                                                                                                                                    cggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg              <  2:90415/150‑1 (MQ=255)
                                                                                                                                                tgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGAGAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgcc  <  1:1260288/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                             
GTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356058‑1356351

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