Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,981,145	A→C	I156M (ATT→ATG)	otsA ←	trehalose‑6‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,981,145	0	A	C	100.0%	98.2 / NA	28	I156M (ATT→ATG)	otsA	trehalose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/8);  total (20/8)

ATCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATAA  >  NC_000913/1981070‑1981220                                                                            |                                                                            aTCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATg                                                                        <  2:348149/81‑1 (MQ=255)    aTAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATg                                                                        >  2:348189/1‑78 (MQ=255)         cacaAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCagaaag                                                                 >  2:348175/1‑80 (MQ=255)          acaAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCagaaaga                                                                <  2:348150/80‑1 (MQ=255)            aaaGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAAc                                                             >  1:348179/1‑81 (MQ=255)             aaGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAAcc                                                            <  2:348151/81‑1 (MQ=255)                      aagcaagGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGat                                                    >  1:348196/1‑80 (MQ=255)                        gcaagGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGattat                                                 >  1:348176/1‑81 (MQ=255)                           agGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCa                                              >  1:348180/1‑81 (MQ=255)                                   tatGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGt                                      >  2:348177/1‑81 (MQ=255)                                    atGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGtt                                     >  2:348172/1‑81 (MQ=255)                                      gTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTg                                   >  1:348178/1‑81 (MQ=255)                                         ggCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCg                                 >  1:348184/1‑80 (MQ=255)                                         ggCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCg                                 >  1:348188/1‑80 (MQ=255)                                              cgcgTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTc                           <  1:348152/81‑1 (MQ=255)                                                cgTTGAAGATTTCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCAt                         >  2:348182/1‑81 (MQ=255)                                                         tttCCGGTGTCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATg                <  2:348153/81‑1 (MQ=255)                                                                tgtCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAg         <  1:348154/81‑1 (MQ=255)                                                                 gtCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAgg        >  1:348195/1‑81 (MQ=255)                                                                  tCGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGt       >  1:348187/1‑81 (MQ=255)                                                                   cGGGAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGt       >  2:348191/1‑80 (MQ=255)                                                                     ggAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGa     >  1:348181/1‑80 (MQ=255)                                                                     ggAAAGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGAt    <  1:348155/81‑1 (MQ=255)                                                                      gtaaGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGAt    >  2:348183/3‑80 (MQ=255)                                                                       aaaGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATaa  >  1:348202/1‑81 (MQ=255)                                                                       aaaGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATaa  >  1:348185/1‑81 (MQ=255)                                                                       aaaGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATaa  <  1:348156/81‑1 (MQ=255)                                                                       aaaGGCATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATaa  >  2:348203/1‑81 (MQ=255)                                                                            |                                                                            ATCATAATCACAAAGCTGTTCAAGCAAGGTGTCATATGTCGGCAGCGCGTTGAAGATTTCCGGTGTCGGGAAAGGAATATGCAGAAAGAAACCAATGCGATTATTCACTCCCCGTTTGCGTAATTCATGCGCAAATGGCAACAGGTGATAA  >  NC_000913/1981070‑1981220

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