Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,170	T→C	S262P (TCG→CCG)	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	71,170	0	T	C	95.3%	119.6 / ‑1.9	43	S262P (TCG→CCG)	araC	DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (22/19);  total (23/20)
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.80e-01

CTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTG  >  NC_000913/70931‑71410                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  cTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCATATCGCCACCGTCGTGCGCAATGTTGGTTTTGCCGATCCACTCTATTTCCCGCCAGTAttt                                                                                                                                                                                                                                        >  2:123737/1‑251 (MQ=255)     tCAGTACATCAGCGAACACTTGGCAGACAGCAATTTTGATATCGCCAGCGCCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACAGCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACTCATTAGTCAGGCGAAGCTGCTTTTGATCTCTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTaaa                                                                                                                                                                                                                                     <  2:267394/250‑1 (MQ=255)           cATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGc                                                                                                                                                                                                                               >  2:461753/1‑250 (MQ=255)                            cAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTACGCTGGCGCGAGGACCACCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGCTGCCTATCGCAACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTACAAAATGCACCGGGGCCAGCCcgcgc                                                                                                                                                                                                             >  2:203880/1‑248 (MQ=255)                             aGACAGCAATGTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGGCGCCGGCGCGTCTGGCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGTTTTTGACGATCAACTGTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                                                                                                                                                                                            <  1:214153/251‑1 (MQ=255)                                      ttttGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGTGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGTTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGGTGGTATTGACGATGAACTCTATTTACCGAGAGTATTTAAAACATGCACCGGGGACAGCCCGAGCGAGGTTCg                                                                                                                                                                                                     >  2:57913/1‑249 (MQ=255)                                               cGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                                                                                                                                                                                          <  1:454656/251‑1 (MQ=255)                                                      gTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTCAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGCTCAACTCTATTTCCCGCGAGTATTTAAAAACTGCACCGGGGCCAGCCCGAACGACTTTCGTGCCGGTTGTGAAGaaa                                                                                                                                                                                    >  1:453452/1‑250 (MQ=255)                                                         gCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTAAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCAACGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAgg                                                                                                                                                                                >  2:493762/1‑250 (MQ=255)                                                           acaGCATGTTTGCTGGTCGCCGTCGCGTCTGTCACATCGTTGCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCTCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGa                                                                                                                                                                              <  2:221937/251‑1 (MQ=255)                                                                 ggTTTGCTTGTCGCCGGCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGGCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCTCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATGTAAAAAATGCACCGGTGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGGGAatgatg                                                                                                                                                                        <  2:487076/250‑1 (MQ=255)                                                                   ttttcTTGTCGCCGTCGCGTCTGTCGCATCTTGTGCGCGAGCAGTTAGAGATGAGCGTCTTAAGCTTGCGCGAAGACCAACGCATTAGTCAGGCGAAGCTGCTTTTTAGCACTTCCCGGATTCCTACCGCCACCGTCGGTCGCAATGTAGGTTTTGACGATCAACTCGATTTCCCGCGAGTTGTTAAAAAATGATCCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAATAAATAGCCAATGATGt                                                                                                                                                                       <  2:24428/246‑1 (MQ=255)                                                                         tGTCGCCGTTGCGTTTTTCACATCTTTTCCGCCAGCAGTTCGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAATCTTCTTTTGAGCATTACCCGGATGCCTATCGTCACCGGCGGTCGCAATGTGGGTTTTGACGATCAACTCTATTTCCCGCGAGCATTTAAAAAATGCACTGGGGCCATCCCGAGCGAGTGTCGTGCCGGTTCTGAAGAAAAAGTGAATCAGGTAGCCGt                                                                                                                                                                 <  2:299327/250‑1 (MQ=255)                                                                         tGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                                                                                                                                                                >  1:261090/1‑251 (MQ=255)                                                                                      tCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAAt                                                                                                                                                  >  1:12027/1‑251 (MQ=255)                                                                                       cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGATTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACAGGGGCCAGCCCGAGCTAGGTTCGTGACGGTTGTGACGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAAtt                                                                                                                                                 >  2:123836/1‑251 (MQ=255)                                                                                           cacaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGAAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGt                                                                                                                                              >  1:128463/1‑250 (MQ=255)                                                                                                          agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAAtt                                                                                                                              >  1:11595/1‑251 (MQ=255)                                                                                                          agcagTTAGGGATTAGCGTCTGAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAAtt                                                                                                                              <  2:233174/251‑1 (MQ=255)                                                                                                           gcagTTAGGGATGAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACGACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCAGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTg                                                                                                                             <  1:118155/251‑1 (MQ=255)                                                                                                                    gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTg                                                                                                                     >  1:207665/1‑250 (MQ=255)                                                                                                                           gacTTAAGCTGGCTCGAGGACCAACGCTTTAGAAATGCGAAGCTGCTTCTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTGTTGACGATCAACTCTATTTCCCGCGAGTATTTACAAAATGCACCGGGGCCAGCCCGAGCGGGGTTCGTGCCGGTTGTGAGGAAAATGAGTATGATGTAGCCTTAAAGT‑GTGCATAATTGGGAAAGAATCAAACAATTGCCGGCtttgcggaat                                                                                                              >  2:143430/3‑242 (MQ=255)                                                                                                                                    tGGCGCGAGGGCCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGCTCAACTCTATTTCCCGCGAGTATTTAAAACATGCACCGGGGACAGCCCGAGCGAGTTTCGTGCCGGGTGAGAAGACAAAGGGAAGGCGTTAGCCGACAAGATGG‑CCTCATTGGAAACTAATCAGAATATTGACAGCTTGACGGAGTAGcctacggt                                                                                                    >  2:439668/1‑244 (MQ=255)                                                                                                                                       cgcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTg                                                                                                 >  2:137257/1‑251 (MQ=255)                                                                                                                                                 aaCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCt                                                                                       <  1:205529/251‑1 (MQ=255)                                                                                                                                                 aaCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCt                                                                                       >  1:136931/1‑251 (MQ=255)                                                                                                                                                                cGAAGCTGCTTTTGAGCACTCCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGa                                                                        <  2:308492/251‑1 (MQ=255)                                                                                                                                                                 gAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGGCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTTACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGAc                                                                       <  2:402591/251‑1 (MQ=255)                                                                                                                                                                  aaGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACa                                                                      <  1:137257/251‑1 (MQ=255)                                                                                                                                                                    gctgctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAgg                                                                    >  1:964/1‑251 (MQ=255)                                                                                                                                                                    gctgctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGGGAATGATGTAGCCGGCAGGTTGT‑CATAATTGGTAACGAATCAGACAAATGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTAGTCAATTTGAccgg                                                                    >  2:43859/1‑248 (MQ=255)                                                                                                                                                                        ctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGcac                                                                 >  2:423231/1‑250 (MQ=255)                                                                                                                                                                            tGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTa                                                             <  1:423231/250‑1 (MQ=255)                                                                                                                                                                                        ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCt                                                <  2:449931/251‑1 (MQ=255)                                                                                                                                                                                             cctacctcccccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTAGTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAac                                           <  2:401268/241‑1 (MQ=255)                                                                                                                                                                                                  cGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACtt                                       <  1:461753/250‑1 (MQ=255)                                                                                                                                                                                                     cACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTttatt                                   <  1:123737/251‑1 (MQ=255)                                                                                                                                                                                                      cccGTCGGTCGCAATGTTGGTTTTGACGATCAACGCTATTTCCCGCGAGTATTTAAAAAATGCACCGGTGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTttatta                                  <  1:364779/250‑1 (MQ=255)                                                                                                                                                                                                      aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGCGCGAGTTTCGTGCCGGGTGTGACGACAAAGTGAGTGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAGTCAGACAATTGACGGCTTGACGTAGTAGCCTAGGGTTTGCAGAGTCCCTGCTTCGGCCATTTGGCAGGCCCATTATGCAATCATTGCTGGCACACTttatt                                   >  1:324115/1‑250 (MQ=255)                                                                                                                                                                                                         gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACttt                                      <  1:163907/244‑1 (MQ=255)                                                                                                                                                                                                         gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGGTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGaccccttt                                      >  2:163907/1‑237 (MQ=255)                                                                                                                                                                                                                        ggTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTAtt                >  1:404894/1‑251 (MQ=255)                                                                                                                                                                                                                                      aCTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCggtg  >  1:352834/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  CTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGT‑CATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTG  >  NC_000913/70931‑71410

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