Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,618,445	A→G	V155A (GTA→GCA)	marC ←	UPF0056 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,618,445	0	A	G	100.0%	159.1 / NA	43	V155A (GTA→GCA)	marC	UPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/22);  total (21/22)

ATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTG  >  NC_000913/1618298‑1618591                                                                                                                                                    |                                                                                                                                                   atacatacCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGccc                                                                                                                                                 >  2:159890/1‑151 (MQ=255)        ccAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCCCCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCa                                                                                                                                            <  1:218530/149‑1 (MQ=255)                aaTCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTg                                                                                                                                  >  2:380213/1‑151 (MQ=255)                          gACGGGAGATCGCTTCAATGCCCCCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACg                                                                                                                        <  2:318344/151‑1 (MQ=255)                          gACGGGAGATCGCTTCAATGCCACCCTTGCCCCCCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACg                                                                                                                        <  2:467906/151‑1 (MQ=255)                              ggAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACt                                                                                                                    <  2:699647/151‑1 (MQ=255)                              ggAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACt                                                                                                                    <  1:437523/151‑1 (MQ=255)                                  aTCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTgg                                                                                                                >  1:302218/1‑151 (MQ=255)                                     gCTTCAATGCCCCCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACg                                                                                                             <  2:308752/151‑1 (MQ=255)                                     gCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACg                                                                                                             >  1:145497/1‑151 (MQ=255)                                              ccACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCta                                                                                                     >  2:525601/1‑150 (MQ=255)                                              ccACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGAGGAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCtaa                                                                                                    <  1:429838/151‑1 (MQ=255)                                                 cccTTGCCCACCAGCCGCATAATTGCGCGGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCtaataa                                                                                                 <  1:217288/151‑1 (MQ=255)                                                 cccTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCtaata                                                                                                  <  2:50492/150‑1 (MQ=255)                                                 cccTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCtaata                                                                                                  >  2:258113/1‑150 (MQ=255)                                                    ttGCCCACCAGCCGCATAATTCCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGGTGAACTCTGACGCACTGTGGACGCCGAGCTaataat                                                                                                >  1:670266/1‑149 (MQ=255)                                                       cccACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCg                                                                                           <  1:145793/151‑1 (MQ=255)                                                           ccaGCCGCATAATTGCGCCGGAGCTGCTTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTa                                                                                        <  1:170901/150‑1 (MQ=255)                                                           ccaGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTa                                                                                        <  1:752287/150‑1 (MQ=255)                                                           ccaGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTAt                                                                                       <  2:413543/151‑1 (MQ=255)                                                                     aaTTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATACCCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGa                                                                             <  1:481854/151‑1 (MQ=255)                                                                             cGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAGGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGt                                                                     >  1:256176/1‑151 (MQ=255)                                                                                 gCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACtt                                                                 >  2:114060/1‑151 (MQ=255)                                                                                      gTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCAt                                                            >  2:126345/1‑151 (MQ=255)                                                                                             tCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAg                                                      >  2:158860/1‑150 (MQ=255)                                                                                                 cacaAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGc                                                  >  2:587671/1‑150 (MQ=255)                                                                                                    aaaaTAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGcaca                                               >  1:669655/1‑150 (MQ=255)                                                                                                       aaaaCCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAgg                                           <  1:326226/149‑1 (MQ=255)                                                                                                       aTAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAgg                                           >  2:133124/1‑151 (MQ=255)                                                                                                           ccGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAat                                       <  2:475693/151‑1 (MQ=255)                                                                                                             gCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAatat                                     <  1:380578/151‑1 (MQ=255)                                                                                                                 ccaGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTgg                                  <  2:391605/150‑1 (MQ=255)                                                                                                                  caGGAAAAAGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCa                                <  1:469220/151‑1 (MQ=255)                                                                                                                      aaaaaGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTgg                            >  2:388874/1‑151 (MQ=255)                                                                                                                      aaaaaGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTgg                            >  1:32098/1‑151 (MQ=255)                                                                                                                       aaaaGATCAGCGGCGGCGCAACCATAAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGAGCGGCAGTAATTGGCATTGCCAGCGGCACAAAGGAAATCTTGGCACAggg                           >  1:414000/1‑151 (MQ=255)                                                                                                                        aaaGATCAGCGGCGGCGCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGt                          <  2:295156/151‑1 (MQ=255)                                                                                                                                  ggcggcgCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTcc                 >  1:184448/1‑150 (MQ=255)                                                                                                                                     ggcgCAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGct             >  2:488299/1‑151 (MQ=255)                                                                                                                                       cgcAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGctct           <  2:277682/151‑1 (MQ=255)                                                                                                                                        gcAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCtt          <  1:629397/151‑1 (MQ=255)                                                                                                                                         cAACCATCAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTc         >  1:353146/1‑151 (MQ=255)                                                                                                                                                cAGTGCCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTg  >  2:455050/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGGCGCAACCATCAGTACCCAATCAGCGAAGGTTGAACTCTGACGCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTG  >  NC_000913/1618298‑1618591

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