Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,148	0	G	A	25.5%	105.2 / 12.6	51	E28K (GAA→AAA)	fabR	DNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (34/4);  new base A (1/12);  total (35/16)
Fisher's exact test for biased strand distribution p-value = 1.35e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCC  >  NC_000913/4160900‑4161388                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  gctttgTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa                                                                                                                                                                                                                                                  >  1:321201/1‑250 (MQ=255) gctttgTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa                                                                                                                                                                                                                                                  >  1:297048/1‑250 (MQ=255)   tttgTCTGATTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGAGGCAGGTTTCCGCTGTACGTAACAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTCTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCCGTACTTTTGGCAAGGA‑TTCAGACATCGTCATGGGCGTAAGAGCGCAACAAAAAGaaa                                                                                                                                                                                                                                                >  1:171136/1‑250 (MQ=255)            ttCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTCCTTTTGTCAAGGA‑TTCAGCCACCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                                                                                                                                                                                                                                       <  1:222418/250‑1 (MQ=255)                   ttttGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtg                                                                                                                                                                                                                               >  2:247955/1‑251 (MQ=255)                               ggtttATTGGGATTTTATCGTCGATAATCTCCGTTTTAGCGGCAGGTTTTCGGTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCCGTAAATATGGTTTTTTGCGTTACCGT‑TCATTCACAAAACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTTGCAAGGG‑GTCAGTAATCGTGGTGTGCGTAAGAGCGCAACCAAAAAAAAAAACCCCCCGCTCTCTGGTGGAAGCCGcc                                                                                                                                                                                                                   >  2:219117/3‑250 (MQ=255)                                  ttattgtttttttATCAACGATCCCTTCCGTTTAAGCGGCAGTTTTCCGCTGTCCGTAAAAGACCCGG‑CCAAAAAAT‑TGCAGTAAATTTTTTTTTTGGCGTTCCCGT‑TCAGTCAAAATATTGGA‑GAAATCCATTATGTTCTTTCTCTGT‑TATAGTGCCAGCAGTCCTTTTGGAAAGGA‑TTCAGACATCGTAAGGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAtt                                                                                                                                                                                                                 <  1:397406/242‑1 (MQ=255)                                               aTCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg                                                                                                                                                                                                   >  1:389359/1‑251 (MQ=255)                                                  gacgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                                                                                                                >  1:358611/1‑251 (MQ=255)                                                  gacgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                                                                                                                >  1:293037/1‑251 (MQ=255)                                                  aatggaggacgaccaatatccGGCAGTTTGCCGCGGTAGGAAAAAGACCGGG‑CAAAAGAAT‑TGCAGAAAATATGTTTTATGGCGTGCCCGT‑CAATTCCCAATACTGGA‑GCAATCCAGTATGTCCATCCTCTGG‑TAAAGTCCCAGCAGTACTTTTGGGAAGTA‑TGCAGACATCGTGATGGGCTTAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCt                                                                                                                                                                                                 <  1:163757/230‑1 (MQ=255)                                                   acgaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                                                                                                                                                                                                >  2:53741/1‑250 (MQ=255)                                                     gaTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGAAAGGA‑TTCAGACATCGTGAGGGGCGTAAGAGCGCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAAc                                                                                                                                                                                             <  2:389359/251‑1 (MQ=255)                                                      aTACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACg                                                                                                                                                                                            >  2:220506/1‑251 (MQ=255)                                                        aCTCTCCGTGTAAGCGGCAGTTTTCGGCTGTACTAAAAAAAACCGG‑C‑AAAAAATGTGCATAAAATATTTTTTTTTGCGTT‑GCGTGTAATCCAAAATA‑TGGAAGCAATCCGTTATGTTCATCCGTGGG‑TATAGGGCAAGCAGTACTTTGGGCAAGGT‑TTCAGACGGGGGGAGGGGCGTAAGGGCGCAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa                                                                                                                                                                                          <  2:38812/251‑1 (MQ=255)                                                                   aaGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc                                                                                                                                                                               >  2:231723/1‑251 (MQ=255)                                                                     gcgggtggTTTCCGTTTTACGAAAAAGACCCGG‑CAAAAAAAT‑TGCATTAAAAATTTTTTATTGCGTTACCTT‑TAATTCAAAATACTGGA‑GCACTCAAGTTTTTTCATTCGCTGG‑TATGGTGCCAGCAGT‑CCTTTGGGAAGGAATTCAGACTTCGGGGGGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCt                                                                                                                                                                             <  1:44157/245‑1 (MQ=255)                                                                              tccgtcggTACGTAAAAGA‑CCGGACAAAAGATT‑TGTAGTAAATATGTTTAATGGCGTTACCGT‑TCATTCAAAATACTGGA‑GAAACCCGGTGTTTTCATTCTCGGG‑TATAGGGCCGGCATTATTGGTGGCAAGGA‑TTCAGACATCGTGTTGGGCGTAAGAGCGTAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGc                                                                                                                                                                     <  1:170925/243‑1 (MQ=255)                                                                                               gAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGa                                                                                                                                                     >  2:291749/1‑249 (MQ=255)                                                                                                   cgggcaaaaGAAT‑TGCGGTAAATATGTTTTTTTGCTTTCCCGT‑TCATTCAAAATACGGGA‑GCAAACAAGTTTGTTCATTTTCTGG‑TATAGTGCCAGAAGTACTTTTGGAAAGGA‑TTCAGACATGGGGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCgg                                                                                                                                                <  2:181044/245‑1 (MQ=255)                                                                                                                  gCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCAcc                                                                                                                                  >  2:28921/1‑251 (MQ=255)                                                                                                                         atatGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCttttt                                                                                                                          >  2:15553/1‑251 (MQ=255)                                                                                                                                ttATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGc                                                                                                                    >  2:184118/1‑250 (MQ=255)                                                                                                                                 tATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCAt                                                                                                                  >  1:202140/1‑251 (MQ=255)                                                                                                                                      cGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGAATTGCTCCCACC‑TCTTTTTATCGGCAttttcg                                                                                                             >  2:239986/1‑248 (MQ=255)                                                                                                                                      cGTTACAGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCcg                                                                                                             >  1:444655/1‑251 (MQ=255)                                                                                                                                          cccGT‑TCATTCACAATACGGGA‑GCAATCCCGTATGTTCTTTCTCTGT‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGAc                                                                                                         <  1:231723/250‑1 (MQ=255)                                                                                                                                           ccGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACg                                                                                                        >  2:41170/1‑251 (MQ=255)                                                                                                                                             gT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGt                                                                                                       >  1:142691/1‑250 (MQ=255)                                                                                                                                              t‑tcattcaCAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGt                                                                                                       >  1:173942/1‑249 (MQ=255)                                                                                                                                                  attcaCAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAA‑CCCCCCGGTCGGTGGGGGGAGCCGCCTTTAACCAATTTAGTGCTTAACCCAGCTTCGGCAGCCTGAGGTTGCGTGGAGTGGCGCCTGAAGCGGGCAATGCTCCCACCCTCttttaacggccattccgccgcgtagacgg                                                                                                  >  1:271464/1‑226 (MQ=255)                                                                                                                                                   ttcaCAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGAGGTAGACg                                                                                                  >  2:459406/1‑250 (MQ=255)                                                                                                                                                       caATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg                                                                                             >  1:274068/1‑251 (MQ=255)                                                                                                                                                       aaaTACGGGA‑GCAATCCAGTTTGTTCATTC‑CTGGGTATAGTGCCAGCAGTATTTTGGGCAAGGA‑TTAAGACATCGGGAGGGGCGTAAGAGCGAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg                                                                                             <  1:210424/250‑1 (MQ=255)                                                                                                                                                        aaTACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCAGCGACGTAGACGAAcggg                                                                                            >  2:181498/1‑248 (MQ=255)                                                                                                                                                             gtgA‑GCAATCCAGTTTGTTCATTCTCTGG‑TATAGTGCAAGCATTCTTTTTGGCAAGAA‑TTCAGCCATCGTGATGGCCGTAAGGGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt                                                                                        <  1:35450/248‑1 (MQ=255)                                                                                                                                                                  gCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCTAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAcc                                                                                    >  1:305939/1‑250 (MQ=255)                                                                                                                                                                        cAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtg                                                                             >  1:166273/1‑251 (MQ=255)                                                                                                                                                                               ttCATTCTCTGG‑TATAGTGCAAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAACACCCGCCGTTCGCTGGTGGAAGCAGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCACCC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                                                                        >  2:99908/1‑249 (MQ=255)                                                                                                                                                                                 cATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCCCC‑TCTTTTTATCGGCATTTCCGCTACGTAGACGAACTGGGTCTGTCCATGGTTGATGAGAGCg                                                                    >  2:429976/1‑251 (MQ=255)                                                                                                                                                                                        cggttataGTGCCGGCAGTACTTTGGGCAAGGA‑TTCAGACTTCGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGtt                                                                 <  1:291749/244‑1 (MQ=255)                                                                                                                                                                                          g‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGc                                                           >  2:358372/1‑251 (MQ=255)                                                                                                                                                                                                           aCTTTTGGAAAGGA‑TGCAGACATCGTGAGGGGCGAAAGAGCGAAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATgc                                            <  1:38211/250‑1 (MQ=255)                                                                                                                                                                                                                tGGCAAGGA‑TTAAGACATCGGGAGGGGCGTAAGAGCGCAACAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAg                                       <  1:41170/250‑1 (MQ=255)                                                                                                                                                                                                                     aGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGt                                 >  1:404763/1‑251 (MQ=255)                                                                                                                                                                                                                      ggA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGACCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCACCTCATGCGCCAGGCGCGTc                                >  2:97611/1‑251 (MQ=255)                                                                                                                                                                                                                               acatggggaggggCGTAAGAGACCAACAAAAAAAAAAACCCCGCCGGTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATc                        <  1:201929/241‑1 (MQ=255)                                                                                                                                                                                                                                 acggtgaggggtgaaaGAGCACAAAAAAAAAAAAATCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAAGGCTACGCCAACTCATGCGCCATGCGCGTCAGCGTATCGc                      <  2:145229/237‑1 (MQ=255)                                                                                                                                                                                                                                      gATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCaa                   >  1:278725/1‑249 (MQ=255)                                                                                                                                                                                                                                              aaggatcacaaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAATGCGGGAgt         <  1:308531/239‑1 (MQ=255)                                                                                                                                                                                                                                               aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAgtg        >  2:210242/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gagCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGtgtg      >  1:314083/1‑251 (MQ=255)                                                                                                                                                                                                                                                     aaaaaaaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATcc  <  2:9421/246‑1 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  GCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGG‑CCAAAGAAT‑TGCAGTAAATATGTTTTATTGCGTTACCGT‑TCATTCACAATACTGGA‑GCAATCCAGTATGTTCATTCTCTGG‑TATAGTGCCAGCAGTACTTTTGGCAAGGA‑TTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACC‑TCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCC  >  NC_000913/4160900‑4161388

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