Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,018,148 C→A 9.6% R31S (CGT→AGT)  mshA_4 → D‑inositol‑3‑phosphate glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,018,1480CA 9.6% 127.5 / 4.5 49R31S (CGT→AGT) mshA_4D‑inositol‑3‑phosphate glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base A (3/2);  total (27/22)
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.14e-01

TGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGACTTGCATGTCGCTTCAACGGGTTTTTCAC  >  KT2440/2018014‑2018277
                                                                                                                                      |                                                                                                                                 
tGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCAccc                                                                                                                                   >  3:6133/1‑135 (MQ=255)
          gTACAGTCGAGGGATACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACgg                                                                                                                         <  2:114306/135‑1 (MQ=255)
                  gAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGAATTTTTTTAAAACGTATTTCCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  5:281667/135‑1 (MQ=255)
                  gAGGTTTACAAGGCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  3:158636/135‑1 (MQ=255)
                   aGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  5:200299/134‑1 (MQ=255)
                   aGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  6:233907/134‑1 (MQ=255)
                   aGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  3:53816/134‑1 (MQ=255)
                   aGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  7:13676/134‑1 (MQ=255)
                   aGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTAACCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  8:230189/134‑1 (MQ=255)
                     gTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCag                                                                                                              <  4:84915/135‑1 (MQ=255)
                      tttACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCCTTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATAGCCGAAGGCACCCGTGAGCACGGCGTGCCATCaga                                                                                                             >  6:188523/1‑135 (MQ=255)
                      tttACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 >  7:215390/1‑131 (MQ=255)
                      tttACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAAt                                                                                                                 <  8:215390/131‑1 (MQ=255)
                      tttACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCaga                                                                                                             >  2:205779/1‑135 (MQ=255)
                      tttACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCag                                                                                                              <  4:174273/134‑1 (MQ=255)
                      tttACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCAGTGAGCACGGCGTGCAATCaga                                                                                                             >  6:262948/1‑135 (MQ=255)
                      tttACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCAGTGAGCACGGCGTGCAATCaga                                                                                                             >  6:258493/1‑135 (MQ=255)
                      tttACAAGTCGGTCAACGAACGATGATAAAGGTTCTGCATTTTTTTAAAACGTAATACCCGGACTCCATGGGTGGCATCGCTCAGCTCATTTTTCAGATCGCCGAAGGCCCCCGTGAGCACGGCGTGCCATCaca                                                                                                             >  2:285161/1‑133 (MQ=255)
                      tatACAAGTCGGTCAAGGAACGATGATCAAGGATCAGCAGTTTTTTAAAACGTATTCACCGGACTCCATGGGTGGCATCCAACAGTGCATTTTTCAGATAGAAGAAGGCACCAGAGATAACGCCGTGCCATCaga                                                                                                             >  4:101386/3‑135 (MQ=255)
                        tACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAgg                                                                                                           >  2:4590/1‑135 (MQ=255)
                          cAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGt                                                                                                          <  2:131008/134‑1 (MQ=255)
                          cAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGt                                                                                                          <  1:226907/134‑1 (MQ=255)
                          cAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCAGTGAGCACGGCGTGCAATCAGAGGt                                                                                                          <  5:258493/134‑1 (MQ=255)
                          cAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCAGTGAGCACGGCGTGCAATCAGAGGt                                                                                                          <  5:262948/134‑1 (MQ=255)
                           aaGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTg                                                                                                         >  1:142335/1‑134 (MQ=255)
                                   cAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATAcc                                                                                                  >  3:146110/1‑133 (MQ=255)
                                   cAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTc                                                                                                >  6:21196/1‑135 (MQ=255)
                                     aGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAg                                                                                              >  6:98871/1‑135 (MQ=255)
                                        aaCGAAGATCCAGGATCTGCCTTTTTTTAAACCGTATTACCCGGACTCCCCGGGTGCCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGGGCAATCAGAGGTGCTATACCTCAGccc                                                                                           <  6:58772/135‑1 (MQ=255)
                                               atCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTggg                                                                                    >  6:176547/1‑135 (MQ=255)
                                               atCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTggg                                                                                    >  2:277450/1‑135 (MQ=255)
                                                  aaGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCt                                                                                 >  1:108303/1‑135 (MQ=255)
                                                   aGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCt                                                                                 <  2:108303/134‑1 (MQ=255)
                                                                 tttAAATCGTATTAGCCGGCCTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAAcc                                                                   <  1:113696/134‑1 (MQ=255)
                                                                   tAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCg                                                                  >  1:62442/1‑133 (MQ=255)
                                                                          tATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGc                                                         <  2:62442/135‑1 (MQ=255)
                                                                                  ggACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGAc                                                 <  5:98871/135‑1 (MQ=255)
                                                                                       ccATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCAcc                                            >  3:282828/1‑135 (MQ=255)
                                                                                       ccATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCAcc                                            >  2:100043/1‑135 (MQ=255)
                                                                                             gTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCCCCGCTCCa                                      >  4:189296/1‑135 (MQ=255)
                                                                                                 cATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCACGCt                                  >  8:63727/1‑135 (MQ=255)
                                                                                                 cATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCt                                  >  3:229053/1‑135 (MQ=255)
                                                                                                   tCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTgg                                >  8:137306/1‑135 (MQ=255)
                                                                                                     gAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGAc                              >  4:141489/1‑135 (MQ=255)
                                                                                                            tCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGACTTGCATg                       <  2:225004/135‑1 (MQ=255)
                                                                                                              aTTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGACTTGCATGTc                     >  4:184104/1‑135 (MQ=255)
                                                                                                              aTTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGACTTGCATGTc                     >  3:148612/1‑135 (MQ=255)
                                                                                                                        tCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGACTTGCATGTCGCTTCAACgg           >  7:126242/1‑135 (MQ=255)
                                                                                                                        tCGCCGAAGGAACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCTAGCTGGACTTGCTTGTCGCTTCAACgg           >  6:23330/1‑135 (MQ=255)
                                                                                                                                 gCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGACTTGCATGTCGCTTCAACGGGTTTTTCAc  <  2:142335/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                 
TGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAAGCTGGACTTGCATGTCGCTTCAACGGGTTTTTCAC  >  KT2440/2018014‑2018277

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