Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,505,955 A→G 8.3% intergenic (+271/‑481) rsmI → / → mraZ Ribosomal RNA small subunit methyltransferase I/Transcriptional regulator MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,505,9550AG 8.3% 139.5 / 4.2 48intergenic (+271/‑481)rsmI/mraZRibosomal RNA small subunit methyltransferase I/Transcriptional regulator MraZ
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base G (2/2);  total (25/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 = 4.76e-01

AAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCCGTCGTAGAGGAATG  >  KT2440/1505821‑1506086
                                                                                                                                             |                                                                                                                                   
aaaGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGcgcaccgca                                                                                                                                     >  3:227101/1‑135 (MQ=255)
 aaGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGcgcaccgcac                                                                                                                                    >  8:7128/1‑135 (MQ=255)
  aGTCCGGGCTCCTCAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACg                                                                                                                                   <  5:51166/135‑1 (MQ=255)
  aGTCCGGGCTCCTCAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACg                                                                                                                                   <  4:109865/135‑1 (MQ=255)
      cGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTg                                                                                                                               >  4:86190/1‑135 (MQ=255)
                  ggCGGAGGGCCAGGTCACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAACAGGTTCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGt                                                                                                                   <  2:144013/135‑1 (MQ=255)
                       aGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTa                                                                                                              <  7:192881/135‑1 (MQ=255)
                          gCCAGGTAACTCCTGGGAGGTGTGAGCCTACGGAAAGTGCCACAGAAAATATACCGCCCACTTTTCCTTTGGGAATCGTAGGTAAGGGTGAAATGGTGCGGTAAGAGCGCACCGCGCGATTGGTAACAG‑TCGTgg                                                                                                                 >  5:263066/1‑135 (MQ=255)
                          gCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAgg                                                                                                            >  6:250558/1‑134 (MQ=255)
                          gCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGATCGTGGCTAgg                                                                                                            >  3:200535/1‑134 (MQ=255)
                               gTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAAcc                                                                                                      >  2:287412/1‑135 (MQ=255)
                                      cggggTGGCGTGAGCCTACGTAAAGTGCCACTGAAAATA‑ACCGC‑‑‑GACAGCAGTTCGG‑‑G‑GCCGGTAAGGGTGAAAAGGGGCGGTAAGAGCGCACCGCGCGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCg                                                                                               <  1:128315/133‑1 (MQ=255)
                                        gggAGGCGTGAGTCTACGGAGAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGa                                                                                             >  1:208368/1‑135 (MQ=255)
                                        gggAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGa                                                                                             <  1:232906/135‑1 (MQ=255)
                                                agAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGAGGCTAGGTAAACCCCACTCGGAGCAAGAc                                                                                      >  4:36177/2‑134 (MQ=255)
                                                 gAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCt                                                                                                                                >  2:257873/1‑91 (MQ=255)
                                                 gAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCt                                                                                                                                <  1:257873/91‑1 (MQ=255)
                                                 gAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCg                                                                                               <  2:226799/124‑1 (MQ=255)
                                                 gAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCg                                                                                               >  1:226799/1‑124 (MQ=255)
                                                 gAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGAcc                                                                                     >  5:98634/1‑134 (MQ=255)
                                                 gAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCa                                                                                    <  4:235996/135‑1 (MQ=255)
                                                           aaaGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCCCTCGGAGCAAGACCAAATAggcta                                                                           >  6:231512/1‑131 (MQ=255)
                                                           aaaGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTc                                                                          <  5:78922/135‑1 (MQ=255)
                                                                gCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCa                                                                     >  6:146153/1‑135 (MQ=255)
                                                                 ccACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAg                                                                    >  2:255195/1‑135 (MQ=255)
                                                                    caGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAACAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCg                                                                 >  5:63411/1‑135 (MQ=255)
                                                                      gAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTg                                                               <  1:207014/135‑1 (MQ=255)
                                                                        aaaTA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGc                                                             >  3:62692/1‑135 (MQ=255)
                                                                                        aaGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAAc                                                 >  6:174481/1‑135 (MQ=255)
                                                                                          gCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAAcc                                                <  7:7128/134‑1 (MQ=255)
                                                                                          gCACTTCGG‑‑G‑GCAGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCg                                               <  2:143748/135‑1 (MQ=255)
                                                                                            aCTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCggg                                             <  4:23318/135‑1 (MQ=255)
                                                                                              ttCGG‑‑T‑GCCGGTAAGGGTGAAAAGGGGCGGTAACAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTa                                           <  6:63411/135‑1 (MQ=255)
                                                                                                  g‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAAcccc                                                                                                    >  4:30608/1‑74 (MQ=255)
                                                                                                  g‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAAcccc                                                                                                    <  3:30608/74‑1 (MQ=255)
                                                                                                     t‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTg                                      >  1:113822/1‑135 (MQ=255)
                                                                                                     t‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTg                                      >  5:258784/1‑135 (MQ=255)
                                                                                                         cGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAACTAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTa                                   <  8:229926/135‑1 (MQ=255)
                                                                                                         cGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTa                                   <  2:38282/135‑1 (MQ=255)
                                                                                                                  tGAAAGGGTGCGGTAAGAGCGCACCGCGCGGCTGGTAACAGTCCGTGGCACGGTAAACTCCACCCGGAGCa                                                                                          <  8:249582/71‑1 (MQ=255)
                                                                                                                  tGAAAGGGTGCGGTAAGAGCGCACCGCGCGGCTGGTAACAGTCCGTGGCACGGTAAACTCCACCCGGAGCa                                                                                          >  7:249582/1‑71 (MQ=255)
                                                                                                                   gAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCa                         <  2:208368/135‑1 (MQ=255)
                                                                                                                      aaGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTg                      >  3:116263/1‑135 (MQ=255)
                                                                                                                      aaGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTg                      >  7:215557/1‑135 (MQ=255)
                                                                                                                       aGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGa                     >  2:19369/1‑135 (MQ=255)
                                                                                                                        ggTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGAt                    <  6:98634/135‑1 (MQ=255)
                                                                                                                        cgTGCGGTAAGAGCGCACCGCAAGGCTGGCTACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGTTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGAt                    <  5:231512/134‑1 (MQ=255)
                                                                                                                             ggTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGACCCGGGTAGGTTGCTAGCGGCGTCCTGTGATGGCcg               >  1:198559/1‑135 (MQ=255)
                                                                                                                              gTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCcgt              >  1:103437/1‑135 (MQ=255)
                                                                                                                                aaGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCcgtcg            <  4:62692/135‑1 (MQ=255)
                                                                                                                                      gcaccgcacGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCCGTCGTagag       >  4:232015/1‑134 (MQ=255)
                                                                                                                                        accgcacGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCCGTCGTAGAGGaa    <  1:287412/135‑1 (MQ=255)
                                                                                                                                          cgcacGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCCGTCGTAGAGGAATg  >  1:188996/1‑135 (MQ=255)
                                                                                                                                          cgcacGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCCGTCGTAGAGGAATg  >  4:34691/1‑135 (MQ=255)
                                                                                                                                             |                                                                                                                                   
AAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATA‑ACCGC‑‑‑CTAAGCACTTCGG‑‑T‑GCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGGCCCGCGTCGGAACCGGGTAGGTTGCTAGAGGCGTCCAGTGATGGCCGTCGTAGAGGAATG  >  KT2440/1505821‑1506086

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: ‑