Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,331,634	C→T	A378V (GCG→GTG)	proP →	proline/glycine betaine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,331,634	0	C	T	100.0%	190.4 / NA	59	A378V (GCG→GTG)	proP	proline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/26);  total (33/26)

GCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCA  >  NC_000913/4331501‑4331745                                                                                                                                      |                                                                                                                gCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGt                                                                                                            >  2:1343323/1‑139 (MQ=255)   ccTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTc                                                                                                          >  2:1285289/1‑139 (MQ=255)    cTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                         <  2:704914/139‑1 (MQ=255)      gTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCg                                                                                                       <  2:59911/139‑1 (MQ=255)                 cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                        >  2:356687/1‑127 (MQ=255)                 cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                        <  1:356687/127‑1 (MQ=255)                 cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCg                                                                                            >  1:1058618/1‑139 (MQ=255)                                     aTTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcgg                                                                        >  2:210008/1‑139 (MQ=255)                                         aCAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                    <  2:570527/139‑1 (MQ=255)                                               aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                         <  1:383718/96‑1 (MQ=255)                                               aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                         >  2:383718/1‑96 (MQ=255)                                               aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAatat                                                              >  1:76981/1‑139 (MQ=255)                                                cGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGATGGCGGCGGCATTTAATAtt                                                             >  1:1227411/1‑139 (MQ=255)                                                    aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGAt                                                                                                              <  1:479289/86‑1 (MQ=255)                                                    aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGAt                                                                                                              >  2:479289/1‑86 (MQ=255)                                                    aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                      <  1:423776/126‑1 (MQ=255)                                                    aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                      >  2:423776/1‑126 (MQ=255)                                                    aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCgg                                                         <  1:505423/139‑1 (MQ=255)                                                     tCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                           <  2:356774/120‑1 (MQ=255)                                                     tCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                           >  1:356774/1‑120 (MQ=255)                                                     tCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGt                                                        >  2:19617/1‑139 (MQ=255)                                                           tGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGt                                                  >  2:409030/1‑139 (MQ=255)                                                             aTTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTg                                                <  2:76981/139‑1 (MQ=255)                                                                  tGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTg                                                    >  2:653921/1‑130 (MQ=255)                                                                  tGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTg                                                    <  1:653921/130‑1 (MQ=255)                                                                                 gCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTg                            >  2:882646/1‑139 (MQ=255)                                                                                  cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                         >  2:549295/1‑93 (MQ=255)                                                                                  cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                         <  1:549295/93‑1 (MQ=255)                                                                                  cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTgg                           >  2:591941/1‑139 (MQ=255)                                                                                         tGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcgg                                                                        <  2:551723/87‑1 (MQ=255)                                                                                         tGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcgg                                                                        >  1:551723/1‑87 (MQ=255)                                                                                           aTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctggc                  >  2:976530/1‑139 (MQ=255)                                                                                            tCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctggct                 <  1:976530/139‑1 (MQ=255)                                                                                               ttAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGt              >  1:1017341/1‑139 (MQ=255)                                                                                                      cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc       >  2:237269/1‑139 (MQ=255)                                                                                                      cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc       >  2:742568/1‑139 (MQ=255)                                                                                                       tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCTAAAGCt      >  1:1181242/1‑139 (MQ=255)                                                                                                       tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt      >  1:251227/1‑139 (MQ=255)                                                                                                        ttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTc     <  2:1058618/139‑1 (MQ=255)                                                                                                         tACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg    <  1:742568/139‑1 (MQ=255)                                                                                                           cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa  <  1:237269/139‑1 (MQ=255)                                                                                                              gCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                   <  1:1027587/39‑1 (MQ=255)                                                                                                              gCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                   >  2:1027587/1‑39 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                         <  1:32237/63‑1 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                         >  2:32237/1‑63 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                    <  1:311003/68‑1 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                    >  2:311003/1‑68 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGt                                           <  2:935762/93‑1 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGc   >  2:524107/1‑133 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGc   <  1:524107/133‑1 (MQ=255)                                                                                                                gTTATGGCTTCTACCTTGCCAGTGATGTTACCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGt                                           >  1:935762/1‑93 (MQ=255)                                                                                                                         tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGcc                                   <  1:590709/92‑1 (MQ=255)                                                                                                                         tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGcc                                   >  2:590709/1‑92 (MQ=255)                                                                                                                          cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                         >  1:1015012/1‑53 (MQ=255)                                                                                                                          cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                         <  2:1015012/53‑1 (MQ=255)                                                                                                                          cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTgg                           >  2:608938/1‑99 (MQ=255)                                                                                                                          cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTgg                           <  1:608938/99‑1 (MQ=255)                                                                                                                                tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                           >  1:249200/1‑45 (MQ=255)                                                                                                                                tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                           <  2:249200/45‑1 (MQ=255)                                                                                                                                      |                                                                                                                GCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCA  >  NC_000913/4331501‑4331745

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