Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG50.7% 99.2 / 123.4 70T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (2/32);  new base G (35/0);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.19e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-14

GTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356058‑1356349
                                                                                                                                                |                                                                                                                                                     
gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTtaaa                                                                                                                                                >  2:775133/1‑150 (MQ=255)
gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTtaaa                                                                                                                                                >  2:1207722/1‑150 (MQ=255)
gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAt                                                                                                                                                >  2:25478/1‑151 (MQ=255)
 tCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAAt                                                                                                                                                <  1:1021458/150‑1 (MQ=255)
 tCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATg                                                                                                                                               >  2:974464/1‑151 (MQ=255)
   aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGcc                                                                                                                                             <  2:739973/151‑1 (MQ=255)
       gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTaaaatcccgaa                                                                                                                                         >  2:134987/1‑143 (MQ=255)
        cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGAt                                                                                                                                         <  2:1231213/150‑1 (MQ=255)
            cGAGGCGCGCGGTAAAAAAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCg                                                                                                                                    <  2:1527286/151‑1 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAAGGCCGGAAAGGGGCGTTCAAcc                                                                                                                         >  2:1026639/1‑150 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCAGATAGGCGGCGTtcctcc                                                                                                                         >  2:332100/1‑150 (MQ=255)
                       gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGCGGGGGTAAATGCCAGATAGGGGACGGTCATcc                                                                                                                         >  2:678574/1‑150 (MQ=255)
                        tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATccc                                                                                                                        <  1:1188931/150‑1 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                   <  2:553050/151‑1 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                   <  2:1379004/151‑1 (MQ=255)
                            aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                   <  2:1618224/151‑1 (MQ=255)
                                     gccgTTTCAGCCGTACCATCATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGCGGTAAATGCCAAATAAGGGGCGCTTATCCCCGGGCGCAGCGGTa                                                                                                          >  2:1583612/1‑151 (MQ=255)
                                        gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAAATCCAGGTAGGCGGCGGTTACCCCCGGAAGCAGCGGTaaaa                                                                                                       >  2:451627/1‑150 (MQ=255)
                                         tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAAc                                                                                                       <  2:168452/150‑1 (MQ=255)
                                          ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGc                                                                                                     <  1:451477/151‑1 (MQ=255)
                                          ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGc                                                                                                     <  1:285951/151‑1 (MQ=255)
                                                  tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAAGCGGCGTTCATTCCCGGACGCAGCGGTCAACGCAAGCGCgg                                                                                             >  1:1460321/1‑151 (MQ=255)
                                                  tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAAAGCCAGAAAGGCGGCGTTTCTACCCGGACGCCGCGCTAAACGCAAGAGCgg                                                                                             >  1:1131988/1‑151 (MQ=255)
                                                    ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAATCCCCATAGGGGGCGGTTCTCCCAGGACGCAGCCTTAAAAGCAAACGggggg                                                                                           >  2:822870/1‑150 (MQ=255)
                                                    ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAAAGGCAGAAAGGGGGCGGTCCTCCCCGGACGCCGCGGTTAAAGCaaaagggggg                                                                                           >  2:1621383/1‑143 (MQ=255)
                                                    ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAAAGCCAGAAAGGGGGCGTTTATTCCCGGACGCAGCGGTTAAAGCAAGCGCgggg                                                                                           >  1:680035/1‑150 (MQ=255)
                                                    ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                           <  1:1185389/151‑1 (MQ=255)
                                                    ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAATGCCCGATAGGCGCCGCTCATCCCCGGACGCCGCGGTTAAAGCAAACGggggg                                                                                           >  2:1615114/1‑150 (MQ=255)
                                                     cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtc                                                                                          <  1:1614335/151‑1 (MQ=255)
                                                      ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGCTAAAAGGCAGAAAGGCGGAGTTTATTCCCGGGCGCAGCGGTTAAAGCAAGAGgggggtg                                                                                         >  2:800629/1‑148 (MQ=255)
                                                      ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAAGGGACGGTCATCCCCGGGCGCAGCGGTAAACGCAAGCGCGGggcg                                                                                         >  2:9280/1‑151 (MQ=255)
                                                       taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAATGCCCGACAAGGGGCGGTTATCCCCCGGCGCAGCGGTAAAAGCCAACGCGGggcgg                                                                                        >  1:1608082/1‑150 (MQ=255)
                                                       taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGAAAGGCGACGGTCATCCCCGGACGCAGCGTTTAACGCAAACGCggggggg                                                                                        >  2:1085972/1‑147 (MQ=255)
                                                        aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCAGAAAAGCGCCGTTCATCCCCGGACGCAGCGGTAAACGCAAACGCGGggcgt                                                                                        >  2:1366046/1‑150 (MQ=255)
                                                          taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTTAATTGCAGATAGGGGCCGGTTCACCCAGGGCGCCGCGGTAAACGCCAAAGCGGGTCGTCa                                                                                      >  1:883790/1‑150 (MQ=255)
                                                              tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaa                                                                                 <  1:752449/151‑1 (MQ=255)
                                                               ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAATGCCAGAAAGGGGGCGGTCATCCCAAGACGCAGCGGTTAAAGCAAACGCGGGGCGGCCCGCaaa                                                                                >  2:1563962/1‑151 (MQ=255)
                                                               ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCCGATAGGGGAAGTTTATCCCAGGGCGCAAGGGTAAAAGCCAGCGCGGGGCGGCACGCaaa                                                                                >  2:521047/1‑151 (MQ=255)
                                                                    gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAACCCCAATATGCGACGGTCATCCCCCCGCCCGGCGCTAAACCCAAACGCGGGGGGTCACCCCAAAAggg                                                                           >  2:564116/1‑150 (MQ=255)
                                                                     cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAAAGGCAGATAGGCGGCGCTCATCCCAGGACGCCGCGTTCAACCCAAGCGCGGGGGGTCCCGCCAAAtaggt                                                                          >  1:939145/1‑147 (MQ=255)
                                                                     cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCCAATAGGCGGCGGTCATCCCCGGACGCCGCGGTTAAAGCAAACGCGGGGCGTCACGCAAAAAggga                                                                          >  2:1690331/1‑149 (MQ=255)
                                                                           gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGCCCGAAAGGGGACCTTTATCCCCGGGCGGCGCGGTAAACACAAGCGCGGGGGGTCACGCAAAAATGATAggggg                                                                    >  1:1406053/1‑147 (MQ=255)
                                                                           gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctg                                                                     <  2:1611107/150‑1 (MQ=255)
                                                                            cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggc                                                                   <  1:1656777/151‑1 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAATTCCCAAAAGGCGACGTTCATCCCCGGGCGCCGCGTTAAACGCAAGCGCGGGGGGGCACCCCAAAAGGATAGccggctggg                                                               >  2:945198/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCCGTGGTTAATGCCAGAAAGGCGGCCGTCCTCCCCCGGCGCCGCGGTTAACACCAGCGCGGGGCGTCCCCCCAAAAGGATAgcctgcctg                                                                >  1:227476/1‑143 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGAAAGGGGGCGTTCATCCCCGGACGCAGCGGTTAAAGCAAACGCGGGGCGTCACGCCAAAAGGGAAGatggcggg                                                                >  2:1575834/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGCGGGGTTTAAAGCCAGAAAGGGGACGGTCTTCCCAGGACGCAGCGGTAAACGCAAGCGCGGGTGGGCACGCCAAAAGGGGAGctggcgggg                                                               >  1:1270224/1‑150 (MQ=255)
                                                                                ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGCGGGGGTTAATGCCCGCTAAGGGGCGGTTATCCCCGGACGCAGCGGTTAAAGCCAACCCGGGGCGTCCCGCCAAAAAGAGAGctgggtggg                                                               >  1:1252954/1‑150 (MQ=255)
                                                                                 gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCAAAAAAGCGGCGTTTCTCCCCAGACCCCGCGTTTAAAGCCAGCGCGGGGCGGCACGCAAAAATGGTAgctgggcgggg                                                              >  1:1374484/1‑145 (MQ=255)
                                                                                  ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                             <  1:353208/151‑1 (MQ=255)
                                                                                  ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                             <  2:410354/151‑1 (MQ=255)
                                                                                        aTCGCCGGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcc                                                       <  2:905423/151‑1 (MQ=255)
                                                                                        aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcc                                                       <  1:1541555/151‑1 (MQ=255)
                                                                                            ccAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggccaggc                                                   <  2:1192303/151‑1 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTAAAAGCCCAATAGGGGACGTTCATCCCCGGACGCCGCGTTAAAACCAAACGCGGGGGGTCCCGCAAAAAGGAAAGCCGGGTGGGCGcagggccgggccga                                               >  2:1215901/1‑141 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:903251/148‑1 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  1:858294/148‑1 (MQ=255)
                                                                                                tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGGTAAATGGCAGATAGGCGGCGGTCATTCCATGGCCCAGCGCTTAAAGCCAGCGCGGCGCGTGACACAAAATCGAAAGATTGCGGGGAGCaggcgcggcccag                                               >  1:1554347/1‑148 (MQ=255)
                                                                                                  ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                             <  1:18112/151‑1 (MQ=255)
                                                                                                  ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                             <  1:438608/151‑1 (MQ=255)
                                                                                                  ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                             <  1:334026/151‑1 (MQ=255)
                                                                                                                    cGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGt                              <  2:1643615/148‑1 (MQ=255)
                                                                                                                      ggATTATTTAGCGGCGGTTTTGCGGTGGTTAATGCCAAAAAGGCGGCCGTTATCCCCGGACGCAGCGGTTAACGCCAGCGCGGGGCGGCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCCGGCCAGGACGGCGCCATCCCCGGCCAggg                         >  1:76252/1‑150 (MQ=255)
                                                                                                                      ggATTATTTAGCGGCGGTTTGGCGGGGGTTAATGCCAGATAGGGGACGTTCATCCCCGGACGCCGCGGTAAAACCCAACGCGGGGCGGCCCGCCAAATGGATAACGTGGTTGCAGGAGGCCAGGCCCGAACGGCGCAGTCCCCCGTCCggg                         >  1:899660/1‑150 (MQ=255)
                                                                                                                        attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTTGCATTCCCCGGTCAGGAg                        <  2:1092343/150‑1 (MQ=255)
                                                                                                                        attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAg                        <  1:1550959/150‑1 (MQ=255)
                                                                                                                            tttAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCg                   <  1:526902/151‑1 (MQ=255)
                                                                                                                                    cggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg            <  1:256571/150‑1 (MQ=255)
                                                                                                                                             cGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgc    <  1:1316219/150‑1 (MQ=255)
                                                                                                                                               gtgtTAAATGCCAGATAGGCGACGTTCATCCCAGGGCGCAGCGTTACACGCAAGCGCGGGTCGTCACCCAAAATGGATAGCTTGCTCGCACCATGCCAGGACAGAACGTCCCATTCCCCGGTCAGGAGTTTCTACAGACGTCCGGTgcgg  >  1:1012055/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356058‑1356349

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