Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,080,039	(G)5→4	100%	coding (735/1509 nt)	putP →	proline:sodium symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,080,035	0	G	.	56.2%	8.4 / 143.3	64	coding (731/1509 nt)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base G (18/10);  new base . (12/24);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

GCTTTGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTAT  >  NC_000913/1079936‑1080128                                                                                                    |                                                                                              gCTTTGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTg                                                                                               <  2:1195320/100‑1 (MQ=255)    ttGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTgggg                                                                                            <  1:1032949/100‑1 (MQ=255)     tGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTggggc                                                                                           <  2:1604768/100‑2 (MQ=255)     tGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTggggc                                                                                           <  2:1168209/100‑2 (MQ=255)     tGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTgggag                                                                                           >  1:1392935/1‑98 (MQ=255)           cTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGtt                                                                                    <  1:1153816/100‑1 (MQ=255)           cTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGtt                                                                                    >  1:976674/1‑100 (MQ=255)           cTCGCTGGAAGTGATCAAACAAAACAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGGTGGGGGCTggg                                                                                      <  2:350022/99‑1 (MQ=255)              gCTGGAAGTGATCAAAGAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTAc                                                                                  <  1:441230/99‑1 (MQ=255)              gCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGttgc                                                                                  >  1:2235614/1‑98 (MQ=255)              gCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTAc                                                                                  <  2:1591658/100‑1 (MQ=255)              gCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACt                                                                                 <  1:651107/100‑1 (MQ=255)                 ggAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCg                                                                              >  1:479383/1‑100 (MQ=255)                 ggAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCg                                                                              <  1:2105068/100‑1 (MQ=255)                  gAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCgg                                                                             <  1:1367526/100‑1 (MQ=255)                       gATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAg                                                                         >  1:987980/1‑100 (MQ=255)                         tCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGATACTACGGGCAGcc                                                                       >  2:233022/1‑100 (MQ=255)                         tCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCg                                                                      >  1:1508303/1‑100 (MQ=255)                         tCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCACACTGCTGGGTTGGGGGCTGGGTTACTTCGGGCAGcc                                                                       >  2:1613654/1‑100 (MQ=255)                          caaacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCg                                                                      >  2:364166/1‑100 (MQ=255)                          caaacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCg                                                                      <  1:185106/100‑1 (MQ=255)                           aaacaaaAGAGCATCGAAAACGTTGATATGGTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGca                                                                    <  2:2266001/100‑1 (MQ=255)                           aaacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGc                                                                     >  2:893571/1‑100 (MQ=255)                            aacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGca                                                                    >  2:1694096/1‑100 (MQ=255)                            aacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGaac                                                                   <  2:795239/100‑1 (MQ=255)                              caaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACAt                                                                 >  2:1548058/1‑100 (MQ=255)                               aaaaGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACAtt                                                                >  2:2322111/1‑100 (MQ=255)                                aaaGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTc                                                               >  1:2301339/1‑100 (MQ=255)                                     gCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGTCGCACATTCTGgc                                                           >  1:781844/1‑100 (MQ=255)                                      cATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGgc                                                           >  1:526809/1‑99 (MQ=255)                                         cGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGcgcg                                                        >  2:1813630/1‑99 (MQ=255)                                         cGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGtt                                                      <  2:1324280/100‑1 (MQ=255)                                            aaaCGTTGATATGCTCAAAGGTCTGAACTTGGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGGCGCACATTCTGGCGCGTTTTa                                                   <  2:1344178/100‑1 (MQ=255)                                              aCGTTTATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTAt                                                  <  1:309748/99‑1 (MQ=255)                                               cGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTTTGTTACTTCTGGCAGCCGCACATTCTGGCGCGTTTTATg                                                 <  2:1956367/100‑1 (MQ=255)                                               cGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATgg                                                <  1:1808396/100‑1 (MQ=255)                                                    atatGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGgcggcg                                           >  1:338138/1‑100 (MQ=255)                                                          tCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGcggcgg                                          <  2:2180735/96‑1 (MQ=255)                                                          tCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATtct                                     <  2:976663/100‑1 (MQ=255)                                                           cAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATtctc                                    <  2:1216802/100‑1 (MQ=255)                                                            aaaGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATtctc                                    <  2:146357/100‑1 (MQ=255)                                                            aaaGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCGGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATtctc                                    >  2:1324374/1‑100 (MQ=255)                                                              aGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCTCGTTTTATGGCGGCGGATTCTcacc                                 <  1:1114477/100‑1 (MQ=255)                                                              aGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTcacc                                 <  2:1192685/100‑1 (MQ=255)                                                              aGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTcacc                                 <  1:757742/100‑1 (MQ=255)                                                                gTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTcaccac                               <  2:1125062/100‑1 (MQ=255)                                                                  cTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACcaca                              >  2:425267/1‑100 (MQ=255)                                                                    gAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGc                            >  2:16028/1‑100 (MQ=255)                                                                     aaCTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCAt                          <  2:2144582/100‑1 (MQ=255)                                                                     aaCTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCAt                          <  2:177077/100‑1 (MQ=255)                                                                          tgttgCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCATCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTc                      <  2:1048497/99‑1 (MQ=255)                                                                             tgCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCAt                   >  1:2195385/1‑100 (MQ=255)                                                                               ccATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATgcg                >  1:546517/1‑100 (MQ=255)                                                                                cATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATgcg                <  1:893582/100‑1 (MQ=255)                                                                                   tATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGcgtc            <  1:2380170/100‑1 (MQ=255)                                                                                    aTCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGcgtcg           >  1:1918194/1‑100 (MQ=255)                                                                                     tctcACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGcgtcg           >  2:1508925/1‑100 (MQ=255)                                                                                     tctcACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGcgtcgt          >  2:535223/1‑100 (MQ=255)                                                                                     actcACTGATGGTATGGGGGCTGGGATACTTCGGGCAGCCGCACGTTCTGGTGCGTATTATGGCGGCGGATTCTCACCACAGCATTGACCGTGCGcgtcg           >  1:2393002/2‑100 (MQ=255)                                                                                      ctcACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTa         <  2:479380/100‑1 (MQ=255)                                                                                       tcACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTa         <  1:1508942/100‑1 (MQ=255)                                                                                        cACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTAt        >  2:602206/1‑100 (MQ=255)                                                                                          cTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAg     <  1:143620/100‑1 (MQ=255)                                                                                          cTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGAGGATTATCACCACAGCATTGTCCATGAGCGTAGTATTAg     >  1:2293931/1‑100 (MQ=255)                                                                                           tGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAg     <  2:987969/100‑1 (MQ=255)                                                                                           tGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGt    >  1:2145757/1‑100 (MQ=255)                                                                                            gATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGt    >  2:411792/1‑100 (MQ=255)                                                                                            gATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGt    <  1:381141/100‑1 (MQ=255)                                                                                             aTGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTAt  <  2:926887/100‑1 (MQ=255)                                                                                                    |                                                                                              GCTTTGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTAT  >  NC_000913/1079936‑1080128

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