Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	943,967	C→T	P185L (CCC→CTC)	dmsB →	dimethyl sulfoxide reductase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	943,967	0	C	T	100.0%	222.1 / NA	69	P185L (CCC→CTC)	dmsB	dimethyl sulfoxide reductase subunit B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/36);  total (33/36)

ACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATC  >  NC_000913/943724‑944212                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      aCTGCCATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATg                                                                                                                                                                                                                                                  <  2:526618/249‑1 (MQ=255) aCTGCCATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATg                                                                                                                                                                                                                                                  <  2:169624/249‑1 (MQ=255)      cATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAg                                                                                                                                                                                                                                           <  1:448322/251‑1 (MQ=255)      cATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGAGCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACATTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAg                                                                                                                                                                                                                                           <  2:324956/251‑1 (MQ=255)        tatGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATAATACAAACTGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGCTGAATCATGTCCGCTGCGGGCGCTGGATTTCCGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAgcc                                                                                                                                                                                                                                         >  1:77659/1‑251 (MQ=255)            gCTTGCCCGTATAGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGccc                                                                                                                                                                                                                                     <  2:158597/251‑1 (MQ=255)                      aTGGCGCACCTCACTATAACTAAACGAAAGGCCATATGGCCAAATGCGATGGTTGTTATGGCCGTGTTGCGGAGGGTAAAAAGTCGTTCTGTCTTGAATCTTGTCTGCTGCGGGTGCTTGATTTCTGCCCTATCGACGAGCTTCTTAAAAATCATGGCGATCTTGTTGCCGTGTCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGTCCGGGGAt                                                                                                                                                                                                                           <  2:114702/251‑1 (MQ=255)                        ggCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATAc                                                                                                                                                                                                                         <  1:140619/251‑1 (MQ=255)                            cACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTTATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACt                                                                                                                                                                                                                     >  1:427879/1‑251 (MQ=255)                            cACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACt                                                                                                                                                                                                                     <  2:29400/251‑1 (MQ=255)                               cggaaTATAACGAAGCGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCAGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTGACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCACGACCGGGGATACCGCTGGc                                                                                                                                                                                                                  <  2:466418/248‑1 (MQ=255)                                gCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCt                                                                                                                                                                                                                 >  1:84120/1‑251 (MQ=255)                                 cAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCt                                                                                                                                                                                                                 >  1:133322/1‑250 (MQ=255)                                 cAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTa                                                                                                                                                                                                                <  1:331530/251‑1 (MQ=255)                                     ataACGAAACGCAAGGCCATTTGACCAAATGCGATGGTTGTTATGACCGTGTTCCGGAGGGAAAAAAGACGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGCTATCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTATGATCAAACTCCATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTAtatg                                                                                                                                                                                                            >  2:364345/1‑248 (MQ=255)                                     ataACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCTGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGAGCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACGACTTTCTATCTg                                                                                                                                                                                                            <  2:657063/251‑1 (MQ=255)                                       aacgaaacgaAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTgg                                                                                                                                                                                                           <  1:254345/250‑1 (MQ=255)                                             acgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAAccc                                                                                                                                                                                                    <  1:537944/251‑1 (MQ=255)                                               gaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGa                                                                                                                                                                                                  <  1:133245/251‑1 (MQ=255)                                                aaaGGACATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCGTGTCCGCTGCTGGCTCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAACCCCGaa                                                                                                                                                                                                 <  1:108215/251‑1 (MQ=255)                                                      cATATGCCCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGCTCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCACTAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGgt                                                                                                                                                                                           >  2:265164/1‑251 (MQ=255)                                                        tatGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGgtgt                                                                                                                                                                                         >  1:189814/1‑251 (MQ=255)                                                             ccAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGACGGAGGtgtg                                                                                                                                                                                        >  2:314873/1‑247 (MQ=255)                                                                  tGCGAGGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCTGATCTGTGTTGAATCTTGTCCGCTGCGGGCGATGGATATCGGACTTATCGAAGAGCTGCGTAAAAAACGTGGCGATATGGCGGCCGTTGCGCCGTTGCAGCCAGCTCACTTAACAACACCGAACATTGTGAAGAAACTCAATGCCAATAGCCGCCCGACCGGGGCGATCACTGGCTATCTGGCAAACCCGCAGGAGGTGAGCGATgggt                                                                                                                                                                                >  2:256220/1‑249 (MQ=255)                                                                        ggTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGAt                                                                                                                                                                         >  1:52760/1‑251 (MQ=255)                                                                         gttgttATGACCGTGTTGCGGAGGGCAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGAGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGCTGGGAAGTGGATg                                                                                                                                                                        >  1:420550/1‑251 (MQ=255)                                                                           tgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATgg                                                                                                                                                                       <  1:240572/250‑1 (MQ=255)                                                                                tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGCATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGCGTGAGATGGGAAGTGGATGGCATGa                                                                                                                                                                  >  2:495215/1‑250 (MQ=255)                                                                                 gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAAt                                                                                                                                                                >  2:582410/1‑251 (MQ=255)                                                                                  aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATg                                                                                                                                                               >  2:119849/1‑251 (MQ=255)                                                                                    cGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGc                                                                                                                                                             >  1:447541/1‑251 (MQ=255)                                                                                     gtgtTGCGGAGGGTAAAAAGCCGATCTGTGTAGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGcc                                                                                                                                                            <  1:239117/251‑1 (MQ=255)                                                                                           cGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGCGATTGGAAGTGGATGGCATGCATGGCCGCtgat                                                                                                                                                      <  2:227876/251‑1 (MQ=255)                                                                                            ggAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCtgatg                                                                                                                                                     <  2:399766/251‑1 (MQ=255)                                                                                            ggAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCtgatg                                                                                                                                                     >  1:41855/1‑251 (MQ=255)                                                                                               gggTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATc                                                                                                                                                  >  1:270054/1‑251 (MQ=255)                                                                                                                    tGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTCAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATTTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGGGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAAtgt                                                                                                                             >  1:630423/1‑251 (MQ=255)                                                                                                                         cgtGTCCGCTGCGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTTCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATCCCACTGGCTATCTGGCAACCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTTGGTCAATGTGTAg                                                                                                                         <  2:48741/248‑1 (MQ=255)                                                                                                                          ttGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCAGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGc                                                                                                                        >  1:585059/1‑250 (MQ=255)                                                                                                                           tGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGc                                                                                                                        >  2:358990/1‑249 (MQ=255)                                                                                                                               cGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAggtgg                                                                                                                  <  1:204776/251‑1 (MQ=255)                                                                                                                                   gCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATTGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGtttt                                                                                                              <  2:189660/251‑1 (MQ=255)                                                                                                                                    ggggCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATGCCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGGGGTTTTa                                                                                                             <  2:314630/250‑1 (MQ=255)                                                                                                                                                        tATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCTAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCgctg                                                                                         >  1:435776/1‑251 (MQ=255)                                                                                                                                                         aTCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCgctgc                                                                                        <  1:528404/251‑1 (MQ=255)                                                                                                                                                         aTCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCgctg                                                                                         <  1:234286/250‑1 (MQ=255)                                                                                                                                                            gacgaGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCa                                                                                     <  1:548203/251‑1 (MQ=255)                                                                                                                                                               gaGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCACTGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGCAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTTTGGCTTTGGCGCTGCTCaa                                                                                    >  2:626029/1‑249 (MQ=255)                                                                                                                                                                    gCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGAc                                                                             >  1:413392/1‑251 (MQ=255)                                                                                                                                                                         aaaaaCATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCg                                                                        >  1:546589/1‑251 (MQ=255)                                                                                                                                                                         aaaaaCATGGCGATCTGGCGGCAGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAATCAATGCCAATAGCCGCCCGACCGGGGATACCACAGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGGGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGGGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCg                                                                        >  2:419176/1‑251 (MQ=255)                                                                                                                                                                           aaaCATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCg                                                                        >  1:122629/1‑249 (MQ=255)                                                                                                                                                                           aaaCATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCg                                                                        >  1:551491/1‑249 (MQ=255)                                                                                                                                                                                        tGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCcagc                                                          >  1:610733/1‑250 (MQ=255)                                                                                                                                                                                               cgttgcgccgttgcCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGtt                                                  >  2:170994/1‑251 (MQ=255)                                                                                                                                                                                                gttgcgccgttgcCGCGAGCTCACTTTACCAAACCGAATATTGTGATAACACTCACTGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATATGGCAAACACGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGACTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGtt                                                  >  1:84181/1‑250 (MQ=255)                                                                                                                                                                                                   gcgccgttgcCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCg                                              <  2:172140/251‑1 (MQ=255)                                                                                                                                                                                                     gccgttgcCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGGGGATGGCATGAATGGCCGCTGATGATCGTCACGGTGTTCGGGCAATGTGTAGCAGGTGGTTGTATCGTTCTGGCTGGGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCAAGCGGCGTGTTATCGcc                                            <  2:630423/251‑1 (MQ=255)                                                                                                                                                                                                           tccGCGAGCTCACTTTACCAAACCGAATATTGTGATCACACTCAATGCCAATAGCCGCCCGACCGGGGATACCAATGGCGATATGGCGAACCCGAAGGAGGTGTGAGATGGGAAGTGGCTGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGAGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATg                                      >  2:523041/2‑251 (MQ=255)                                                                                                                                                                                                                 aGCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGt                                <  1:582410/251‑1 (MQ=255)                                                                                                                                                                                                                  gCTCACTTTACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGtt                               <  2:141359/251‑1 (MQ=255)                                                                                                                                                                                                                          tACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCg                                                                               <  2:313773/195‑1 (MQ=255)                                                                                                                                                                                                                          tACCAAACCGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCg                                                                               >  1:313773/1‑195 (MQ=255)                                                                                                                                                                                                                             cAAAACGAATATTGTGATCAAACTCAATGCCAATAGCCGCCAGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGAGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCAAGGTCTTCGGGCAATGGGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGAGTTATCGCCTGCATGTTTGGTTTATGGGTGCTg                    >  2:618878/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGAATATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATggg               <  2:270054/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaTATTGTGATCAAACTCAATGCCAATAGCCGCTCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAATTGGATGGCATGAATGTCCGCTGGTGAGCTTCACGGTCGTCGTGCAATGTGTAGCAGGGGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCATCAGCGTGTGATCGCCTGCAGGTGTGGTTTATGGGTGCTGATGGGCa             <  2:268978/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaTATTGTGATCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCa             <  1:215079/251‑1 (MQ=255)                                                                                                                                                                                                                                            taTCAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCAATGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACTGTCTTCGGGCAATGTGTAGCAGATGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGTATGTTTGGTTTATGGGTGCTGATGGGCACTGGCtt      <  2:427879/249‑1 (MQ=255)                                                                                                                                                                                                                                               cAAACTCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATc  <  2:610733/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      ACTGCCATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATC  >  NC_000913/943724‑944212

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