Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	696,470	C→T	noncoding (35/75 nt)	glnX ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,470	0	C	T	97.9%	143.7 / ‑6.8	48	noncoding (35/75 nt)	glnX	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/21);  minor base A (0/1);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTT  >  NC_000913/696331‑696597                                                                                                                                            |                                                                                                                                gTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGa                                                                                                                              >  2:855751/1‑143 (MQ=255)    gaaTCGGACGTTGAAGGATAAGGCGTGGCAACATCGGATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGGGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATc                                                                                                                           <  1:552611/141‑1 (MQ=37)         ccaCATTGAAGGATAAGACGTGGCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTg                                                                                                                      <  1:759276/141‑1 (MQ=255)          gACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGc                                                                                                                     >  1:1030373/1‑143 (MQ=255)                    gATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCtt                                                                                                           <  1:68710/143‑1 (MQ=255)                      tAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGCAATGCCGGAATTAGAATCCGGTGCCTTACAACTTgg                                                                                                         >  1:1092851/1‑143 (MQ=255)                        aGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCg                                                                                                       <  2:135395/143‑1 (MQ=255)                          aCGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGAt                                                                                                     >  1:716922/1‑143 (MQ=255)                          aCGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGAt                                                                                                     >  2:503216/1‑143 (MQ=255)                              gtCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTGTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGGATTCGAATCCGGGGCGTTGCCGCCTGGCGAGAccc                                                                                                 >  1:461682/1‑143 (MQ=37)                               tCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATAcccc                                                                                                >  1:819041/1‑143 (MQ=255)                               tCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATAcccc                                                                                                >  2:49216/1‑143 (MQ=255)                                     tCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGCTACGATGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTgg                                                                                          <  1:1106216/143‑1 (MQ=255)                                        cATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGACTTAGACTCCGGTGCCTTACCGCTTGGCGATACCCCAACTGAGTg                                                                                       >  2:858729/1‑143 (MQ=255)                                          ttCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCa                                                                                     <  2:922640/143‑1 (MQ=255)                                           tCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAGGGTGCAc                                                                                    >  2:751458/1‑143 (MQ=255)                                               cATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGGCTTACCGCTTGGCGATACCCCAAGTGGGGGCACttgc                                                                                >  1:214114/1‑141 (MQ=255)                                                aTTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCACACAGGCTGCACTTACa                                                                               >  2:236670/1‑143 (MQ=25)                                                 ttGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACaa                                                                              >  1:813488/1‑143 (MQ=255)                                                 tggaatgagCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACaa                                                                              <  1:810032/141‑1 (MQ=255)                                                    agtgaaCGCAGAAGAGGAAGACGCTCGCCGAAGCGGGCTGTTGTACTGTGGCCACAGTACGAGGACCCGAACCGCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGt                                                                           <  2:259251/141‑1 (MQ=18)                                                    aatgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAGGCTGGAGGTACAAGGt                                                                           >  2:581335/1‑143 (MQ=21)                                                      tgaaCGCAGAAAGGCAAACAGCCCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTaa                                                                         <  2:971005/143‑1 (MQ=255)                                                           gCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTc                                                                    <  2:469259/143‑1 (MQ=255)                                                             agagAAGGAAAAAGGCCGCCGAAGCGAGCTTTGGTAATGTGGCTGGGGTACGAGGGTTCGAACCTCGGAATGCCGGAATTAGGATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCtt                                                                  <  2:406908/143‑1 (MQ=21)                                                             aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCtt                                                                  <  1:819231/143‑1 (MQ=255)                                                             aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGGGCACTCACACGGTAAGCGCCtt                                                                  >  1:721297/1‑143 (MQ=21)                                                             aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACACCAACAAGTTGGACTTACAAGGTAcgcgcatt                                                                  >  2:1166258/1‑139 (MQ=18)                                                                   gCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTAGGTGCACTTACAAGGTAAGCGCCTTGAATaa                                                            >  2:563738/1‑143 (MQ=255)                                                                    cAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATaaa                                                           <  1:49216/143‑1 (MQ=255)                                                                       ggaGCTCGCCGAAGCGAGCTTTTTGAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGACATAAAAATCTGGAGTTATCAATGGTGTCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTg                                                        <  2:878056/141‑1 (MQ=21)                                                                             cGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTggg                                                  <  1:855751/143‑1 (MQ=255)                                                                             cGCCGAAGCGAGCTTGTGTAATGAGGCTGGGGTACGAGGATCCGAACCTCGGAATGCCGGAATTAGACTCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTggg                                                  <  1:1198006/143‑1 (MQ=25)                                                                              gCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCCCTTACAAGGTAAGCGTCTTGAATAAATTGGCTgggg                                                 >  1:492834/1‑143 (MQ=255)                                                                                      gAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGa                                         >  1:884339/1‑143 (MQ=255)                                                                                        gCTTTTTTAATTTGGCTGGGGTACGAGGATCCGAACCTCGTAATGCCGGAATTAGACTCCGGTGCCTTACCGCTTGGTGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGAtt                                       <  2:451235/143‑1 (MQ=25)                                                                                         cTTTTTGAATGTGGGTGGGGTAGGAGGAGTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTCCCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTc                                      <  1:984655/143‑1 (MQ=25)                                                                                          ttttttAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAACCGGCTGGGGTACGAGGATTCg                                     >  1:457558/1‑143 (MQ=255)                                                                                          ttgtttAATGTGGCTGGGGTACGAGCCTTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTCAATAAATTGGCTGGGGTACGAGGATTCg                                     <  2:1013230/140‑1 (MQ=25)                                                                                           tttttAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAAGTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTAACAAGGGAAGCGGCGTGACTAAACTTGCAGGGGGACGAGGATTCGa                                    >  1:558367/1‑143 (MQ=18)                                                                                            ttttAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGaa                                   >  1:613998/1‑143 (MQ=255)                                                                                                  tgtgGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTGCAAGGTAAGCGGCTTGACTAAATTGGCTGGGGTACGTGGCTTCGaccggcgg                             >  1:1128464/1‑136 (MQ=21)                                                                                                    tgGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGaa                           <  1:727706/143‑1 (MQ=255)                                                                                                      gCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGCATg                         >  1:135789/1‑143 (MQ=255)                                                                                                               cGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTACGCGTCTTGAATAAATTGGCTGGGGTACGCGGATTCGAACCTCGGAATGCCGGAATCa                >  1:141722/1‑143 (MQ=255)                                                                                                                        gAACCTCGGAATGCCGGAATTAGAAGCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt       <  1:1045469/143‑1 (MQ=255)                                                                                                                        gAACCTCGGAATGCCGGAAGTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt       <  1:755148/143‑1 (MQ=255)                                                                                                                             tCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGTGTACGAGGATTCGAACCTCGGAATGGCGGAATCAGAATCGGGTGACtt  >  1:472381/1‑143 (MQ=37)                                                                                                                                            |                                                                                                                                GTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTT  >  NC_000913/696331‑696597

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