Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	405,224	G→A	P141S (CCA→TCA)	proC ←	pyrroline‑5‑carboxylate reductase, NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	405,224	0	G	A	100.0%	132.4 / NA	42	P141S (CCA→TCA)	proC	pyrroline‑5‑carboxylate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/18);  total (24/18)

GCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGGGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGC  >  NC_000913/405115‑405353                                                                                                              |                                                                                                                                  gCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAgg                                                                                                      >  2:267704/1‑139 (MQ=255)       gaagaaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATcc                                                                                                <  2:1733836/139‑1 (MQ=255)           aaCCGCTCACACCGACCACTGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTgg                                                                                                                 >  2:1635774/1‑118 (MQ=255)           aaCCGCTCACACCGACCACTGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTgg                                                                                                                 <  1:1635774/118‑1 (MQ=255)           aaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGc                                                                                            >  1:944108/1‑139 (MQ=255)           aaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGc                                                                                            >  1:339196/1‑139 (MQ=255)           aaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTAAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGc                                                                                            >  2:1241369/1‑139 (MQ=255)                cTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTaa                                                                                       >  1:1553993/1‑139 (MQ=255)                 tCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAAc                                                                                      <  1:1241369/139‑1 (MQ=255)                       cGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGt                                                                                                     <  2:1472819/118‑1 (MQ=255)                       cGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGt                                                                                                     >  1:1472819/1‑118 (MQ=255)                          ccaccGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACg                                                                                                          <  1:657120/110‑1 (MQ=255)                          ccaccGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACg                                                                                                          >  2:657120/1‑110 (MQ=255)                          ccaccGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAAc                                                                                      >  1:1705262/1‑130 (MQ=255)                          ccaccGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAAc                                                                                      <  2:1705262/130‑1 (MQ=255)                                        atcGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATgg                                                               >  2:848474/1‑139 (MQ=255)                                           ggCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGgcgc                                                            >  2:435915/1‑139 (MQ=255)                                              tCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGACATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGa                                                         <  2:508509/139‑1 (MQ=255)                                                aGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATa                                                       >  1:1423605/1‑139 (MQ=255)                                                aGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATa                                                       >  1:608131/1‑139 (MQ=255)                                                  cAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAAt                                                     >  2:650560/1‑139 (MQ=255)                                                            cGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCa                                           >  2:712055/1‑139 (MQ=255)                                                               ttCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATgg                                        >  1:962777/1‑139 (MQ=255)                                                                       agcagcGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAgcgc                                <  2:1423605/139‑1 (MQ=255)                                                                        gcagcGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAgcgcg                               >  1:653824/1‑139 (MQ=255)                                                                           gcGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCgg                            >  1:1312100/1‑139 (MQ=255)                                                                            cGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAg                                                                                                       <  2:670659/63‑1 (MQ=255)                                                                            cGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAg                                                                                                       >  1:670659/1‑63 (MQ=255)                                                                             ggAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCAt                                                                                          >  2:1070407/1‑75 (MQ=255)                                                                             ggAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCAt                                                                                          <  1:1070407/75‑1 (MQ=255)                                                                                 aaTATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGc                      >  2:73026/1‑139 (MQ=255)                                                                                     ttCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGt                  <  2:944108/139‑1 (MQ=255)                                                                                     ttCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGt                  <  1:848474/139‑1 (MQ=255)                                                                                     ttCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATAGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGt                  <  2:339196/139‑1 (MQ=255)                                                                                          cacaTCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGAt                                                        <  2:158234/96‑1 (MQ=255)                                                                                          cacaTCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGAt                                                        >  1:158234/1‑96 (MQ=255)                                                                                             aTCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTg          >  2:1649639/1‑139 (MQ=255)                                                                                                  cGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATg                                                                >  2:647689/1‑80 (MQ=255)                                                                                                  cGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATg                                                                <  1:647689/80‑1 (MQ=255)                                                                                                  cGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACAcc     <  1:73026/139‑1 (MQ=255)                                                                                                     tATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTgc  <  2:1312100/139‑1 (MQ=255)                                                                                                     tATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTgc  <  1:267704/139‑1 (MQ=255)                                                                                                              |                                                                                                                                  GCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGGGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGC  >  NC_000913/405115‑405353

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