Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	70,581	C→T	V65V (GTC→GTT)	araC →	DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	70,581	0	C	T	92.1%	106.2 / 1.0	39	V65V (GTC→GTT)	araC	DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (24/11);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

GATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGA  >  NC_000913/70335‑70817                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              gATGCAATATGGACAATTGGTTTCTTCTATGAATGGCGGGAGTATGAAAAGTATGGCTGCAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCCACGGTTATCTCCGTTTTTTTATCGACCGACCCCTGGGGATGAAAGGTTTTATTCT‑CAATCTCACCCTTCGCGGTCAGGGGGGGGGGAAAAATTAGGGACGAGAATTTGTTTGcc                                                                                                                                                                                                                                           >  1:136549/1‑251 (MQ=255)                     tttCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGAAGCGCCAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAAGCCGATTGAAGCCAAAGGTTATCTCGGTTTTTTTTTCGAACGACCGCTGGGAATGAAAGATTATATTCTTCAAT‑TTACCAATTGCGGGCATGGGGTGGTGAAAAATCCGGGGGGAGAATTTTTTTTCCGACCGGGGGATATTTTGCTg                                                                                                                                                                                                                       >  1:191411/1‑251 (MQ=255)                             ctGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc                                                                                                                                                                                                                  <  2:604919/248‑1 (MQ=255)                               gAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                                                                                                                                                                                               <  1:11023/249‑1 (MQ=255)                                                                          ccccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTGAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATg                                                                                                                                                                  <  1:163747/251‑1 (MQ=255)                                                                                         tACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGGTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCCGGGGGTGGTGAAAAATCAGGGCCGCGAATTTTTTTGCCCACCGGGTGATATTTTACTGTTCCCGCCAGGAGAGTTTACTAACTACGGTCGTCCTCCGGAGGCTCGCGAATGGTAtacccagggggt                                                                                                                                                   >  2:613841/1‑240 (MQ=255)                                                                                          actagtttaatgctgctcgggttgCGGGTTTAACGCCGATTGAGTCCAACGGTTATCTCGATTTTTTTACCGACCGCCCGCTGGGAATGAAAGGTTATAGTCT‑CAATCTCACCATTCGCGTACAGGTGGTGGTGAAAACTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGt                                                                                                                                                   <  1:109031/227‑1 (MQ=255)                                                                                                      cccATCTGGGGGCGGGTTTAACGCCGATTGAGGCCAACGGTTTTCTCGTTTTTTTTATCCCCCGACCTCTGGGAATGAAAGGTTATTTTCT‑CAATCCCACCATTCCCGGTCCGGGGGTGGTGAAAAATCGGGGACGAAAATTTGTTTGCCTACCGGGTGATTTTTTGTTTTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCcg                                                                                                                                      <  1:24774/251‑1 (MQ=255)                                                                                                                ggCGGGTTTAACGCCGATTGAGGCAAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACt                                                                                                                            <  1:463968/251‑1 (MQ=255)                                                                                                                   gggTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGc                                                                                                                         >  1:174484/1‑251 (MQ=255)                                                                                                                    ggTTTAACGAGGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTAT‑CAATCTAACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTCCGGTCGTCATCCGGAGGCTCGCGACTGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGc                                                                                                                         >  2:469344/1‑250 (MQ=255)                                                                                                                         aaCGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAAt                                                                                                                   >  1:447482/1‑251 (MQ=255)                                                                                                                                aTTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTa                                                                                                             >  1:91434/1‑250 (MQ=255)                                                                                                                                    aGGCCAAGGGTTATCTCGATATTTTTATCGTCCTACCGCGGGGAGTGAAGGGTTATATTCA‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGTGACGAGAATTTGTGTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCTTCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑Tgg                                                                                                       <  1:464318/251‑1 (MQ=255)                                                                                                                                                         tttttATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACg                                                                                  >  1:152960/1‑251 (MQ=255)                                                                                                                                                              aTCGACCGACAGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTTTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACggggtt                                                                             <  1:692991/251‑4 (MQ=255)                                                                                                                                                                 gaccgaccGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCCGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGCATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATAGCTTAAC‑TGACCGTCAATATTTGCCAATCCggtgttctt                                                                          >  2:479319/1‑244 (MQ=255)                                                                                                                                                                     gaccgCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGc                                                                      >  1:17939/1‑251 (MQ=255)                                                                                                                                                                      cccgCTGGGTAGGAAAGGTTATTTTCT‑CAATCTCACCTTGCCGGGTCAGGGGGGGGTAAAAAATCAGGGACGATAATTTGTTTACCGACCGGGTGATATTTTGCTTTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGcc                                                                     <  2:476797/250‑1 (MQ=255)                                                                                                                                                                         gCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTTAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACAAGTGGGTTTACTTTCGTCCGCGCGCATACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCgg                                                                  >  1:356394/1‑251 (MQ=255)                                                                                                                                                                         gCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGCATGGCTTAAC‑TGGCCGTCAATATTTGCCCATACgg                                                                                 >  2:402967/1‑236 (MQ=255)                                                                                                                                                                         gCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACgg                                                                                 <  1:402967/236‑1 (MQ=255)                                                                                                                                                                           tGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGAt                                                                >  1:9993/1‑251 (MQ=255)                                                                                                                                                                               aaTGAAGGGTTATATTCT‑CAATTCCACCATTCGCGGTCAGGGGTTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAg                                                            <  1:469344/251‑1 (MQ=255)                                                                                                                                                                               aaTGAAAGGTTATATACA‑CAAGCTCACCATTCGCGGTCAGGGGGTGGTAAAAAAACAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCTCTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAg                                                            <  1:431184/251‑1 (MQ=255)                                                                                                                                                                                            atTCT‑CAATCTCACCATTCGCGGTCAGGGGGAGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGc                                                >  1:295647/1‑250 (MQ=255)                                                                                                                                                                                             ttCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTACCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGATTAAC‑TGGCAGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAgcgccaccgccgcct                                              >  1:50450/1‑240 (MQ=255)                                                                                                                                                                                              tct‑cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATg                                                                                                                  >  1:462089/1‑183 (MQ=255)                                                                                                                                                                                              tct‑cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGCTGAAGCGCACCAGCCGCAtt                                             >  1:182478/1‑251 (MQ=255)                                                                                                                                                                                               ct‑cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCCCCAGCCGCAttt                                            >  1:202501/1‑251 (MQ=255)                                                                                                                                                                                                     tatcACCATACGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGACTTTGTTTGCCGAACGGATGATATTTTACTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATTGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCATCTGGAATGAAGGTCTTAACGTGGC‑GTCAATTTTTGACAAAACGGGGCTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTAAGCg                                       >  2:533456/3‑251 (MQ=255)                                                                                                                                                                                                        cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAcc                                    >  1:574027/1‑251 (MQ=255)                                                                                                                                                                                                             ttCGCGGTCAGGGGGTGGTGAAAATTCAGGGACGAGATTTTTTTGGCCTACGGGGTGATATTTTGTTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                                <  1:320954/250‑1 (MQ=255)                                                                                                                                                                                                                  ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGc                           >  1:29091/1‑250 (MQ=255)                                                                                                                                                                                                                  ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGc                           >  2:705362/1‑250 (MQ=255)                                                                                                                                                                                                                      aGGGGGTGGTGAACAATCAGGGACGAGAATTTTTTTGCCGACCGGGTGATATTTTGATGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGGATCACCAGTGGGTTTACTTTCGACCGCGCGCCTACTGGCATGCATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAACGCCCCCACCGCATTTCACCGCCCTGTTTTGGCCAATc                      >  2:580660/1‑251 (MQ=255)                                                                                                                                                                                                                             ggtgAAAAATAAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCCTCAC‑TGGCCGTCACTATTTGCCCATACGGGGTTTTTTCGCACGGATGAAGCGCCCCAGCCGACTTTCAGCGACCTGTTTGGGCAAATCATTAACg               >  2:317415/1‑251 (MQ=255)                                                                                                                                                                                                                                        aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCCGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAggg    >  2:601780/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTACCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGCGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCAGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGa  >  2:415032/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              GATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCT‑CAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAAC‑TGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGA  >  NC_000913/70335‑70817

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