Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	674,768	G→T	R6S (CGC→AGC)	leuS ←	leucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	674,768	0	G	T	97.8%	118.7 / ‑3.8	46	R6S (CGC→AGC)	leuS	leucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (23/22);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAAT  >  NC_000913/674633‑674906                                                                                                                                        |                                                                                                                                           ccATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCCATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGc                                                                                                                                     >  2:18202/1‑143 (MQ=255) ccATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGATCTCGTATTCGGCAACTTCAGGTGTGTCTATTGCAGCTCAATGAAGCTGCAATACAGATTCTATCTCTTCCCGGCTGTCTTGc                                                                                                                                     >  2:829784/1‑143 (MQ=255)        tAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCa                                                                                                                              >  1:174146/1‑143 (MQ=255)           tCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATgg                                                                                                                           <  1:399817/143‑1 (MQ=255)            cGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCGGTATTGCTCTTGCATGGc                                                                                                                          <  2:766695/143‑1 (MQ=255)             gACCAGAAGGATAGTGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTTTTCGGTTATTTCAAATGTGCGCTTGTCATCGGGATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCa                                                                                                                         <  1:270713/143‑1 (MQ=255)              aCCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCCAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGTGTTCTTTCTCTTCCGGGTTGTTTTGCTTTTACATGACAg                                                                                                                        >  1:254659/1‑143 (MQ=255)               ccAGATGGATAGCGATGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGGGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCCATCTATTACGGGTTTTATTGCTCTTGCATGGCAGc                                                                                                                       <  2:744734/143‑1 (MQ=255)                      gATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTCCTTCAAATGTGCGCTTCTCCTCCCATTTAACCTGTACTTTAAATTCTATCTCAACCAAACTGAATTGCTCTTGCATGGCACCCAGTCGt                                                                                                                >  1:414013/1‑143 (MQ=255)                          gggAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGAGCTGTATTGCTCTTGCATGGCAGCCAGTgggtcac                                                                                                            >  2:203500/1‑138 (MQ=255)                              tgCATTGACAGGCAGTAATGGTTGTCGTAGCGGTCGGCGTCGGGTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGtttt                                                                                                        <  1:271459/142‑1 (MQ=255)                                      cAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACgg                                                                                                <  2:688572/143‑1 (MQ=255)                                      cAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCAGCACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACgg                                                                                                <  1:313565/143‑1 (MQ=255)                                           aGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGTCTAcc                                                                                           >  2:133248/1‑142 (MQ=255)                                             tAATGCTTCTCTTTGCTCTCGTCTTCGGGTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACaaa                                                                                         <  1:362668/143‑1 (MQ=255)                                             tAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACaaa                                                                                         <  1:362788/143‑1 (MQ=255)                                                 aCTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCGAGTGGGCCTGTTTTCAGTACGGCTGCacgagta                                                                                     >  2:494442/1‑137 (MQ=255)                                                    tctcTTTGCTGTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCAGCCCAATGAAGCTGTACTTTGGATTCTAGCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCg                                                                                  <  2:362782/143‑1 (MQ=255)                                                       ctTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTg                                                                               <  2:242241/143‑1 (MQ=255)                                                         ttGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAg                                                                             <  2:174146/143‑1 (MQ=255)                                                             tctcGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGTCTACAAATGGAGCGTCGAggtgg                                                                         >  2:650771/1‑143 (MQ=255)                                                                cGTCTTCGGTTACTTCAAGTGTGCTCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGttt                                                                      <  2:667242/143‑1 (MQ=255)                                                                  tCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTCTTGCTGTTGCATGGCAGCCCGTGGTCCTGTTTTCAATACGGCTACAAATGTAGAGTTGAGGTGGtttgt                                                                    >  1:327427/1‑141 (MQ=255)                                                                        gTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAACTATAGCGTTGAGGTGGTTTTTCAGATc                                                              >  1:578300/1‑143 (MQ=255)                                                                        gTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGttttttttttg                                                              >  2:440803/1‑137 (MQ=255)                                                                        gTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACCGCTACAAATGTAGCGTCGGGGTGGGTTTTccgctc                                                              >  2:459690/1‑138 (MQ=255)                                                                        gTTACTTCAAATGTGCGCTTCTCAGCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGGTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATc                                                              <  1:43262/143‑1 (MQ=255)                                                                          tACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTGGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCg                                                            <  1:564343/143‑1 (MQ=255)                                                                           agttCAAATGTGGGCTTATCGTCCCATGGAGGCTCTACTTTGGGTTCTATCTCTTGCGTGCTGTATTGCTCTGGGATGGCAGTCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCgc                                                           <  2:277691/141‑1 (MQ=255)                                                                             ttCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTCCACATGTCGCGTTGTGGTGGGTTTTCAGATCCgcat                                                         >  1:635519/1‑143 (MQ=255)                                                                             ttCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAACAGTAGCGTAGCGGTGGTTTTTCAGCTCCccat                                                         >  2:582973/1‑143 (MQ=255)                                                                             ttCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCgcat                                                         >  1:415131/1‑143 (MQ=255)                                                                             ttCAAATGTGCCGTTCTTAGCCGGTTTGAGCTGCACGTTGGATTCTCTCTCTTCCGGCCTATATTGCTCTTGCATTGCAGCGAGTAGTCGTGTTGTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCgcat                                                         <  2:845567/143‑1 (MQ=255)                                                                                   tgttCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTAACGCGTCCGGGCTGTAGTGCACTTGCAGGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCgcatagcata                                                   <  1:135005/143‑1 (MQ=255)                                                                                     tgCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTTTTGCATGGCTGTCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCatagcatagc                                                 <  1:94390/143‑1 (MQ=255)                                                                                              cATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCCGTGGTCCTGTTTTCAATACGGCTCCTCATGTAGCGTTGGCGTGGTTTTTCCAATCCACATAGCATATCCCAAACGTc                                        >  1:169008/1‑143 (MQ=255)                                                                                                tCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCCTGGCAGCCAGTGGTCCTGTTTTCGATACGGCTCCAAATGTAGCGGTGAGGTGGATGTTAAGATCTGAAAAGGATAGACTaaataagag                                      >  1:237073/1‑137 (MQ=255)                                                                                                  ccAATGAAGCTGTACTTTGGATTCTATCTCTTCCGAGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAACGTAGCGTTGAGGTGGCTTTTCAGATCCGCATAGCATAGCCCAAACGTCcgaa                                    >  2:162995/1‑141 (MQ=255)                                                                                                        aaGCTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGAGGGCCTGTGTACAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCaaaa                              <  1:284089/143‑1 (MQ=255)                                                                                                           cTGTACTTTGGATTCTATCTCTTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGAGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGGCCGCGTCAAAACAg                           >  2:578006/1‑143 (MQ=255)                                                                                                                  ttGGGTTCTATCTCTGCCGGGCTGTATTGCGCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTGGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTcg                    <  2:25241/143‑1 (MQ=255)                                                                                                                  ttGGATTCTATCTCTTCCGGCCTGTATTGCCCTTGCATGGCAGCCAGTGGCCCTGCTTTCAACAACGCTACCAATGTCACGTTGAGGAGGTTTTTCAGATCCTCATAGCAGAGCCCAAACGCCCGCGTCGACACACCCATTca                    >  2:233293/1‑142 (MQ=255)                                                                                                                              tctTCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGGTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACCCGACGt        >  1:834906/1‑143 (MQ=255)                                                                                                                                ttCCGGGCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCCAAACAGGCTTTCGCGCACTCGACGTTg      >  2:818550/1‑143 (MQ=255)                                                                                                                                   agggCTGTATTGCTCTGGCAAGGCAGCCAGTGGGCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGaac   <  1:161505/142‑2 (MQ=255)                                                                                                                                    gggCTGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTTGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAat  <  2:414013/143‑1 (MQ=255)                                                                                                                                        |                                                                                                                                           CCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAAT  >  NC_000913/674633‑674906

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