Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,643,899	G→A	G428G (GGC→GGT)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,899	0	G	A	100.0%	163.5 / NA	51	G428G (GGC→GGT)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/28);  total (23/28)

GTGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACA  >  NC_000913/3643764‑3644034                                                                                                                                        |                                                                                                                                        gTGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCt                                                                                                                                      <  2:1359339/139‑1 (MQ=255)     ccGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGc                                                                                                                                  <  1:1533939/139‑1 (MQ=255)          cTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTca                                                                                                                             >  1:1585502/1‑139 (MQ=255)              tGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                                 <  2:858907/131‑1 (MQ=255)              tGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                                 >  1:858907/1‑131 (MQ=255)                 aaGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCatccatcc                                                                                                                      <  2:1585502/139‑1 (MQ=255)                           tCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCccgcc                                                                                                            >  2:473175/1‑139 (MQ=255)                                   tcgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGt                                                                                                    >  2:527327/1‑139 (MQ=255)                                      tgagtgaACAGCGCCGGTTTACCATTCACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                                 <  2:12013/139‑1 (MQ=255)                                                gcgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                       <  2:1409348/139‑1 (MQ=255)                                                 cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATAc                                                                                      >  1:176097/1‑139 (MQ=255)                                                                 aCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             >  1:195207/1‑116 (MQ=255)                                                                 aCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             <  2:195207/116‑1 (MQ=255)                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGc                                                                                                                                  >  1:1155800/1‑75 (MQ=255)                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGc                                                                                                                                  <  2:1155800/75‑1 (MQ=255)                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCgcgc                                                                                                        <  2:444190/101‑1 (MQ=255)                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCgcgc                                                                                                        >  1:444190/1‑101 (MQ=255)                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTAcc                                                                  >  2:1127408/1‑139 (MQ=255)                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTAcc                                                                  >  2:644562/1‑139 (MQ=255)                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTAcc                                                                  >  1:941762/1‑139 (MQ=255)                                                                      gCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcg                                                                 >  1:1034268/1‑139 (MQ=255)                                                                        ggTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcgcg                                                               <  2:1034268/139‑1 (MQ=255)                                                                        ggTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcgcg                                                               <  2:632703/139‑1 (MQ=255)                                                                          ttGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCa                                                             >  2:206071/1‑139 (MQ=255)                                                                               gTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCg                                                                                <  1:905078/115‑1 (MQ=255)                                                                               gTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCg                                                                                >  2:905078/1‑115 (MQ=255)                                                                                   cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                    <  1:206071/139‑1 (MQ=255)                                                                                   cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                    <  1:280017/139‑1 (MQ=255)                                                                                   cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                    <  2:615552/139‑1 (MQ=255)                                                                                     aGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAttt                                                  >  1:599464/1‑139 (MQ=255)                                                                                        cAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCg                                               <  1:59339/139‑1 (MQ=255)                                                                                          aaTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTc                                             >  2:1360304/1‑139 (MQ=255)                                                                                           aTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCa                                            <  1:527327/139‑1 (MQ=255)                                                                                             aCGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                          <  2:42853/139‑1 (MQ=255)                                                                                                cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCcgcg                                                                                                         >  2:1200656/1‑73 (MQ=255)                                                                                                cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCcgcg                                                                                                         <  1:1200656/73‑1 (MQ=255)                                                                                                       tttttGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCgaag                                <  1:1283755/139‑1 (MQ=255)                                                                                                                 gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcg                                                                 <  2:695642/96‑1 (MQ=255)                                                                                                                 gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcg                                                                 >  1:695642/1‑96 (MQ=255)                                                                                                                       tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAttt                                                  >  2:968812/1‑105 (MQ=255)                                                                                                                       tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAttt                                                  <  1:968812/105‑1 (MQ=255)                                                                                                                       tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                    <  2:640668/103‑1 (MQ=255)                                                                                                                       tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                    >  1:640668/1‑103 (MQ=255)                                                                                                                         aGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGc              <  2:176097/139‑1 (MQ=255)                                                                                                                         aGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGc              <  1:644562/139‑1 (MQ=255)                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             >  2:637594/1‑55 (MQ=255)                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             <  1:637594/55‑1 (MQ=255)                                                                                                                                 gCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATc      <  1:1360304/139‑1 (MQ=255)                                                                                                                                     cTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGa                                                                         <  2:78933/68‑1 (MQ=255)                                                                                                                                     cTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGa                                                                         >  1:78933/1‑68 (MQ=255)                                                                                                                                     cTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCCTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACa  >  1:1509491/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        GTGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACA  >  NC_000913/3643764‑3644034

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