Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,775,361	Δ1 bp	coding (938/1149 nt)	mtlD →	mannitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,775,360	0	G	.	100.0%	156.3 / NA	37	coding (937/1149 nt)	mtlD	mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/17);  total (20/17)

AAAGATGATGTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTCGGTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCAATGCACTTCCGCAGTG  >  NC_000913/3775282‑3775432                                                                               |                                                                         aaaGATGATGTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTCg                                                                          <  1:673940/79‑1 (MQ=255)  aaGATGATGTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTCgtac                                                                       >  1:673995/1‑78 (MQ=255)    gatgatGTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGc                                                                    >  1:673998/1‑81 (MQ=255)      tgatgTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTg                                                                  >  2:674011/1‑81 (MQ=255)        atgTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTgg                                                                 <  1:673944/80‑1 (MQ=255)         tgTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGaa                                                               <  1:673945/81‑1 (MQ=255)           tAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAata                                                             >  1:673983/1‑81 (MQ=255)             gagCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAatat                                                            <  2:673945/80‑1 (MQ=255)             gagCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATg                                                           >  1:674004/1‑81 (MQ=255)               gcgcGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGt                                                         <  2:673946/81‑1 (MQ=255)                  cgTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTg                                                      >  2:674012/1‑81 (MQ=255)                    tAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGc                                                     >  1:673999/1‑80 (MQ=255)                     aGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGTCACTGCTC‑GTACGCTGGAATATGGTCTGCca                                                   >  2:674004/1‑81 (MQ=255)                      ggCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCcac                                                  <  1:673949/81‑1 (MQ=255)                      ggCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCcac                                                  <  1:673948/81‑1 (MQ=255)                        ccGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACAt                                                <  2:673949/81‑1 (MQ=255)                         cGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACAt                                                >  2:673993/1‑80 (MQ=255)                          gTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATa                                               >  2:674015/1‑80 (MQ=255)                             aGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATaaaaa                                           >  1:674011/1‑81 (MQ=255)                               ccACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAAcc                                         <  2:673950/81‑1 (MQ=255)                                     cGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGAt                                     <  2:673951/79‑1 (MQ=255)                                      gTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTg                                   >  2:674013/1‑80 (MQ=255)                                       tAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGaa                                 <  1:673953/81‑1 (MQ=255)                                         aaCTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAgg                               <  2:673954/81‑1 (MQ=255)                                         aaCTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAg                                <  1:673954/80‑1 (MQ=255)                                          aCTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAgg                               >  1:674018/1‑80 (MQ=255)                                           cTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTa                             >  1:674024/1‑81 (MQ=255)                                              aGTGCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTAtt                           <  1:673956/80‑1 (MQ=255)                                                 gCTGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCg                       >  1:674008/1‑81 (MQ=255)                                                   tGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCt                     >  2:673956/1‑81 (MQ=255)                                                   tGGCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCt                     >  1:674009/1‑81 (MQ=255)                                                     gCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTg                    <  2:673957/80‑1 (MQ=255)                                                     gCGACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGc                   <  2:673958/81‑1 (MQ=255)                                                       gACCGTCTGATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCaa                 >  1:674032/1‑81 (MQ=255)                                                               gATCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCAATGCACTTc         >  2:674024/1‑81 (MQ=255)                                                                aTCAAGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCAATGCACtt          <  1:673960/79‑1 (MQ=255)                                                                   aaGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCAATGCACTTCCg       >  1:674017/1‑79 (MQ=255)                                                                   aaGCCACTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCAATGCACTTCCGCa     <  1:673961/81‑1 (MQ=255)                                                                        aCTGCTC‑GTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCAATGCACTTCCGCAGTg  >  1:674035/1‑79 (MQ=255)                                                                               |                                                                         AAAGATGATGTAGAGCGCGTAGGCCGTCAGCCACTGCGTAAACTGAGTGCTGGCGACCGTCTGATCAAGCCACTGCTCGGTACGCTGGAATATGGTCTGCCACATAAAAACCTGATTGAAGGTATTGCCGCTGCAATGCACTTCCGCAGTG  >  NC_000913/3775282‑3775432

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