Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,080,039	(G)5→4	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	.	92.2%	211.2 / 12.8	51	coding (731/1509 nt)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base . (25/22);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCTTTGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTATGA  >  NC_000913/1079936‑1080130                                                                                                    |                                                                                                gCTTTGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTg                                                                                                 >  2:936977/1‑100 (MQ=255)   tttGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTggg                                                                                               >  1:50714/1‑100 (MQ=255)    ttGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTgggg                                                                                              <  2:829348/100‑1 (MQ=255)       gTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCt                                                                                           >  1:1201205/1‑100 (MQ=255)         gACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTggg                                                                                        <  1:1183705/100‑1 (MQ=255)          aCTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTGTGTTGCCATTATCTCCCTGATGGGTT‑GGGGCTGGGt                                                                                       >  2:428529/1‑100 (MQ=255)           cTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGtt                                                                                      >  1:63032/1‑100 (MQ=255)           cTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGtt                                                                                      >  1:238694/1‑100 (MQ=255)             cGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTAc                                                                                    >  1:1522167/1‑100 (MQ=255)                  gAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCgg                                                                               <  2:788119/100‑1 (MQ=255)                     gTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCa                                                                            <  1:493064/100‑1 (MQ=255)                     gTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCa                                                                            <  2:206308/100‑1 (MQ=255)                     gTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTCGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCa                                                                            <  2:874985/100‑1 (MQ=255)                       gATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGc                                                                          >  1:199881/1‑100 (MQ=255)                       gATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGc                                                                          >  2:160086/1‑100 (MQ=255)                       gATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGc                                                                          >  1:1272801/1‑100 (MQ=255)                         tCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCg                                                                        >  1:638643/1‑100 (MQ=255)                          caaacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGc                                                                       <  1:1402252/100‑1 (MQ=255)                           aaacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGc                                                                       <  2:148103/100‑1 (MQ=255)                            aacaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGcac                                                                     >  1:1432134/1‑100 (MQ=255)                             acaaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGcaca                                                                    <  1:221760/100‑1 (MQ=255)                              caaaAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTGCGGGCAGCCGCACAt                                                                   >  1:1418079/1‑100 (MQ=255)                                     gCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGgcg                                                            <  2:238693/100‑1 (MQ=255)                                     gCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGgcg                                                            <  2:63032/100‑1 (MQ=255)                                       aTCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTATTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGcgcg                                                          >  1:1110269/1‑100 (MQ=255)                                         cGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGt                                                         >  2:1131389/1‑100 (MQ=255)                                           aaaaCGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGttt                                                       <  2:542640/100‑1 (MQ=255)                                           aaaaCGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGtttt                                                      >  1:58340/1‑100 (MQ=255)                                                    atatGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGgcggcg                                             <  1:936985/100‑1 (MQ=255)                                                       tGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGAt                                          >  1:1247549/1‑100 (MQ=255)                                                        gCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGAtt                                         <  1:1517980/100‑1 (MQ=255)                                                        gCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGAtt                                         >  1:1178244/1‑100 (MQ=255)                                                         cTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATtc                                        >  1:802737/1‑100 (MQ=255)                                                         cTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATtc                                        >  1:1292631/1‑100 (MQ=255)                                                             aatgTCTGAACTTTGTTGCCATTATCTCACTGATGGGTT‑AGCGCTAGGTTACTTCGGGCAGTCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTca                                     <  2:1142947/96‑1 (MQ=255)                                                                       cTTTGTTGCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTg                          >  2:1048932/1‑100 (MQ=255)                                                                         ttgttgCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGt                         >  2:125898/1‑99 (MQ=255)                                                                           gttgCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCa                      >  2:1501174/1‑100 (MQ=255)                                                                           gttgCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCa                      >  1:628738/1‑100 (MQ=255)                                                                            ttgCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCAt                     >  1:586072/1‑100 (MQ=255)                                                                              gCCATTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTTACCACAGCATTGTCCATgc                   >  1:1209215/1‑100 (MQ=255)                                                                                 aTTATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGcgcg                <  2:687501/100‑1 (MQ=255)                                                                                   tATCTCACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGcgtc              >  2:396193/1‑100 (MQ=255)                                                                                     tctcACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGcgtcgt            <  1:897265/100‑1 (MQ=255)                                                                                     tctcACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGcgtcgt            >  2:260536/1‑100 (MQ=255)                                                                                        cACTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTAtt         <  2:254239/100‑1 (MQ=255)                                                                                         aCTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTa        <  2:1513108/100‑1 (MQ=255)                                                                                          cTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAg       <  1:260537/100‑1 (MQ=255)                                                                                          cTGATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAg       <  2:1036356/100‑1 (MQ=255)                                                                                            gATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTa     <  2:1429839/100‑1 (MQ=255)                                                                                            gATGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTa     <  1:1023235/100‑1 (MQ=255)                                                                                              tGGGTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTATg   >  2:1017546/1‑100 (MQ=255)                                                                                              tGGGTT‑GGGGCTGGGATACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTATg   <  1:1546033/100‑1 (MQ=255)                                                                                               gggTT‑GGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTCTGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTATGa  <  1:428530/100‑1 (MQ=255)                                                                                                    |                                                                                                GCTTTGGTGACTCGCTGGAAGTGATCAAACAAAAGAGCATCGAAAACGTTGATATGCTCAAAGGTCTGAACTTTGTTGCCATTATCTCACTGATGGGTTGGGGGCTGGGTTACTTCGGGCAGCCGCACATTCTGGCGCGTTTTATGGCGGCGGATTCTCACCACAGCATTGTCCATGCGCGTCGTATTAGTATGA  >  NC_000913/1079936‑1080130

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