Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,770	T→C	F140F (TTT→TTC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,770	0	T	C	95.8%	142.0 / ‑3.8	49	F140F (TTT→TTC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/27);  total (20/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.56e-01

TGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTT  >  NC_000913/71532‑71999                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      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)                                                                                                                                                                                                                           ggagctgccgttgtcaaaaTTCATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCACCGTTTTACTGCCTGCCTGGGATTCGGGCGGGAGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGtt  <  2:104046/235‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      TGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCT‑GCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTT  >  NC_000913/71532‑71999

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: ‑