Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,152,125	T→C	G116G (GGA→GGG)	pphC ←	protein‑serine/threonine phosphatase PphC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,152,125	0	T	C	100.0%	133.5 / NA	38	G116G (GGA→GGG)	pphC	protein‑serine/threonine phosphatase PphC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/20);  total (18/20)

AAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGTCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTAAGCCCCAATTCCCCGCCCTGC  >  NC_000913/2151997‑2152254                                                                                                                                 |                                                                                                                                  aaGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCAGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAAtgt                                                                                                                        <  1:145938/140‑1 (MQ=255)      cGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAgg                                                                                                                   >  2:678788/1‑140 (MQ=255)          aTGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCaaa                                                                                                               >  2:13476/1‑140 (MQ=255)                 cATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGa                                                                                                        <  2:634915/140‑1 (MQ=255)                          ccGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAAt                                                                                               >  1:564178/1‑140 (MQ=255)                           cGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAAtt                                                                                              >  2:369292/1‑140 (MQ=255)                                        agagGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGc                                                                                 >  1:550120/1‑140 (MQ=255)                                            gCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCg                                                                             <  1:334346/140‑1 (MQ=255)                                            gCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCg                                                                             >  2:11135/1‑140 (MQ=255)                                                       ccGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCg                                                                  <  1:631665/140‑1 (MQ=255)                                                           gACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGtt                                                                >  2:860749/1‑138 (MQ=255)                                                           gACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGtt                                                                <  1:860749/138‑1 (MQ=255)                                                             cccAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACg                                                            >  2:213410/1‑140 (MQ=255)                                                              ccAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGa                                                           <  2:550120/140‑1 (MQ=255)                                                                   tCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGt                                                      <  1:688734/140‑1 (MQ=255)                                                                   tCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGt                                                      <  2:353281/140‑1 (MQ=255)                                                                     tACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCa                                                    <  1:678788/140‑1 (MQ=255)                                                                      accaccACGCCGCCATCGCCATTCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAg                                                   >  2:613052/1‑140 (MQ=255)                                                                       ccaccaCGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCa                                                    <  2:614583/138‑1 (MQ=255)                                                                       ccaccaCGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCa                                                    >  1:614583/1‑138 (MQ=255)                                                                         accacGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTAc                                                <  2:827836/140‑1 (MQ=255)                                                                           cacGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCa                                              <  1:369292/140‑1 (MQ=255)                                                                             cgccgcCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGt                                                      >  1:5020/1‑130 (MQ=255)                                                                             cgccgcCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGt                                                      <  2:5020/130‑1 (MQ=255)                                                                               ccgccATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATAtt                                          >  2:695586/1‑140 (MQ=255)                                                                                 gccATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTtg                                        <  1:790021/140‑1 (MQ=255)                                                                                     tCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGc                                    <  2:200870/140‑1 (MQ=255)                                                                                      cGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCg                                   <  1:213410/140‑1 (MQ=255)                                                                                       gCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGa                                  <  1:68383/140‑1 (MQ=255)                                                                                          aaTCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGa                               <  1:695586/140‑1 (MQ=255)                                                                                            tCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGAca                             >  2:682659/1‑140 (MQ=255)                                                                                                cATAATCAAGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGAcatcat                         >  2:760729/1‑140 (MQ=255)                                                                                                        aGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCTGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTa                       >  1:180347/1‑134 (MQ=255)                                                                                                        aGGTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCTGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTa                       <  2:180347/134‑1 (MQ=255)                                                                                                         ggTGCCGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTAAGCCCCa                >  1:345208/1‑140 (MQ=255)                                                                                                             ccGTCAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTAAGCCCCAATTc            >  1:516099/1‑140 (MQ=255)                                                                                                                 cAGGCGATGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTAAGCCCCAATTCCCCg        <  2:564178/140‑1 (MQ=255)                                                                                                                       aTGAAATCAGCCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTAAGCCCCAATTCCCCGCCCTGc  >  2:406829/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                  AAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATGAAATCAGTCCCAGAAATGTGCAGGCAAAATCGCGAACCGCTAATTCTTTGGCTTCCGCCTCGGCAAAAAGTCGTTGACGAATGGTCAGTACCATATTTGTGGCGAGGACATCATTAAGCCCCAATTCCCCGCCCTGC  >  NC_000913/2151997‑2152254

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