Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,4320TG59.6% 3.5 / 51.3 104L159R (CTG→CGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base T (10/28);  new base G (38/18);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 110 total aligned reads displayed.

ACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTG  >  NC_002947/4160302‑4160562
                                                                                                                                  |                                                                                                                                   
ccctgTACCGTGCTGATGGGCGTGTGGGGGGGAATGAACGGGTGCATGCGCGACCGCGATACCGCCTTGCTGACGGCCGCGGGCAGCGCTAGCTGCGGCGCCGCCGGGGTGCTGGCATTCGAATGGGCGCTgcgct                                                                                                                                <  2:1804154/132‑2 (MQ=255)
aCCTGTACCCTGCTGAGCGGCCTGTGGTGTGTCATTAATGTTTTAAAGCTCGACCGCTATACCGCCTTGCTGCCCGCCGCCGGCAGCGCCTTGGGGGGCGCCGCCGCGGGGGTGGCATTCGAATCggcgcggcgcc                                                                                                                                <  2:2765969/136‑2 (MQ=255)
aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTgcgcc                                                                                                                                >  2:2428767/1‑135 (MQ=255)
aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCa                                                                                                                                <  1:1965366/136‑1 (MQ=255)
aCCAGTACCCTGCTGATCGGCCTGGGGGGTGGCATGAAAGTTTTCAAGCTCGACCGCTATACCGCCTTGCTGACCGCCGCCGGCAGCTCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCACTGCGCa                                                                                                                                <  1:3846376/136‑1 (MQ=255)
  cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGTGCTGCGCAgc                                                                                                                             >  2:3943871/1‑136 (MQ=255)
  cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGcgccgc                                                                                                                             >  2:1020361/1‑136 (MQ=255)
  cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAgc                                                                                                                             >  1:954411/1‑136 (MQ=255)
  aaGTACCCGGCTGAGCGGCGTGGGGGGTGGAATCAATGTGTTCAAGCGCGACCGCGAAACCGCCTGGCTGCCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAgc                                                                                                                             <  2:2351381/135‑1 (MQ=255)
     tACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGGCGCGcc                                                                                                                          >  1:2670087/1‑136 (MQ=255)
       cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGGGCGCCGcccccc                                                                                                                        >  2:1587071/1‑136 (MQ=255)
          tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCCTGCTGGCATTCGAATCGGGGCGGCGGACccccccccc                                                                                                                     >  1:23450/1‑135 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGGGGCGccccccccaa                                                                                                                  >  1:2876918/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGCCGccccccccccca                                                                                                                  >  1:1361701/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGCCGccccccccccaa                                                                                                                  >  1:2557199/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGCCCCGccccccccaa                                                                                                                  >  1:3915852/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGGCGCGccccccccca                                                                                                                  >  1:4330262/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCGGcccccccccccc                                                                                                                  >  1:4446926/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCGcccccccaaa                                                                                                                  >  1:3280408/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCACccccccccccaa                                                                                                                  >  1:1823031/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAGaa                                                                                                                  >  1:1796691/1‑136 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCCCGcccccccaaa                                                                                                                  >  2:3523703/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCCCCGcccccccccaaa                                                                                                                  >  1:3103387/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATAGGCGGAGCGCcccccccccccaaa                                                                                                                  >  1:4685683/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAATCGAATCGGCGCGGCGGACCGcccccccaaa                                                                                                                  >  1:3060416/1‑132 (MQ=255)
             tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGAGCTGGCATTCCAATCGGCGCAGCGACGCGgcccccccaa                                                                                                                  >  2:3027755/1‑132 (MQ=255)
               aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCCGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGCGCCCCGcccccccccccaca                                                                                                                >  1:1922618/1‑130 (MQ=255)
               aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGCCGCGCCCCCCCCaaaa                                                                                                                >  2:3336178/1‑136 (MQ=255)
               aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCGGCGgccccccccaca                                                                                                                >  1:3575165/1‑130 (MQ=255)
               aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCACCGCCCCCCCCaaaa                                                                                                                >  1:156609/1‑136 (MQ=255)
               aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGGCGCGCCCCCCCCaaaa                                                                                                                >  1:200441/1‑136 (MQ=255)
                  ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCCGCCCCCCCCCCaaaaaaa                                                                                                             >  1:382472/1‑133 (MQ=255)
                  ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCAGCGCCGCGCCCCCCAGaaaaacc                                                                                                             >  2:1021567/1‑133 (MQ=255)
                  ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCaaaaaaccc                                                                                                             >  2:942789/1‑133 (MQ=255)
                  ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCACCCCCCCCCCAaaaaacc                                                                                                             >  2:2599722/1‑133 (MQ=255)
                    ccTGAGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCCGCCCCGCCGCCCTGCTGGCATTCGAATCGGCGCGGCCACGCACCCCCCCCaaaaaacgc                                                                                                           >  2:528415/1‑131 (MQ=255)
                       gtggtgtgtCATTAATGTTTTCAAGCTCGCCGGCTAAACGGCGTTGCTGACCGCCGCGGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGGTGGCTTTCGAATCGGGGCTGCGCTCCGCCCCCCATAAAAGCGCCATg                                                                                                        <  2:2314840/136‑1 (MQ=255)
                       gtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGGGCCGCGCCCCCCAGAAAAGCGCCaag                                                                                                        >  2:4643363/1‑134 (MQ=255)
                       gtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCAGCGCCCCCCAAAAAAGCCCCagg                                                                                                        >  2:4437544/1‑134 (MQ=255)
                         ggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCCCCCCCCCAAAAAACCCCAAggg                                                                                                      >  2:870095/1‑135 (MQ=255)
                         ggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAACGCCATGGc                                                                                                      >  2:4404687/1‑136 (MQ=255)
                         ggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGGGCTGGCATTCGCATCGGCGCGGCGGCGCCGCCCCCACACAAACGcccaggc                                                                                                      >  2:2857641/1‑131 (MQ=255)
                            gtggCATTAAAGTGTTCAACCCCGACCGCGAACCCGCCTTCCTCCCCGCCGCCGGCCGCGCCCTCTGCGGCGCCGCCGCCGTGCGGGCATTCGATCCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGt                                                                                                   <  2:2235802/136‑1 (MQ=255)
                            gtggCATGATTGTGTGCACGCGCGCCCGCCATCCCGCGTGGCGCCCCGCCGCGGCCCGCGCCTGCTGCGGCGCCGCCGCGGGGCTGGCATTCGACGGGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGt                                                                                                   <  2:164876/136‑1 (MQ=255)
                            gtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGACGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGt                                                                                                   <  1:4678263/136‑1 (MQ=255)
                                cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCGCCCCCCCCCCAAAAAACGCCCAGGGAGTCggg                                                                                               >  1:4612958/1‑135 (MQ=255)
                                cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGGCCAGCGCCCCCCAAAAAAACGCCCTAGCAgggcgc                                                                                               >  1:1954996/1‑131 (MQ=255)
                                cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGAGCTGGCATTCGAATCGGCGCAGCGCACCCCCCCCCCAAAAAGCCCCAAAGCACTCGGa                                                                                               >  2:2306609/1‑136 (MQ=255)
                                cATGAAAGTGTTCAAGCTCGACAGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGTAGCGCCCCCCATAAAAGCGCCAAGGCAGtagca                                                                                               >  2:872071/1‑132 (MQ=255)
                                 attaaaGTGTCCACGCGCCCCCGCGATCCCGCTTTGCCCCCCGCCGCCGCCCGCGCCATCGGCGGCGCCCCCGCGGTGCTGCAATTCTAAGCGGCGGGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                              <  2:1768025/133‑1 (MQ=255)
                                 aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGaa                                                                                              >  2:3566035/1‑136 (MQ=255)
                                     aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGGGCGCCGCCGCCCCCCAAAAAACGCCAAGGGAGTGGGGAggca                                                                                          >  2:293038/1‑133 (MQ=255)
                                     aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCCGAAAACCGCCAAGGCAGGCGGGAGCgg                                                                                          >  2:3468635/1‑135 (MQ=255)
                                      tgttttCACGCTCGACCGCTATACCGCTTGGCGGACCCGGCCCGGCAGCGCCTGCTGCGGCGCCGCCGCCGGGCTGTAATTCTAAGCGGCTCTGCGCTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                         <  2:2720092/132‑1 (MQ=255)
                                      tgtgtTCAACCTCGCCGGCTATCCCGCCTTGCTGACCGCCGCCGCCTGCGCCTTGTGGGGCGCCGCGGCGGTGGTGGAATTCCAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                         <  2:3820433/135‑1 (MQ=255)
                                      tgtgcTCACGCGCGCCCGCGATCCCGCCTTGCCCCCCGCCGCCGGCCGCGCCCTCGGCGGCGCCGCCGGCGTGGTGTAATTCGACCCGGCGCGGCGCCCCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                         <  2:3327530/135‑1 (MQ=255)
                                      aGTTTTCAAGCTCTCCCGCTAAACCGCCTTGGTCACCCCCGCCGGCAGCGCCTTCGGGGGCGCCGCCGCCGTGCTGGCATTCTAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                         <  2:3583257/136‑1 (MQ=255)
                                      aGTGTTCACGCGCGACCGCGATACCGCCTTGCGGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGGTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                         <  2:1720330/136‑1 (MQ=255)
                                           tCAAGCTCCCCCCCGATCCCGCTTTGCTGCCCCCCGCCGGCAGGGCCATGTGCGGCGCCGCCGCGGTGGTGTCATTCCAAGCGGCGCGGCGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                    <  2:2910334/136‑1 (MQ=255)
                                             aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCGGCGCCCCGCCCCCCCAAAAAGCGCCCAGGCCGTCGGCAGGGTAGttgggg                                                                                  >  2:3969755/1‑131 (MQ=255)
                                             aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCCGCGCCCCCCCCAAAAACCCCCAGGGACGCGGAAGCGGAggggggt                                                                                  >  1:2268517/1‑130 (MQ=255)
                                             aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCCGCGCCCCCCCCAAAAAACCCCAAGCAGGCGGAAGCGGAgtggtgg                                                                                  >  1:959430/1‑135 (MQ=255)
                                             aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCAGCGCCCCCCCAAAAAACGCCAAGGGAGGCGGGAGGGGAGTGgtgt                                                                                  >  1:2703704/1‑136 (MQ=255)
                                             aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCCCCCCCCCCAAAAACGCCCAGGCAGGCGGACGCCTAggtgggg                                                                                  >  1:1364019/1‑130 (MQ=255)
                                                  cGACCGCGATCGCGCCTGGCTGCCGGCCGCCGCCAGCCCCATCGGCGGCGCCGCCGCGGTGCTGGAACTCCACGGGGGCGTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                             <  2:4207468/136‑1 (MQ=255)
                                                       gcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCGCCCCCCAAAAAAAAGCCCAGGCGCGCGGGAGCGTAGTGGTGTTggggggcgt                                                                        >  2:3659428/1‑130 (MQ=255)
                                                       gcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAGCCCCCAGGCAATCGGAAGGGTATTGCTGGTCGGttactt                                                                        >  2:3094289/1‑131 (MQ=255)
                                                           tcccGCCTTGCTGCCCGCCCCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCATTCTAAGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                    <  2:4579434/134‑1 (MQ=255)
                                                           tACCGCGTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCGGTGCTGGATTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                    <  2:1697542/136‑1 (MQ=255)
                                                           ccccgccgTGCCGCCCGCCGCCGGCCGCGCCCGCGGCGGCGCCGCCGCCGGTGTGCAATTCGACGCGGCGCGGCCCACCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                    <  1:1712140/134‑1 (MQ=255)
                                                           ccccGCGTTGCTGCCCGCCGCGGGCAGCGCCATCTGGGGCGCCGCCGCGGTGCGGGCATTCGAAGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                    <  2:3766324/134‑1 (MQ=255)
                                                           acccGCCTTGCTGCCCGCCCCCGGCAGCGCCAGCTGGGGCGCCGCCGCCGTGCTGGCATTCTACTCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                    <  2:1023357/134‑1 (MQ=255)
                                                                   tGCTGACCGCCACCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGTGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                    >  2:2699281/1‑128 (MQ=255)
                                                                      tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCCGCGCCCCCCCAAAAAGCCCCCAGGCAGTCGGAAGCGTAGTGGTGGTAGGGTCCGTGTCGATTTTACTTTa                                                         >  1:1838582/1‑136 (MQ=255)
                                                                          cgccgccgGCCGCGCCAGGGCCGGCGCGGCCGCCGGGCTTGCATAGGACGCGGCGCTCCCCCCCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                     <  1:2569172/136‑1 (MQ=255)
                                                                          cgccgccgGCAGCGCCTGCGGCGGCGCCGCCGCGGTGCTGGCATTCGAAGCGGCGCGGCGCCCCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                     <  1:4154252/136‑1 (MQ=255)
                                                                               ccgGCAGCGCCAGCGGCGGCGCCGCCGCCGTGCTGGATTTCGACGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGc                                                <  2:3743973/136‑1 (MQ=255)
                                                                                             gggggcgCCGCGGCGGTGCTGGCATTCGACGCGGCGCGTCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                  <  2:4004587/133‑1 (MQ=255)
                                                                                             gggggcgCCGCCGCGGGGCTGGAATGCTACGCGCCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                  <  2:3982240/133‑1 (MQ=255)
                                                                                             ggcggggCCGCGGCGGGGGTGGCTTTCTATCCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                  <  2:2546666/135‑1 (MQ=255)
                                                                                             ggcggcgCCGCCGCCGTGCTGGCTTTCGAGGGGGCGCGGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                  <  1:4626034/135‑1 (MQ=255)
                                                                                             ggcggcgCCGCCGCCGTGCTGGCTTCCCCCGCGGCGCGGCCCCCCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                  <  1:1236444/135‑1 (MQ=255)
                                                                                               gggcgCCGCGGCGGTGGTGGCTTTCGAAGGGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                <  1:1746202/135‑1 (MQ=255)
                                                                                               gggcgCCGCCGGGGTGCTTGCCTTCGCCGCGGCCCGGCGCCCCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                <  1:1020361/135‑1 (MQ=255)
                                                                                               gggcgCCGCCGCCGTGCTGGCATTCTTATCGGCGCGGCGGAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                <  2:3060416/135‑1 (MQ=255)
                                                                                               gggcgCCGCCGCCGGGGTGGGATTCGAAGCGGCGCGGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                <  2:3854353/135‑1 (MQ=255)
                                                                                                gccgacgccgcggGGCTGGAATTCTATGCGGCTCTGGTCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAATCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACACACTGGCGATCAATGCCggctggc                               <  2:2253383/131‑1 (MQ=255)
                                                                                                  cgccgccgcCGTGCTGTCTTCCGACCCGGCCCGGCGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGgctggctg                             <  2:4007259/136‑1 (MQ=255)
                                                                                                   gccgccgccGGGCTGGAATTCTATCGGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGc                            <  2:382472/136‑1 (MQ=255)
                                                                                                     cgccgccgTGCTGCTATTCGCCGCGGCGCGGCGCCCCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCAt                          <  1:3855071/136‑1 (MQ=255)
                                                                                                       ccgccgTGCTGGCATTCTAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                        <  2:493199/136‑1 (MQ=255)
                                                                                                        ggccgGGGTGGCTTTCGATGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                       <  2:2400070/135‑1 (MQ=255)
                                                                                                                        gAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg       <  2:2773664/136‑1 (MQ=255)
                                                                                                                        caagcggcgcggcgCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg       <  2:286844/132‑1 (MQ=255)
                                                                                                                          aTCGGCGCTGCGCAGCGCCCCCCAAAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTTCCATGTAGATGTTACTTTAAACACTGGCGGTAAAATCCGGCTGGCTGCATCAGGGTACCAAGGGCGCCggg     >  2:3517485/1‑136 (MQ=255)
                                                                                                                             ggcgcggcgcTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg  <  2:1800458/136‑1 (MQ=255)
                                                                                                                             ggcgcggcgcCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg  <  2:4239127/136‑1 (MQ=255)
                                                                                                                             ggcgcggcgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg  <  2:801557/136‑1 (MQ=255)
                                                                                                                             ggccgtGCGCTGCCACCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg  <  1:838181/131‑1 (MQ=255)
                                                                                                                             ggCGCTCCGGCGCCCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg  <  1:248880/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
ACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTG  >  NC_002947/4160302‑4160562

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: