Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,764	T→A	A138A (GCT→GCA)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,764	0	T	A	95.8%	144.4 / ‑2.9	50	A138A (GCT→GCA)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (18/28);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTC  >  NC_000913/71520‑71996                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         ggCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAAttc                                                                                                                                                                                                                                    <  2:126273/251‑2 (MQ=255)  gCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTAGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAAttc                                                                                                                                                                                                                                    >  1:267865/1‑249 (MQ=255)             tggcagggattgtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCT‑TCCGGTGGCAAAATTCATTACGCCGAat                                                                                                                                                                                                                        >  1:177059/25‑251 (MQ=255)                       ggtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACGCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCGGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGc                                                                                                                                                                                                               <  1:319090/236‑1 (MQ=255)                                      gggCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgccac                                                                                                                                                                                               >  1:248788/1‑249 (MQ=255)                                       ggcgactggttctccgtatggctggGCTGGCGTTTTAAAAAGCCGTTGCATCGGTGGTCATTGCGGAAGAAAAACAAAGCATTACTCGATAAATCTGAACATGCGTTGCATCAACACAGCATGGTCACCATTCTGGTCGGTCGCTTGGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcaacc                                                                                                                                                                                              <  2:174077/233‑5 (MQ=255)                                       ggCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgccacc                                                                                                                                                                                              <  1:432338/251‑4 (MQ=255)                                          gactgaaactccttctGGATGGGCTTGCGTTTTAAAAAGCCGCTGCATCGCAGGTCATTTCTGAAGAAAATCAAAGCATTACTAGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACTATTCTGGTCGGTCTCTTTGTTGGCACGACGCGTACGCTGGTGCCAATGGTGGCGGGGATGCTGGATCTTGCCGG‑GGCGAAATTCATTACCCAGAATATTTTCGGCGGCctgctgctgccaccctt                                                                                                                                                                                           >  2:307781/13‑240 (MQ=255)                                              ggatctccttctGGCTGGTGTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTCCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGAATCAACACAGCATGTTCACCATCCTGGTCGGTCGCTTTGTTGGCCCGATGCGGCCGCTGGTGCCAATGGTGGTGGGAATGCTGGATCT‑TCCGGTTGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTa                                                                                                                                                                                        <  1:408613/242‑1 (MQ=255)                                              ggatctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGtttt                                                                                                                                                                                         <  1:384339/241‑1 (MQ=255)                                                                gctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTGCTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCg                                                                                                                                                                                                                 <  1:31735/205‑1 (MQ=255)                                                                  tGGCGTTTTAAAAAGCCGATGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGAGTCCGCTGGTGCCAATGGTGGCGGGGATGGTGGATCT‑GCCGGTGGCAAAATTCAATACGCCGTCTGTTTTCGGCTGTCTGCTGTGGGCACTGTTTTACTTGCTgactggggattgggg                                                                                                                                                                   >  2:306171/1‑237 (MQ=255)                                                                      gTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcgggg                                                                                                                                                               >  1:543892/1‑250 (MQ=255)                                                                        tttAAAAAGCCGTTGCATCGATGGTCATTTCTGAAGAAAAACAAAACATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑TCCGGTGGCAAAATTCATTACGCCGAATATTAACGGCTGCCTGCTGTGGACACCGTTTTTCTTCCTGCCTGGGCATCTggcgggggg                                                                                                                                                             >  2:396293/1‑248 (MQ=255)                                                                         ttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGGTTCTGGCGGGcgcg                                                                                                                                                            >  1:326382/1‑251 (MQ=255)                                                                          tAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGcgcg                                                                                                                                                            <  1:396293/250‑1 (MQ=255)                                                                          tAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GACGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGACTGGGATTCTGGCGGGCgcgg                                                                                                                                                           >  1:396221/1‑251 (MQ=255)                                                                                        aTCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTcc                                                                                                                                              >  1:321187/1‑250 (MQ=255)                                                                                            cTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGctgcct                                                                                                                                                                                                          <  1:52392/186‑1 (MQ=255)                                                                                            cTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCACGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATTCTGGATCT‑GCCGGTGGCAAAATTCATTACACCGAATATTATCGGctgcct                                                                                                                                                                                                          >  2:52392/1‑186 (MQ=255)                                                                                                 aagttCTGAAGAAAACCAAAGCATTACTCAGTAAATCGGAACAGCCGTGTCATCAACATAGCATGTCCACCAGTCTGTGCGGTCGTTTGGTTGTGCCGAGGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCG‑GCCGGTGGCAAAATTCCTGACGCCGACTATTATCGCCTGCCTGCTGTGGCCACCGGTTTACTTCCTGCCTTGGATCCGGGCGGGCGCGTCGATCGATATTCCAGCCGGaag                                                                                                                                     <  2:396719/247‑2 (MQ=255)                                                                                                  aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTCCAGCCGGAATGc                                                                                                                                   >  1:684186/1‑251 (MQ=255)                                                                                                      cTGAAGAAAAACAAAGCATTACTCGATAAACCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGACTGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATCTTATCGGCTGCCTGCTGTGGCAACCGTTTTAATTCCTGCCTGGGATTTGGGCGGGGGCGGCGAGCCATATTACCGCCGGGATGCag                                                                                                                                 >  2:153498/1‑249 (MQ=255)                                                                                                       tGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGAATGCAGAGc                                                                                                                              >  1:13277/1‑251 (MQ=255)                                                                                                        gaagaaAAACAAAGCATTACTCGATAACACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGGGGCAAAACTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGATTTACTTCCTGCCTGGGATTCTGGCGGGGGCGGCGAGCGGTATTCCAGCCGGAATGCAGAGCg                                                                                                                             >  1:276772/1‑251 (MQ=255)                                                                                                        aaagaaAAACAAAGGATTACGCTATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCGGGTCGGTCGCTTTGTTGGCCCGAGGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATAGTCCAGCCGGAATGCAGAGCg                                                                                                                             <  2:686909/250‑1 (MQ=255)                                                                                                          agcaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGATGCTGGATCT‑GCCGGTGGCAAAATTCATTAAGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGGGCGGCGAGCGATATTcaatccgggatgaagagcggg                                                                                                                           >  2:120214/4‑231 (MQ=255)                                                                                                                                 aaaaCTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTgg                                                                                                     <  1:20845/250‑1 (MQ=255)                                                                                                                                  aaaCTGAACTTGCGTTGCAGCAACCCAGCATGTGCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGTGACCGCTGGTTCTAATGGTGGCGGGAAGGCTGGATCT‑GCCGGTGGCTAAATTTATTTCGCCATATATTATCGGCTGCTTGCTGTGGCCGCCGGTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATGGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCa                                                                                                   <  1:102577/251‑1 (MQ=255)                                                                                                                                    gttGAACAGGCGTTGCACCAACCCAGCATGTTCCCCCTTCTGGGCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGTCAATGGTGGCGGGAATGTTGGTTCT‑GCCGGTGGCAAAATGCATTACGCTGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGcgac                                                                                                 <  2:224300/249‑4 (MQ=255)                                                                                                                                      tGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCCAAATTCATTACGCCGGATATTATCGGCTGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTCCCGCCGGAATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGcgacgg                                                                                               >  1:38484/1‑246 (MQ=255)                                                                                                                                       gAACATGCGTTGCAGCAACACAGCATGTTCGCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAGGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGcgacggc                                                                                              <  2:353706/251‑7 (MQ=255)                                                                                                                                          cATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGGGGTTAAATGGGTGCTGCTGGCGAcggcggg                                                                                           >  1:371943/1‑250 (MQ=255)                                                                                                                                           aTGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGGTATTCCAGCCGGGAAGCCGAGCGGCGAGTTTTAATGGTTGCTGCTGGcggcggcgggt                                                                                          >  1:670979/1‑241 (MQ=255)                                                                                                                                                    aTCAACCCAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTgg                                                                                 <  1:320239/251‑1 (MQ=255)                                                                                                                                                                 gttccctagttgggtcggtcGTTTTTTTGGCCCGACTCGTCCGCTGGGGCCAATGGTGGCGGGAATGCGGGATCT‑GTCGGTGGCAAAATTCATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCGGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGtggctggc                                                                    <  2:276772/239‑1 (MQ=255)                                                                                                                                                                       cATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGGCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCt                                                              <  1:32503/251‑1 (MQ=255)                                                                                                                                                                                 ggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATg                                                    <  1:360386/251‑1 (MQ=255)                                                                                                                                                                                               cccGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGGTGGCGACGGCGGGAATTGTGTGGGGTGGTGGGTGGCTGTGCTGGCGGtttttggggagggggtaa                                       >  1:384496/1‑234 (MQ=255)                                                                                                                                                                                                 cGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTtggt                                                                           <  1:143259/212‑1 (MQ=255)                                                                                                                                                                                                   gcgcgtccgtgggggtcAATGGGGGGGGGAATGTTGGAGCT‑GCCGGTGGAAAAATTAATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTGAAATGGTTGCTGCTGGCGACGGCGGTATGTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaac                                  <  2:404764/235‑4 (MQ=255)                                                                                                                                                                                                     gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGGGCGGCGGTCGATATTCCTGCCGGGATGCAGAGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGGTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCgggtttggggcgcgggggaaacggc                                  >  1:287141/1‑226 (MQ=255)                                                                                                                                                                                                               gggacattggtggCGGGACTGCTGGATCT‑GCCGGTGGCAAAATTCATTTCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTGCTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATGTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGtttaag                      <  2:686860/244‑4 (MQ=255)                                                                                                                                                                                                                      tggtggCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgc                                                  <  2:186724/216‑1 (MQ=255)                                                                                                                                                                                                                        ggggggggAATGTTGGATCT‑GCCGGTGGCAAAATTCATTCCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATGCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTAttt             <  2:427604/246‑1 (MQ=255)                                                                                                                                                                                                                        ggggCGGGAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGTTTGCGACGGCGGTATTTGTGTGGGTTGGTGGCGGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTAttt             <  1:307781/249‑1 (MQ=255)                                                                                                                                                                                                                          ggCGGGAATGCTGGATCT‑GGCGGGGGCAAAATTCATTACGCCGAATATTAACGGCTGCCTGCTGTGGCCACCGTTTTGCTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGt           <  2:670979/251‑1 (MQ=255)                                                                                                                                                                                                                             gggAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgc        <  1:168389/251‑3 (MQ=255)                                                                                                                                                                                                                              ggAATGCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcg       <  1:621154/251‑4 (MQ=255)                                                                                                                                                                                                                                   gCTGGATCT‑GCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCggcc  <  1:120214/251‑3 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         GGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTC  >  NC_000913/71520‑71996

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