Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,959,758	A→G	L19P (CTG→CCG)	ppiC ←	peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,959,758	0	A	G	100.0%	203.2 / NA	53	L19P (CTG→CCG)	ppiC	peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/28);  total (25/28)

TTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTC  >  NC_000913/3959623‑3959887                                                                                                                                        |                                                                                                                                  ttATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCggag                                                                                                                                >  2:38934/1‑139 (MQ=255)   aTCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGcc                                                                                                                              <  1:944486/139‑1 (MQ=255)          gCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCACGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTtct                                                                                                                       <  2:728055/139‑1 (MQ=255)           ccGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTtctc                                                                                                                      <  1:38934/139‑1 (MQ=255)                   aTCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTa                                                                                                              <  2:962478/139‑1 (MQ=255)                   aTCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTa                                                                                                              <  1:601264/139‑1 (MQ=255)                   aTCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTa                                                                                                              <  1:415637/139‑1 (MQ=255)                    tCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTa                                                                                                              >  2:294556/1‑138 (MQ=255)                    tCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTa                                                                                                              <  1:294556/138‑1 (MQ=255)                    tCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTAc                                                                                                             >  1:686026/1‑139 (MQ=255)                            tGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGatat                                                                                                     >  1:355535/1‑139 (MQ=255)                              gCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGc                                                                                                   >  2:238574/1‑139 (MQ=255)                                  aaTTCACCTCAAGCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCCGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTCCTTTTACAAGGATATGCCGtg                                                                                               >  1:209086/1‑139 (MQ=255)                                   aTTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGtgc                                                                                              >  1:94679/1‑139 (MQ=255)                                       aCCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTa                                                                                                              >  2:482603/1‑119 (MQ=255)                                       aCCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTa                                                                                                              <  1:482603/119‑1 (MQ=255)                                       aCCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGtgctgct                                                                                          >  1:424501/1‑139 (MQ=255)                                             aTCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGtt                                                                                    <  1:141002/139‑1 (MQ=255)                                                  cgccgcGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGCAATCGGCCCCTTTCTTAATCTGCTCGAGAAGATCCGGCGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGc                                                                               >  1:249508/1‑139 (MQ=255)                                                  cgccgcGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTAc                                                                                                             <  1:285137/109‑1 (MQ=255)                                                  cgccgcGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTAc                                                                                                             >  2:285137/1‑109 (MQ=255)                                                    ccgcGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCa                                                                             <  2:424501/139‑1 (MQ=255)                                                    ccgcGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCa                                                                             <  2:383685/139‑1 (MQ=255)                                                     cgcgTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCAt                                                                            >  1:133857/1‑139 (MQ=255)                                                      gcgTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGc                                                                                                                               >  1:1077813/1‑87 (MQ=255)                                                      gcgTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGc                                                                                                                               <  2:1077813/87‑1 (MQ=255)                                                              cTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGcc                                                                   <  2:94679/139‑1 (MQ=255)                                                              cTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGcc                                                                   <  2:1031612/139‑1 (MQ=255)                                                                gATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCtt                                                                 <  1:255784/139‑1 (MQ=255)                                                                  tGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGa                                                               <  1:238574/139‑1 (MQ=255)                                                                      cAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGa                                                           >  2:1087117/1‑139 (MQ=255)                                                                                   tcttcGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACa                                                                                                            <  1:649053/77‑1 (MQ=255)                                                                                   tcttcGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACa                                                                                                            >  2:649053/1‑77 (MQ=255)                                                                                   tcttcGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAggc                                              >  1:470513/1‑139 (MQ=255)                                                                                     ttcGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAggcgg                                            <  2:249508/139‑1 (MQ=255)                                                                                       cGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGcc                                                                                                                              <  1:502658/55‑1 (MQ=255)                                                                                       cGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGcc                                                                                                                              >  2:502658/1‑55 (MQ=255)                                                                                           aGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGAttt                                      >  2:40435/1‑139 (MQ=255)                                                                                             cTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTtc                                                                                                                        >  1:129306/1‑55 (MQ=255)                                                                                             cTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTtc                                                                                                                        <  2:129306/55‑1 (MQ=255)                                                                                                gCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCa                                 <  1:351594/139‑1 (MQ=255)                                                                                                    aaaTCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGAt                             <  1:126190/139‑1 (MQ=255)                                                                                                    aaaTCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGAt                             <  2:355535/139‑1 (MQ=255)                                                                                                              cGTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGc                   <  2:336106/139‑1 (MQ=255)                                                                                                               gTTCTTAATCTGCTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCt                  >  2:229104/1‑139 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCt                                                                  >  2:201592/1‑79 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCt                                                                  <  1:201592/79‑1 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAggcgg                                            >  1:163876/1‑101 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAggcgg                                            <  2:163876/101‑1 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTgtc      >  1:297145/1‑139 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTgtc      >  2:600234/1‑139 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTgtc      >  1:493277/1‑139 (MQ=255)                                                                                                                           cTCGAGAAGATCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTgtc      >  1:672559/1‑139 (MQ=255)                                                                                                                               agaagaTCCGGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTgtcagtc  <  2:508499/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  TTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTCTTCTTTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTC  >  NC_000913/3959623‑3959887

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: ‑