Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,806	G→A	P152P (CCG→CCA)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,806	0	G	A	91.5%	113.3 / ‑1.1	48	P152P (CCG→CCA)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (14/29);  total (16/32)
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.98e-01

GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGT  >  NC_000913/71565‑72051                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      ggatctccttctGGCTGGTGTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTCCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGAATCAACACAGCATGTTCACCATCCTGGTCGGTCGCTTTGTTGGCCCGATGCGGCCGCTGGTGCCAATGGTGGTGGGAATGCTGGATCTTCCGGTTGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTa                                                                                                                                                                                                                                               <  1:408613/242‑1 (MQ=255) ggatctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGtttt                                                                                                                                                                                                                                                <  1:384339/241‑1 (MQ=255)                     tGGCGTTTTAAAAAGCCGATGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGAGTCCGCTGGTGCCAATGGTGGCGGGGATGGTGGATCTGCCGGTGGCAAAATTCAATACGCCGTCTGTTTTCGGCTGTCTGCTGTGGGCACTGTTTTACTTGCTgactggggattgggg                                                                                                                                                                                                                          >  2:306171/1‑237 (MQ=255)                         gTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcgggg                                                                                                                                                                                                                      >  1:543892/1‑250 (MQ=255)                           tttAAAAAGCCGTTGCATCGATGGTCATTTCTGAAGAAAAACAAAACATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTTCCGGTGGCAAAATTCATTACGCCGAATATTAACGGCTGCCTGCTGTGGACACCGTTTTTCTTCCTGCCTGGGCATCTggcgggggg                                                                                                                                                                                                                    >  2:396293/1‑248 (MQ=255)                            ttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGGTTCTGGCGGGcgcg                                                                                                                                                                                                                   >  1:326382/1‑251 (MQ=255)                             tAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGcgcg                                                                                                                                                                                                                   <  1:396293/250‑1 (MQ=255)                             tAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGACGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGACTGGGATTCTGGCGGGCgcgg                                                                                                                                                                                                                  >  1:396221/1‑251 (MQ=255)                                           aTCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTcc                                                                                                                                                                                                     >  1:321187/1‑250 (MQ=255)                                                    aagttCTGAAGAAAACCAAAGCATTACTCAGTAAATCGGAACAGCCGTGTCATCAACATAGCATGTCCACCAGTCTGTGCGGTCGTTTGGTTGTGCCGAGGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCGGCCGGTGGCAAAATTCCTGACGCCGACTATTATCGCCTGCCTGCTGTGGCCACCGGTTTACTTCCTGCCTTGGATCCGGGCGGGCGCGTCGATCGATATTCCAGCCGGaag                                                                                                                                                                                            <  2:396719/247‑2 (MQ=255)                                                     aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTCCAGCCGGAATGc                                                                                                                                                                                          >  1:684186/1‑251 (MQ=255)                                                         cTGAAGAAAAACAAAGCATTACTCGATAAACCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGACTGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATCTTATCGGCTGCCTGCTGTGGCAACCGTTTTAATTCCTGCCTGGGATTTGGGCGGGGGCGGCGAGCCATATTACCGCCGGGATGCag                                                                                                                                                                                        >  2:153498/1‑249 (MQ=255)                                                          tGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGAATGCAGAGc                                                                                                                                                                                     >  1:13277/1‑251 (MQ=255)                                                           gaagaaAAACAAAGCATTACTCGATAACACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAACTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGATTTACTTCCTGCCTGGGATTCTGGCGGGGGCGGCGAGCGGTATTCCAGCCGGAATGCAGAGCg                                                                                                                                                                                    >  1:276772/1‑251 (MQ=255)                                                           aaagaaAAACAAAGGATTACGCTATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCGGGTCGGTCGCTTTGTTGGCCCGAGGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATAGTCCAGCCGGAATGCAGAGCg                                                                                                                                                                                    <  2:686909/250‑1 (MQ=255)                                                             agcaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGATGCTGGATCTGCCGGTGGCAAAATTCATTAAGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGGGCGGCGAGCGATATTcaatccgggatgaagagcggg                                                                                                                                                                                  >  2:120214/4‑231 (MQ=255)                                                                                    aaaaCTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTgg                                                                                                                                                            <  1:20845/250‑1 (MQ=255)                                                                                     aaaCTGAACTTGCGTTGCAGCAACCCAGCATGTGCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGTGACCGCTGGTTCTAATGGTGGCGGGAAGGCTGGATCTGCCGGTGGCTAAATTTATTTCGCCATATATTATCGGCTGCTTGCTGTGGCCGCCGGTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATGGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCa                                                                                                                                                          <  1:102577/251‑1 (MQ=255)                                                                                       gttGAACAGGCGTTGCACCAACCCAGCATGTTCCCCCTTCTGGGCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGTCAATGGTGGCGGGAATGTTGGTTCTGCCGGTGGCAAAATGCATTACGCTGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGcgac                                                                                                                                                        <  2:224300/249‑4 (MQ=255)                                                                                         tGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCCAAATTCATTACGCCGGATATTATCGGCTGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTCCCGCCGGAATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGcgacgg                                                                                                                                                      >  1:38484/1‑246 (MQ=255)                                                                                          gAACATGCGTTGCAGCAACACAGCATGTTCGCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAGGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGcgacggc                                                                                                                                                     <  2:353706/251‑7 (MQ=255)                                                                                             cATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGGGGTTAAATGGGTGCTGCTGGCGAcggcggg                                                                                                                                                  >  1:371943/1‑250 (MQ=255)                                                                                              aTGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGGTATTCCAGCCGGGAAGCCGAGCGGCGAGTTTTAATGGTTGCTGCTGGcggcggcgggt                                                                                                                                                 >  1:670979/1‑241 (MQ=255)                                                                                                       aTCAACCCAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTgg                                                                                                                                        <  1:320239/251‑1 (MQ=255)                                                                                                                    gttccctagttgggtcggtcGTTTTTTTGGCCCGACTCGTCCGCTGGGGCCAATGGTGGCGGGAATGCGGGATCTGTCGGTGGCAAAATTCATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCGGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGtggctggc                                                                                                                           <  2:276772/239‑1 (MQ=255)                                                                                                                          cATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGGCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCt                                                                                                                     <  1:32503/251‑1 (MQ=255)                                                                                                                                    ggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATg                                                                                                           <  1:360386/251‑1 (MQ=255)                                                                                                                                                  cccGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGGTGGCGACGGCGGGAATTGTGTGGGGTGGTGGGTGGCTGTGCTGGCGGtttttggggagggggtaa                                                                                              >  1:384496/1‑234 (MQ=255)                                                                                                                                                    cGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTtggt                                                                                                                                  <  1:143259/212‑1 (MQ=255)                                                                                                                                                      gcgcgtccgtgggggtcAATGGGGGGGGGAATGTTGGAGCTGCCGGTGGAAAAATTAATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTGAAATGGTTGCTGCTGGCGACGGCGGTATGTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaac                                                                                         <  2:404764/235‑4 (MQ=255)                                                                                                                                                        gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGGGCGGCGGTCGATATTCCTGCCGGGATGCAGAGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGGTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCgggtttggggcgcgggggaaacggc                                                                                         >  1:287141/1‑226 (MQ=255)                                                                                                                                                                  gggacattggtggCGGGACTGCTGGATCTGCCGGTGGCAAAATTCATTTCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTGCTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATGTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGtttaag                                                                             <  2:686860/244‑4 (MQ=255)                                                                                                                                                                         tggtggCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgc                                                                                                         <  2:186724/216‑1 (MQ=255)                                                                                                                                                                           ggggggggAATGTTGGATCTGCCGGTGGCAAAATTCATTCCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATGCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTAttt                                                                    <  2:427604/246‑1 (MQ=255)                                                                                                                                                                           ggggCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGTTTGCGACGGCGGTATTTGTGTGGGTTGGTGGCGGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTAttt                                                                    <  1:307781/249‑1 (MQ=255)                                                                                                                                                                             ggCGGGAATGCTGGATCTGGCGGGGGCAAAATTCATTACGCCGAATATTAACGGCTGCCTGCTGTGGCCACCGTTTTGCTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGt                                                                  <  2:670979/251‑1 (MQ=255)                                                                                                                                                                                gggAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgc                                                               <  1:168389/251‑3 (MQ=255)                                                                                                                                                                                 ggAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcg                                                              <  1:621154/251‑4 (MQ=255)                                                                                                                                                                                      gCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCggcc                                                         <  1:120214/251‑3 (MQ=255)                                                                                                                                                                                         ggagctgccgttgtcaaaaTTCATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCACCGTTTTACTGCCTGCCTGGGATTCGGGCGGGAGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGtt                                                      <  2:104046/235‑1 (MQ=255)                                                                                                                                                                                                gcggttgcaatattcattaCGCCGAATATTATCGGTGGCCTGCTGTGGCCGCCTTTTTACTTCCTGCCGGGGATTCGGGGGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGGTTGTtgtg                                               <  2:361390/236‑1 (MQ=255)                                                                                                                                                                                                          aatccattaCTCCGAATATTATCGGCTGCCTGTGGGGGCTACTGTGTTACTTCCTGCCTGGGATTCTGGGGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCg                                     <  2:543892/246‑1 (MQ=255)                                                                                                                                                                                                                   cGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCt                                                                                                                     <  1:340579/162‑1 (MQ=255)                                                                                                                                                                                                                        aaTATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCAt                       <  1:614108/251‑1 (MQ=255)                                                                                                                                                                                                                                cggcggCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGTGGTGAGGTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTggctggct                                                                                                                          <  1:703100/144‑1 (MQ=255)                                                                                                                                                                                                                                  gCTGCCTGCTGGGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggtg             <  2:248788/251‑1 (MQ=255)                                                                                                                                                                                                                                     gcctgcTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgc           <  1:639679/250‑1 (MQ=255)                                                                                                                                                                                                                                             gggCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt  <  1:156471/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      GGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGT  >  NC_000913/71565‑72051

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