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%	106.3 / NA	34	P141S (CCA→TCA)	proC	pyrroline‑5‑carboxylate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/19);  total (15/19)

AACATAAATACGTAGGCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGGGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAA  >  NC_000913/405100‑405358                                                                                                                             |                                                                                                                                       aaCATAAATACGTAGGCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGt                                                                                                                          <  1:953270/139‑1 (MQ=255)          aCGTAGGCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTa                                                                                                                 <  1:941348/139‑1 (MQ=255)              aGGCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGa                                                                                                             >  1:207440/1‑139 (MQ=255)                 cTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGt                                                                                                          <  1:244752/139‑1 (MQ=255)                    gCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCAt                                                                                                       >  2:819925/1‑139 (MQ=255)                      gaagaaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATcc                                                                                                     <  1:110097/139‑1 (MQ=255)                      gaagaaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATcc                                                                                                     <  2:322809/139‑1 (MQ=255)                      gaagaaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATcc                                                                                                     <  1:842629/139‑1 (MQ=255)                          aaCCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGc                                                                                                 >  1:643632/1‑139 (MQ=255)                               cTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTaa                                                                                            >  2:141721/1‑139 (MQ=255)                                        accaccGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCgg                                                                                   >  2:59076/1‑139 (MQ=255)                                          caccGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGa                                                                                 <  2:643632/139‑1 (MQ=255)                                          caccGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGa                                                                                 <  2:301622/139‑1 (MQ=255)                                           accGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTa                                                                                             <  2:182454/126‑1 (MQ=255)                                           accGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTa                                                                                             >  1:182454/1‑126 (MQ=255)                                               ggTGGATCATCGGCGTAGCTATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTAGCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGtt                                                                            <  1:1015606/139‑1 (MQ=255)                                                   gATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGc                                                                        >  2:127822/1‑139 (MQ=255)                                                          ggCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGgcgc                                                                 <  1:59076/139‑1 (MQ=255)                                                          ggCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGgcgc                                                                 <  1:922841/139‑1 (MQ=255)                                                           gCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGcgcg                                                                >  2:563515/1‑139 (MQ=255)                                                             tCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGa                                                              <  1:127822/139‑1 (MQ=255)                                                              cAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGAt                                                             >  1:604235/1‑139 (MQ=255)                                                                 cAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAAt                                                          >  2:521169/1‑139 (MQ=255)                                                                     tACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTc                                                      <  1:970669/139‑1 (MQ=255)                                                                                  ccAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCa                                         >  2:865926/1‑139 (MQ=255)                                                                                            ggAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGa                                                                                 <  1:276521/89‑1 (MQ=255)                                                                                            ggAAAATATTCAGCACATCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGa                                                                                 >  2:276521/1‑89 (MQ=255)                                                                                                           caTCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGt                <  2:895388/139‑1 (MQ=255)                                                                                                            aTCAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTg               >  2:1000730/1‑139 (MQ=255)                                                                                                              cAGCGGTATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGac             >  2:841294/1‑139 (MQ=255)                                                                                                                    tATCTTCTGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTgc       <  2:797686/139‑1 (MQ=255)                                                                                                                         tctGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGccc                                          <  2:752096/99‑1 (MQ=255)                                                                                                                         tctGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGccc                                          >  1:752096/1‑99 (MQ=255)                                                                                                                         tctGAGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCaa  <  1:521169/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                       AACATAAATACGTAGGCTGGCGAAGAACCGCTCACACCGACCACCGGGTGGATCATCGGCTCAGCAATTACTTCCGCTTCGCCAAAGCAGCGGAAAATATTCAGCACATCAGCGGTATCTTCTGGGGTTACCAGCGCGTTTGGCGTTACGGAGGTCATCCCGGCATTAACCAGTGCGGGAGTGTTCGGCATGGCGCGGATAATTTTCCGGTCATGGCCCAGCGCGCGGGCAAGCTGGTCGAGCGTGACACCTGCAGCAA  >  NC_000913/405100‑405358

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