Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	104,426	C→T	L149F (CTT→TTT)	ftsA →	cell division protein FtsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	104,426	0	C	T	100.0%	138.5 / NA	39	L149F (CTT→TTT)	ftsA	cell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/25);  total (14/25)

GATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCAAAGC  >  NC_000913/104348‑104502                                                                               |                                                                             gATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGAttt                                                                            >  2:25335/1‑81 (MQ=255) gATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGAttt                                                                            >  1:523528/1‑81 (MQ=255) gATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTACGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGAttt                                                                            >  1:25352/1‑81 (MQ=255)  aTGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGAtttt                                                                           >  1:25365/1‑81 (MQ=255)   tGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTc                                                                          <  2:25287/81‑1 (MQ=255)   tGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTc                                                                          <  2:25286/81‑1 (MQ=255)         tCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGt                                                                    <  1:25288/81‑1 (MQ=255)              tgctgcATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCgg                                                                <  2:25289/80‑1 (MQ=255)              tgctgcATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGa                                                               <  2:25290/81‑1 (MQ=255)                ctgcATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATg                                                             >  2:108705/1‑81 (MQ=255)                    aTGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAgg                                                         <  2:25291/81‑1 (MQ=255)                     tgtgATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGc                                                        <  2:25292/81‑1 (MQ=255)                       tgATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCaa                                                      <  2:25293/81‑1 (MQ=255)                         aTCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCaaaa                                                    <  2:25294/81‑1 (MQ=255)                           cccGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGt                                                  >  1:25359/1‑81 (MQ=255)                            ccGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGt                                                  <  1:25295/80‑1 (MQ=255)                            ccGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGt                                                  <  1:25296/80‑1 (MQ=255)                                aaGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAcc                                             <  2:25297/81‑1 (MQ=255)                                  gagTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGTGCGTGCGGATGCAGGCAAAAGTGGACCTg                                           <  2:25298/81‑1 (MQ=255)                                   agTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                                          >  1:25364/1‑81 (MQ=255)                                   agTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                                          <  1:25299/81‑1 (MQ=255)                                    gTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGAt                                         <  2:25301/81‑1 (MQ=255)                                    gTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGAt                                         >  2:25368/1‑81 (MQ=255)                                    gTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGAt                                         <  2:25300/81‑1 (MQ=255)                                    gTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGAt                                         <  1:25302/81‑1 (MQ=255)                                              gACTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTca                                >  1:25357/1‑80 (MQ=255)                                               aCTATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTcaca                              <  1:25303/81‑1 (MQ=255)                                                 tATCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAAc                            >  2:25363/1‑81 (MQ=255)                                                   tCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAAc                            <  1:25304/79‑1 (MQ=255)                                                   tCAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGa                          >  2:25361/1‑81 (MQ=255)                                                    cAGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGTAAAAGTGCACCTGATCACATGTCACAACGa                          >  2:25358/1‑80 (MQ=255)                                                     aGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGata                        >  1:25370/1‑81 (MQ=255)                                                     aGGAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGat                         <  2:25305/80‑1 (MQ=255)                                                      ggAAGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGatat                       >  2:25343/1‑81 (MQ=255)                                                        aaGGGATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATgg                     >  2:25356/1‑81 (MQ=255)                                                           ggATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCg                   >  1:25375/1‑80 (MQ=255)                                                           ggATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGa                  <  2:25306/81‑1 (MQ=255)                                                           ggATCAAGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGa                  <  2:25307/81‑1 (MQ=255)                                                           ggATCAAGAATCCGGTAGGATCTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCg                   >  2:25362/1‑80 (MQ=255)                                                                aaGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa              <  2:25308/80‑1 (MQ=255)                                                                aaGAATCCGGTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGaaaaa              <  2:25309/80‑1 (MQ=255)                                                                        ggTAGGATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCa      >  1:25372/1‑80 (MQ=255)                                                                            ggATTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCAAAGc  <  1:25310/80‑1 (MQ=255)                                                                               |                                                                             GATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCGAAAAACATCGTCAAAGC  >  NC_000913/104348‑104502

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