Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,386	0	G	A	34.9%	48.7 / 17.8	43	A89A (GCG→GCA)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base G (24/4);  new base A (0/15);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 2.56e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAA  >  NC_000913/3486137‑3486621                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            gCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCg                                                                                                                                                                                                                                             >  1:450995/1‑250 (MQ=255)                       tGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTg                                                                                                                                                                                                                       >  1:31926/1‑250 (MQ=255)                         gTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGaa                                                                                                                                                                                                                     >  1:242575/1‑250 (MQ=255)                                          acaTTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAttt                                                                                                                                                                                                   >  1:561716/1‑251 (MQ=255)                                                       cccATCCAAGAGCACCTTTATTCACCAGGGTGAAAAAGCGGAACCGCTGTACTACACCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAAAAA‑TGATCCTCCCTTATCTGAACCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGGGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTc                                                                                                                                                                                      <  1:132680/251‑1 (MQ=255)                                                        ccATCCAAAGGCGCGCTTATGCACCAGGGTGAAAAAGGGAAAACCCTGTACTACATCGTAAAAGGCTCTGTGGAAGTGCTGAAAAAAGACGAAGAGGGAAAAAAAT‑TGTTCCTCCCCTATCTGAATTGGGGGTATTTTATT‑GG‑CGAAGGGG‑G‑CC‑GTTTGGAAGAGGG‑CCAGGAAC‑GTAGCGCATGGTT‑ACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCa                                                                                                                                                                                     <  1:78529/251‑1 (MQ=255)                                                         caccaaaaagcccGCTTATTCACCAGGGTAAAAAAGAGGAAACGTGGTACTACTTCTTTAAAGGCTCTGTGGCATTGTTTACAAAATACGAAGAGGGAAAAGAAA‑TGATCCTCTCCTATCTTAATC‑GGGTTATTTTTTTTGG‑CTAACGGG‑G‑CCTGTTTG‑AAGAGGG‑CCGGTAAT‑GTGGCGCAGGGGT‑CTTTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAg                                                                                                                                                                                    <  1:490967/239‑1 (MQ=255)                                                                 agCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATAAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑T‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGACTGGGT‑ACGTGAGAAAACCGCCTGTGAAGTGGCTGCAATTTAGTACAAAAAATTTCGCCATTTGATTCAgacaaacc                                                                                                                                                                             >  2:470870/1‑243 (MQ=255)                                                                          aTTCACCAGGGTAAAAAAGCGGAACCGCAGTACTACATCTTTAAAGTCTGTGTGGCAGGTCGAATCAAAGAAGAAG‑GGGTAAAAAAAATGACCCTCTCCTATTGGAATCGGGGTTATTTTTTG‑GG‑CGAACGGG‑G‑CTTGTTTG‑AAGGGTG‑CAAGGAAC‑GTGGCGCATGGGT‑ATGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTTCAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCt                                                                                                                                                                   <  1:502447/251‑1 (MQ=255)                                                                              aCCAGGGTGAAAAAGCGTAAACGTTGTACAACTACGTTAATGTCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGACCCTCTTCTATTTGAATCAGGGTGATTTTATT‑GG‑CGACCGGG‑G‑CCTTTGTA‑AAGGGGG‑CAAGAAAC‑GTAGCGAATGGGT‑CCGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATg                                                                                                                                                               <  1:92496/251‑1 (MQ=255)                                                                              aCCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATg                                                                                                                                                               >  1:506606/1‑251 (MQ=255)                                                                               caaGGGTAAAAAAGCGGAAACCCTTTACTAAATCGTTAAAGTCTCGGTGGCAGTGTTAATAAAAGACGAAGAGGGAAAAGAAA‑TGACCCCCTCCTATCTAAATCAGGTTTTTTTTTGT‑GG‑CGAAC‑GG‑G‑CCTTTTTGGAAGAGGG‑CCAGGAAC‑GTAGCGC‑TGGGTTACTTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGc                                                                                                                                                              <  1:487712/249‑1 (MQ=255)                                                                                cAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGTAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTAGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCg                                                                                                                                                             >  1:70804/1‑251 (MQ=255)                                                                                  gggTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACAGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAAGTAAACCCGGACATTCTGATGCGtt                                                                                                                                                           <  1:291382/251‑1 (MQ=255)                                                                                             cGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGcac                                                                                                                                                  >  1:255885/1‑249 (MQ=255)                                                                                                            acgtctatAAAGGCTTTGTGGCAGTTCTAAACAAAGACGAAGAGGTTAAAAAAA‑TTCCCCTCTCCTATCAGAATCGGGGTGATTTTATT‑GG‑CGAA‑TGG‑GACTTTTTTA‑AAGAGCG‑CCAG‑AACAGTAGCACAGGGGT‑ACGGGAAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTg                                                                                                                                 <  2:161612/244‑1 (MQ=255)                                                                                                                  ttAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTc                                                                                                                           >  1:493958/1‑251 (MQ=255)                                                                                                                    aaaGGCTCTGTGGCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGA‑CCAGGAAC‑GTAGCGCATGGGT‑ACGGGCGAAAACCGCCTGTGAAGAGGCTGAAATTTCGTAAAACAAATTTCGCCAATTGATTCAGGTAAAACCGGACAATCTGGTGCGTTTGTTTGCACAGGAGGCGCATACTCTGCAAGTccc                                                                                                                         >  2:359306/1‑249 (MQ=255)                                                                                                                      gggctcggggggagggctgaaaaaaGACGAAGGTGTTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGGGTTTTTATT‑GG‑CAAGCTGG‑G‑CCGTTTTG‑AAGGGGG‑CCAGAAAC‑GTAGCGAATGGGT‑ACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACtt                                                                                                                       <  1:381111/229‑1 (MQ=255)                                                                                                                       ggCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTAAAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTTTGCACAGATGGCGAGTCGTCGGAAAGTCCCTTc                                                                                                                      >  2:26608/1‑251 (MQ=255)                                                                                                                       ggCTCTGGGGCAGTGCTGATAAAAGACAAAGAGGTAAAAGAAA‑TGATCCCCTCCTATCTAAATCAGGGTGATTTTATG‑GG‑CGAACTGG‑C‑CCTTTTTG‑AAGGGGG‑CCAGGAAC‑GTAGCGAATGGGT‑ACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATGTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTc                                                                                                                      <  1:353309/251‑1 (MQ=255)                                                                                                                               ggCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGATGG‑CCAAGAAC‑GAAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGACGTGGCTGAAATTTGGTACAAAAAATTTCGCCAGTTGATTCAGGTAACACCGGACATTCTGCTGCGTTTGTCTGCACAGATGGCGCGACGTCTGCACGTAACTTAAgagag                                                                                                                >  1:20953/1‑248 (MQ=255)                                                                                                                                  agtgcgaaacaaaaacaaagatgatAAACAAA‑TTATCCTCTCCTATCTGAATCGGGGTGATTTTATT‑GG‑C‑AACTGGTG‑TCTTTTTG‑AAGTGGC‑CCAGAAAC‑GTCGCGCTTGGTT‑CCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTgg                                                                                                           <  1:375624/227‑1 (MQ=255)                                                                                                                                     gCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGc                                                                                                         >  1:486019/1‑250 (MQ=255)                                                                                                                                       tGATCAAAGACGAAGGGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGCTTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCaa                                                                                                       >  2:529452/1‑250 (MQ=255)                                                                                                                                           aaaaaaggAAGAGGAAAAAGAAA‑GGACCCTCCCCTATCGAAATAGGGGTTATTTTATG‑GG‑GGACCGGG‑G‑CTTGTTGG‑AAGGGGG‑CAAGGAAG‑GTAGCGCAGGGGT‑AGGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTgg                                                                                                  <  1:277534/244‑1 (MQ=255)                                                                                                                                            aaaGAAGACGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATAAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AACAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGGGACGGGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCCGGGAACCCCGGACATTCTGATGCGTTTGTCTGCACCGACGGCGCGACGTCTCCAAGTAACTTCAGAGTAAGTGTGCAACATggg                                                                                                 >  2:530342/1‑250 (MQ=255)                                                                                                                                                aggAAGAGGGTAAAAAAA‑GGATCCTCCCCTATCTGAATAGGGGTGATTTTATG‑GG‑GGAACTGG‑G‑CTTGTT‑G‑AAAAGGGGCCGGAAAC‑GAGGCGCATGGGT‑CCGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTc                                                                                             <  1:498625/249‑1 (MQ=255)                                                                                                                                                          tAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGAt                                                                                                                                                                >  1:426600/1‑174 (MQ=255)                                                                                                                                                          tAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTg                                                                                    >  1:345182/1‑250 (MQ=255)                                                                                                                                                                             cTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGcaa                                                                    >  1:282216/1‑247 (MQ=255)                                                                                                                                                                               aTCTGAAACAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCAAATTGATTCAGGTAAACCCGGACATTCTGATGAGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTAACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCCCAGAct                                                               >  2:234472/1‑251 (MQ=255)                                                                                                                                                                                   gAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTctgc                                                           >  1:108384/1‑251 (MQ=255)                                                                                                                                                                                   gAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGAGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCAATTCAGAGAAGTTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTctgc                                                           >  2:189331/1‑251 (MQ=255)                                                                                                                                                                                             taatttATT‑GGGCGAACGGG‑G‑TTTGTTTG‑AAGCGGG‑CCATGAAC‑GTAGCGCATGGGG‑ACGTGAAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCCCAGACTCTGCTGAATCTg                                                   <  2:132175/247‑1 (MQ=255)                                                                                                                                                                                                     gggg‑CGAACGGG‑G‑CTTGTTTA‑AAGGGGG‑CCAGGACC‑GGAGCGCAGGGGA‑ACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAaacaa                                          <  1:530127/249‑1 (MQ=255)                                                                                                                                                                                                            aaCTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGAc                                    >  1:74155/1‑251 (MQ=255)                                                                                                                                                                                                                 g‑g‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCCTTGCACAGCCTCTGCTGAATCTGGCAAAACAACAAGACGccat                               >  2:467591/1‑248 (MQ=255)                                                                                                                                                                                                                           tG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACAGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTATGATGCGATTGTCTGCACAGATGGCGCGACGTCTGCAAGTCACTTCAGAGAAAGGGGGCAACCTGGCGATCCTCGACGTGACGGGCCGCATGGCACAGACTCTGCTGAATAGGGCAAAACAACCAGACGCTAAGACTCAcc                       >  2:167313/1‑251 (MQ=255)                                                                                                                                                                                                                                    g‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGt               >  1:315431/1‑251 (MQ=255)                                                                                                                                                                                                                                        caggaagggggcgcATGGGT‑ACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTAt             <  1:173397/241‑1 (MQ=255)                                                                                                                                                                                                                                           gaatagtagcgaggggagctgTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGAGGGCGCGTCGGCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGTGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTAt             <  2:292386/229‑1 (MQ=255)                                                                                                                                                                                                                                           aaggtgggcgtaTGGGA‑ACGTGAGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTGCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCa          <  2:68732/240‑1 (MQ=255)                                                                                                                                                                                                                                                   cgcagggtt‑ACGTGCAAAACCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCaaaa  <  1:467591/243‑1 (MQ=255)                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            GCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATT‑GG‑CGAACTGG‑G‑CCTGTTTG‑AAGAGGG‑CCAGGAAC‑GTAGCGCATGGGT‑ACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAA  >  NC_000913/3486137‑3486621

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