Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	944,547	G→A	S172N (AGT→AAT)	dmsC →	dimethyl sulfoxide reductase subunit C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	944,547	0	G	A	97.8%	146.1 / NA	45	S172N (AGT→AAT)	dmsC	dimethyl sulfoxide reductase subunit C
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/19);  minor base T (0/1);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAGTCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCT  >  NC_000913/944300‑944785                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               aGCGGTTCGATCTTGTTGGCCGTAGTCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGTCTGCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATct                                                                                                                                                                                                                                             <  1:231218/251‑1 (MQ=255)         gATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCACTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCg                                                                                                                                                                                                                                     >  2:69228/1‑251 (MQ=255)             ttttttGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCg                                                                                                                                                                                                                                  >  1:221868/1‑250 (MQ=255)             ttttttGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCg                                                                                                                                                                                                                                  >  2:355092/1‑250 (MQ=255)                 ttGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATgg                                                                                                                                                                                                                              <  1:261007/250‑1 (MQ=255)                   gCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCt                                                                                                                                                                                                                            <  2:93656/250‑1 (MQ=255)                           cggcATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACATCATCTGGACGCCGATGGGCGTCTTTCTGGCGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTGGTTGAATCTGGCTGGTGTCGAGGGCTGGGCGATGc                                                                                                                                                                                                                   <  2:182171/251‑1 (MQ=255)                                 cggcggGTTGCTGGCAATGCTGAAAAAGCTGTCGCCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCGTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCAGTGATACCGTTCCGACCTGGTACAGCATCTGGACTCCGATGGGCTTCTTCCTGACGATGGGTATGGGCGGCCCGTTACTGTGTTACCTGTTGTTGATTCTGGCGAGGGTCGATGGCGGGGCgggtcgtctg                                                                                                                                                                                                              <  2:228202/250‑10 (MQ=255)                                      ggTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGcc                                                                                                                                                                                                         >  1:229342/1‑250 (MQ=255)                                           cTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAAt                                                                                                                                                                                                   >  1:179201/1‑251 (MQ=255)                                              gCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTc                                                                                                                                                                                                <  1:80865/251‑1 (MQ=255)                                                      gAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTgg                                                                                                                                                                                        >  2:287066/1‑251 (MQ=255)                                                        aaaaGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCa                                                                                                                                                                                      <  2:82212/251‑1 (MQ=255)                                                                  aCCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTagtg                                                                                                                                                                             >  1:23780/1‑250 (MQ=255)                                                                                  cTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGt                                                                                                                                                            >  1:85192/1‑251 (MQ=255)                                                                                     tgGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGAt                                                                                                                                                         >  1:272596/1‑251 (MQ=255)                                                                                     tgGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGAt                                                                                                                                                         >  2:272702/1‑251 (MQ=255)                                                                                         tGATAGAGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGc                                                                                                                                                       >  2:37563/1‑249 (MQ=255)                                                                                                 aCGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGgcgc                                                                                                                                                >  2:127961/1‑248 (MQ=255)                                                                                                     tGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGAAACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGCTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGCTGCGTCTGCTGCCAGCAATTTCTGTACTGGCCCTGGTAGTGAGTGGCGTGGCGTCAGTGATGCAGGGCGCAGAGCTg                                                                                                                                         >  1:47981/1‑251 (MQ=255)                                                                                                                gTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTc                                                                                                                               >  2:52363/1‑250 (MQ=255)                                                                                                                  cATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCAt                                                                                                                             <  2:306700/250‑1 (MQ=255)                                                                                                                            tGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTgc                                                                                                                   >  1:377273/1‑250 (MQ=255)                                                                                                                                    ggTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGcc                                                                                                           <  1:373140/250‑1 (MQ=255)                                                                                                                                         gtgtgtATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAgcg                                                                                                     <  1:5010/251‑1 (MQ=255)                                                                                                                                            tgtATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGAAGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATATGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTAATAGTTCTGTGCAGCAGGCCGCAGCGCTg                                                                                                  >  1:222763/1‑251 (MQ=255)                                                                                                                                             gtATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTgg                                                                                                 <  2:7764/251‑1 (MQ=255)                                                                                                                                                                gTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCt                                                                              <  2:201777/251‑1 (MQ=255)                                                                                                                                                                       ccTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTa                                                                                                                                                                                            <  2:94568/134‑1 (MQ=255)                                                                                                                                                                       ccTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTa                                                                                                                                                                                            >  1:94568/1‑134 (MQ=255)                                                                                                                                                                                        gcGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTTGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCtttt                                                      <  2:74085/250‑1 (MQ=255)                                                                                                                                                                                        aCGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCtttt                                                      >  1:237829/1‑251 (MQ=255)                                                                                                                                                                                              aTGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGt                                                >  1:370306/1‑251 (MQ=255)                                                                                                                                                                                                    ttcttcCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCt                                          >  1:291748/1‑251 (MQ=255)                                                                                                                                                                                                         ccTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCt                                     <  1:29563/251‑1 (MQ=255)                                                                                                                                                                                                             aCGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTtgtg                                 <  1:268456/251‑1 (MQ=255)                                                                                                                                                                                                                      attggcggcCCGTTACTGGTGTAACTGTTGTTGAATATGGGTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAAGTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGAGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGGCGTTGCCCTGTGCTTGTGTAttcctct                        <  2:102356/248‑6 (MQ=255)                                                                                                                                                                                                                         ggcggcCCGTTACTGGGTTACCTTGTGTTGAATCTGGCTGGCGTCGATGGCGGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTTGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTAGTTTGCGCcaca                     <  1:303842/251‑1 (MQ=255)                                                                                                                                                                                                                              cccGTTACTGGGTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCCGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTa                 >  1:297975/1‑250 (MQ=255)                                                                                                                                                                                                                                   tACTGGGTTACCTGTTGGTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAggg           >  2:152811/1‑250 (MQ=255)                                                                                                                                                                                                                                       gggTTACCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTa         <  1:290237/249‑1 (MQ=255)                                                                                                                                                                                                                                       gggTTACCTGGTGCTGAATATGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGCGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCCTGGGGGTTGCACCACAGCTAAAAGGTTATc       >  2:54127/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCCTGTTGTTGAATCTGGCTGGCGTCGATGGGTGGGCGATGTGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGGGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGATGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAACAGGATATCAGCCt  >  1:254098/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCCTGTTGTTGAATCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCt  >  1:257251/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCCTGTTGTGGAATCGGGCGGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTGCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACGATGGGGCGCTGAGGTCCTGGCGGAGCGTGCTTTTGGCCGTTGCGCTGTGCGTGTGGATTGCACGACAGCTAAAAGGTTATCAGcc   <  2:175559/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               AGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAGTCTGGCTGGCGTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCT  >  NC_000913/944300‑944785

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