Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	941,941	C→T	P328L (CCA→CTA)	dmsA →	dimethyl sulfoxide reductase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,941	0	C	T	100.0%	191.2 / NA	54	P328L (CCA→CTA)	dmsA	dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/28);  total (26/28)

AGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTC  >  NC_000913/941809‑942072                                                                                                                                     |                                                                                                                                    aGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGt                                                                                                                              <  2:1080388/140‑1 (MQ=255)        cccTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGa                                                                                                                       >  1:466406/1‑140 (MQ=255)             ttCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAAt                                                                                                                  <  2:833704/140‑1 (MQ=255)             ttCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAAt                                                                                                                  <  1:577361/140‑1 (MQ=255)                gTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGc                                                                                                               >  2:949753/1‑140 (MQ=255)                gTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGc                                                                                                               <  1:579172/140‑1 (MQ=255)                gTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGc                                                                                                               <  1:1148379/140‑1 (MQ=255)                gTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGc                                                                                                               <  1:1137822/140‑1 (MQ=255)                  ccGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCa                                                                                                             <  2:445540/140‑1 (MQ=255)                  ccGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCa                                                                                                             >  2:456156/1‑140 (MQ=255)                     ggAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACta                                                                                                          >  1:305048/1‑140 (MQ=255)                      gAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACtat                                                                                                         <  1:1129726/140‑1 (MQ=255)                        aCAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATaa                                                                                                       <  2:305048/140‑1 (MQ=255)                                  cACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCt                                                                                                   >  2:2295/1‑134 (MQ=255)                                  cACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCt                                                                                                   <  1:2295/134‑1 (MQ=255)                                         tAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTg                                                                                      <  2:154471/140‑1 (MQ=255)                                          aaCGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGa                                                                                     >  2:300096/1‑140 (MQ=255)                                            cGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAg                                                                                   <  1:1200579/140‑1 (MQ=255)                                             ggTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAgg                                                                                  <  2:1204071/140‑1 (MQ=255)                                               tCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGc                                                                                >  2:123228/1‑140 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTg                                                                                      >  1:462429/1‑132 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTg                                                                                      <  2:462429/132‑1 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAgg                                                                                  <  1:1189293/136‑1 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAgg                                                                                  >  2:1189293/1‑136 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCa                                                                              >  2:1206709/1‑140 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCa                                                                              >  1:1083911/1‑140 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCa                                                                              >  1:821878/1‑140 (MQ=255)                                                 tGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCa                                                                              >  1:177014/1‑140 (MQ=255)                                                       aCGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGc                                                                        <  1:1206709/140‑1 (MQ=255)                                                        cGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCg                                                                       >  2:262359/1‑140 (MQ=255)                                                                   cTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACg                                                            <  1:955104/140‑1 (MQ=255)                                                                   cTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACg                                                            >  2:791682/1‑140 (MQ=255)                                                                    tGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAg                                                                             <  2:1096835/122‑1 (MQ=255)                                                                    tGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAg                                                                             >  1:1096835/1‑122 (MQ=255)                                                                      aaaaCCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCg                                                         >  1:1168397/1‑140 (MQ=255)                                                                           cTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATg                                                    >  2:724997/1‑140 (MQ=255)                                                                             ggtggATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATggg                                                  <  1:724997/140‑1 (MQ=255)                                                                               tggATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGcc                                                <  2:865000/140‑1 (MQ=255)                                                                               tggATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGcc                                                <  2:466406/140‑1 (MQ=255)                                                                                    cAGGCCTTCCTCGATAAATATTGCGTAGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGTCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCa                                           <  2:1083911/140‑1 (MQ=255)                                                                                    cAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCa                                           <  1:300096/140‑1 (MQ=255)                                                                                       gCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAAt                                        <  1:9664/140‑1 (MQ=255)                                                                                       gCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAAt                                        <  2:1168397/140‑1 (MQ=255)                                                                                       gCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCAGCTAGCAAGTGCGACGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAAt                                        <  2:6486/140‑1 (MQ=255)                                                                                         aTTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCa                                      >  2:139449/1‑140 (MQ=255)                                                                                           tCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACt                                    >  1:542710/1‑140 (MQ=255)                                                                                            ccTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTg                                   >  2:353897/1‑140 (MQ=255)                                                                                            ccTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTg                                   >  1:661767/1‑140 (MQ=255)                                                                                            ccTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTg                                   >  1:1083555/1‑140 (MQ=255)                                                                                               cGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGgtg                                <  2:542710/140‑1 (MQ=255)                                                                                                   aaaTACTGCGTTGGCTACGATGAGAACACCCTGCTAGCCATTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTTAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTcc                            <  2:1083555/140‑1 (MQ=255)                                                                                                         tGCGTTGGCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGa                      >  1:1126926/1‑140 (MQ=255)                                                                                                                gCTACGATGAGAAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAatc               >  1:1136488/1‑140 (MQ=255)                                                                                                                          gaAAACCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTgg     >  1:1214137/1‑140 (MQ=255)                                                                                                                             aaCCCTGCTAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGGAGACAAAATCATCAAATTGGCTc  <  1:139449/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    AGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTTATATTCTGGGTGAAGGGCCAGATGGCGTGGCTAAAACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTC  >  NC_000913/941809‑942072

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