Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,875	T→C	G175G (GGT→GGC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,875	0	T	C	96.5%	151.8 / ‑3.7	60	G175G (GGT→GGC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (10/46);  total (13/47)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

TACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGG  >  NC_000913/71640‑72116                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  tacccgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGtt                                                                                                                                                                                                                                     <  1:920128/245‑1 (MQ=255)  actcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACACCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTtgc                                                                                                                                                                                                                                   <  1:505415/247‑1 (MQ=255)    tcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCCGCCGGAATGCCGAGCGGCGGGTTTAAATgggttctg                                                                                                                                                                                                                                 >  1:677428/3‑245 (MQ=255)       aTAAACCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                                                                                                                              <  2:336791/251‑1 (MQ=255)       aTAAAACTGAACATGCGTTGCATCAACACAGCATTTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                                                                                                                              <  1:191263/251‑1 (MQ=255)        aaaaactaaaccggcgctGCATCCACACTGCAGGTGCACCATTCGGGGCGGGCGCTTTGTTGGCCCGACGCTGCCGCGGGTGCCAATGTGGGCGGGAAGGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTgg                                                                                                                                                                                                                             <  2:388327/234‑1 (MQ=255)           aaCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCCGCCGGAATGCAGAGCGGGGGGTTTAAATGGTTGCTGCTGGcgg                                                                                                                                                                                                                          >  1:708062/1‑249 (MQ=255)                aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGGGGGCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCCGCCGGAATGAAGAGCGGCGGGGTTAAATGGTTGGTGCTGGcggaggcg                                                                                                                                                                                                                     >  1:142919/1‑244 (MQ=255)                   aTGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTAGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGATGGGTATGCCGGTGGCAAAATTAATTAAGCCGAATATTATCGGCTGACTGGTGTGGACACCGTTTTAATTCCCGCCTGGGGTTCTGGGGGGGCGCGG‑GGTCGATATTCCAGGCGGCATCCAGAGCGGGGGGTTTTAATGGTtgctggcggggacggaggc                                                                                                                                                                                                                   >  1:415293/1‑236 (MQ=255)                     gctttGCATCAACCCAGCAGTTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtat                                                                                                                                                                                                                 <  2:114978/247‑3 (MQ=255)                     gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGAGATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtatt                                                                                                                                                                                                                <  1:637119/251‑4 (MQ=255)                         tGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGCATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATTTTCCAGCCGGAATGCAGAGCGGCGAGTTTTAATGGTTGCTGCTggggacggcggtatttgtg                                                                                                                                                                                                            >  1:81655/1‑234 (MQ=255)                                    aGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGGATGCAAAGCGGCGAGTTTAAATGGTTTCTTACGGCGACGGCGGTATTTGTGTGGGTttgtgg                                                                                                                                                                                                 >  1:611649/1‑251 (MQ=255)                                          ttCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGGTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTTGTggctggct                                                                                                                                                                                           >  1:269377/1‑251 (MQ=255)                                                      gacggtcgCTTTGTGGGCCCGCCGCGTCCGCGGGTGCCAATGGTGGCGGGAATGGTGGATCTGCCGGTGGCAAAATTCATTGCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCGGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGtt                                                                                                                                                                               <  2:203310/249‑1 (MQ=255)                                                        cggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGGTCGATATTCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGGCGGCGGGATTTGTGTGGGGTGGTGGCTGGCTGTGCTGGCgggttt                                                                                                                                                                             >  1:603635/1‑247 (MQ=255)                                                         ggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTGGGGGGTGGCTGTGCTGGCgggtgtt                                                                                                                                                                            >  1:230587/1‑246 (MQ=255)                                                               tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCCGCCGGGATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTTGGTGGGGTGGGGGCTGGCTGtggtgggggggtttgggcgg                                                                                                                                                                      >  1:117252/1‑237 (MQ=255)                                                               tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGCGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTTGTGGCTGGCTGTGCTGGCgggtttggggggg                                                                                                                                                                      >  1:92557/1‑240 (MQ=255)                                                                 tggtgtcccggcgggtcctcgggTGCAAATGGTTGCGGGATTGTTGGATCTGCCGGTGGCAAAATGCATTACGCCGAATATTATCGTCTGCCGGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCTGCGATTGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGTTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                                                                                                                                     <  2:64203/229‑1 (MQ=255)                                                                 ggtttgcccaccgcgTCCGCTGGTGCCAATGGTGGCGGGAAGGCTGGCTCTGCCGGTGTCCAAATTCACTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGGTTGACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGGTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                                                                                                                                     <  2:394733/239‑1 (MQ=255)                                                                      gCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTTCCGCCGGGATGCAGAGCGGCGAGGTTAAATGGGTGCTGGTGGCGGCGGCGGTTTTTTTGTGGGGTGGGGGCTGGCGGTGCTTGCgggttttgtggtccgggg                                                                                                                                                                 >  2:307855/1‑233 (MQ=255)                                                                             gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                                                                                                                                        <  1:473692/251‑6 (MQ=255)                                                                             gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                                                                                                                                        <  1:149110/251‑6 (MQ=255)                                                                                cccGCTGGTGCCACTGGTGGCGGGAATGCTGGAGCTGCCGGTGGCAAAAGTCATTTCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccgatc                                                                                                                                                     <  2:427747/250‑9 (MQ=255)                                                                                               ggtggCGGGAATGCTGGATCTGCCGGTGGCAAAAGTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTATTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCattat                                                                                                                                       <  1:993768/250‑1 (MQ=255)                                                                                                     gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                                                                                                                                    >  1:118778/1‑167 (MQ=255)                                                                                                     gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                                                                                                                                    <  2:118778/167‑1 (MQ=255)                                                                                                        aaTGCTGGAGCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcgcg                                                                                                                             <  1:865995/251‑1 (MQ=255)                                                                                                         aTGCTGGATCTGCAGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCgg                                                                                                                            <  1:783426/251‑1 (MQ=255)                                                                                                               gttctgtgtgtgtcaaaagtcattaCGGCGAATATTATCGGCTGCCGGCTGTGGTATTCGTTTTTCTTCATGACTGCGATTCTGGCGGG‑CGCGTCGATCGATGGTCCAGCCGGAATTCAGAGTGGCGAGTTGAAATGTTTGCTGATGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGAGGTAATGCAACCGATCGTTTAAGTCATGATTTGTCGCGCGGCCGttt                                                                                                                      <  2:756836/230‑1 (MQ=255)                                                                                                                    gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTTGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGTCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttg                                                                                                                     <  1:242311/247‑1 (MQ=255)                                                                                                                    gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttgt                                                                                                                    <  1:137733/248‑1 (MQ=255)                                                                                                                     ccGGTGGCAAAATTCATTACGCCGAATATTTTCGGCTGCCTGCTGTGGCCACCGTTTTTCTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTtgtg                                                                                                                <  2:648349/251‑1 (MQ=255)                                                                                                                      cGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTgg                                                                                                               <  1:4234/251‑1 (MQ=255)                                                                                                                          ggcaaaaTTCATTACGACGAAGAGTATCGGCTGGCTGCTGTGGCCATCGTTTTACTGCCGGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTg                                                                                                            <  2:183427/246‑1 (MQ=255)                                                                                                                                           cGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGcc                                                                                          <  1:553288/251‑1 (MQ=255)                                                                                                                                               gattagCGGCTGCCTGCTGGGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAGTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACTGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTAGGGCGGAGTGGTAAAGCAACCGATCGTTTAAGTCATTATTTATCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCg                                                                                      <  1:616216/250‑1 (MQ=255)                                                                                                                                                 ttaccGGCTGCCTTCTGTGGCCACCTTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATTGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCTATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCgg                                                                                     <  2:830481/246‑1 (MQ=255)                                                                                                                                                  tatCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGc                                                                                    <  1:73009/250‑1 (MQ=255)                                                                                                                                                          ggctgtggtgGCCTCGGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgcg                                                                           <  2:81655/244‑1 (MQ=255)                                                                                                                                                          gcctgcTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgc                                                                            <  2:344918/250‑1 (MQ=255)                                                                                                                                                          gcctgcTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCTGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGtggtgg                                                                             <  1:924854/249‑1 (MQ=255)                                                                                                                                                                   ggCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGgtggtg                                                                  <  1:841082/251‑1 (MQ=255)                                                                                                                                                                               aCTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCa                                                       <  2:708062/250‑1 (MQ=255)                                                                                                                                                                                   cctgcctgGGATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt                                                  <  1:663762/251‑1 (MQ=255)                                                                                                                                                                                      gcctgggATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGtt                                                 <  2:937636/245‑1 (MQ=255)                                                                                                                                                                                      gcctgggATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCAGCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTg                                                <  2:891673/246‑1 (MQ=255)                                                                                                                                                                                          gggATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCATAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGGGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGcc                                           <  2:799652/251‑1 (MQ=255)                                                                                                                                                                                           ggATTGTGGCGGG‑CGCGTCTATCGATATTCCAGCCGGAATGGAGAGCGGCGAGTTTAAATGGTTGCTGCGGGCGTCGGCGGTATGTGTGTGGTGTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGTCGTTGATTTCTGCCATCGGCGTGGTGGCTCTGGTGGTGTTAAGTCGCCCCCCGTTGATGCCg                                          <  2:533777/251‑1 (MQ=255)                                                                                                                                                                                            gATTCTGGCGGG‑CGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCgg                                         <  1:334793/251‑1 (MQ=255)                                                                                                                                                                                                                    atatTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaaa                  <  1:734050/251‑1 (MQ=255)                                                                                                                                                                                                                      atTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAg                 <  1:356411/250‑1 (MQ=255)                                                                                                                                                                                                                        tccagCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTg               <  2:400324/246‑1 (MQ=255)                                                                                                                                                                                                                         ccagCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTgg              <  2:24512/247‑1 (MQ=255)                                                                                                                                                                                                                         ccagCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTgg              <  1:555721/247‑1 (MQ=255)                                                                                                                                                                                                                              cGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCGGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTgggg        <  2:232588/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaTGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGttt     <  1:266385/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaTGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGttt     <  1:544371/251‑1 (MQ=255)                                                                                                                                                                                                                                    gCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAgg  <  1:235711/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  TACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGG  >  NC_000913/71640‑72116

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