Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,884,163 G→A 12.2% A107V (GCA→GTA)  rumA ← 23S rRNA (uracil‑5)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,884,1630GA12.2% 221.6 / 17.1 82A107V (GCA→GTA) rumA23S rRNA (uracil‑5)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (48/24);  new base A (5/5);  total (53/29)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCATACGCCAA  >  Exported/1884097‑1884232
                                                                  |                                                                     
cGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGc                                                                   <  1:2564668/71‑1 (MQ=255)
 ggCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCa                                                                  >  1:3518171/1‑71 (MQ=255)
   cGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACt                                                                <  1:5000422/71‑1 (MQ=255)
      tAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACttt                                                              <  1:2930917/70‑1 (MQ=255)
      tAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACttt                                                              <  1:3241174/70‑1 (MQ=255)
      tAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACttt                                                              <  1:671632/70‑1 (MQ=255)
        gCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCCCTTTTg                                                            <  1:668480/70‑1 (MQ=255)
        gCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGc                                                           <  1:2400691/71‑1 (MQ=255)
        gCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACtttt                                                             >  1:1819413/1‑69 (MQ=255)
        gCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACTTTTg                                                            >  1:1912843/1‑70 (MQ=255)
        gCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACTTTTg                                                            >  1:230310/1‑70 (MQ=255)
        gCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACTTTTg                                                            >  1:2393985/1‑70 (MQ=255)
        gCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACTTTTg                                                            >  1:4306419/1‑70 (MQ=255)
          cccAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTACCGCACttt                                                              <  1:2688621/66‑1 (MQ=255)
                 aaCATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCgctgct                                                  >  1:2971957/1‑71 (MQ=255)
                 aaCATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCgctgct                                                  >  1:1488176/1‑71 (MQ=255)
                 aaCATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCgctgct                                                  >  1:992253/1‑71 (MQ=255)
                 aaCATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCgctgct                                                  >  1:4669468/1‑71 (MQ=255)
                 aaCATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCgctgc                                                   >  1:1572524/1‑70 (MQ=255)
                   cATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGt                                                <  1:2414434/71‑1 (MQ=255)
                   cATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGt                                                <  1:1450906/71‑1 (MQ=255)
                     tCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         <  1:1090202/60‑1 (MQ=255)
                     tCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCg                                                       >  1:3821525/1‑62 (MQ=255)
                     tCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCg                                                       >  1:3161042/1‑62 (MQ=255)
                     tCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCg                                                       >  1:3154837/1‑62 (MQ=255)
                     tCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCg                                                       >  1:780883/1‑62 (MQ=255)
                      cGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTaaa                                             >  1:875944/1‑71 (MQ=255)
                        gCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         <  1:1418646/57‑1 (MQ=255)
                           aTCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTc                                                        <  1:3774147/55‑1 (MQ=255)
                           aTCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCAc                                        >  1:186852/1‑71 (MQ=255)
                           aTCACTTCAGAGACATCGTGATTCATTAATCGGGCGAGTACCGCACttt                                                              <  1:3003004/49‑1 (MQ=38)
                              aCTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCt                                     >  1:512646/1‑71 (MQ=255)
                                  cAGAGACTTCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCa                                         >  1:1323151/1‑63 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACtt                                                               >  1:1164913/1‑41 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         >  1:2890403/1‑47 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         >  1:4117549/1‑47 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         >  1:4027824/1‑47 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         >  1:3633978/1‑47 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         >  1:3599245/1‑47 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         >  1:3158493/1‑47 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCtt                                                         >  1:2825661/1‑47 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCa                                         >  1:3143274/1‑63 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCa                                         >  1:2558516/1‑63 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCa                                         >  1:3435782/1‑63 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCa                                         >  1:3468699/1‑63 (MQ=255)
                                  cAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCa                                         >  1:2932444/1‑63 (MQ=255)
                                        cATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGtgttgt                           <  1:2804672/71‑1 (MQ=255)
                                         aTCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTgttgtt                          >  1:4525975/1‑71 (MQ=255)
                                         aTCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTgttgtt                          >  1:2664804/1‑71 (MQ=255)
                                          tCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATc                                           <  1:1442597/53‑1 (MQ=255)
                                          tCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATc                                           <  1:1669993/53‑1 (MQ=255)
                                          tCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATc                                           <  1:2586251/53‑1 (MQ=255)
                                               tttCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCt                       <  1:4668443/68‑1 (MQ=255)
                                               tttCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGCAAATCCACGCTGGCGTGttgttg                         <  1:597395/66‑1 (MQ=255)
                                               tttCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGCAAATCCACGCTGGCGTGttgttg                         <  1:555096/66‑1 (MQ=255)
                                               tttCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGCAAATCCACGCTGGCGTGttgttg                         <  1:4353136/66‑1 (MQ=255)
                                                ttCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGa                     >  1:3020820/1‑69 (MQ=255)
                                                ttCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGa                     >  1:4684658/1‑69 (MQ=255)
                                                ttCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGa                     >  1:2740244/1‑69 (MQ=255)
                                                 tCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAg                  >  1:1660240/1‑71 (MQ=255)
                                                 tCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAg                  >  1:2913785/1‑71 (MQ=255)
                                                 tCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAg                  >  1:2213214/1‑71 (MQ=255)
                                                 tCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAg                  >  1:2679409/1‑71 (MQ=255)
                                                 tCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAg                  >  1:5080731/1‑71 (MQ=255)
                                                 tCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAg                  >  1:2434288/1‑71 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAAt                                            >  1:2799467/1‑40 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAAt                                            >  1:3411319/1‑40 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAAt                                            >  1:171219/1‑40 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAAt                                            >  1:2965905/1‑40 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAAt                                            >  1:5142365/1‑40 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCg                                 >  1:3608888/1‑51 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCg                                 >  1:2647762/1‑51 (MQ=255)
                                                      aaTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCg                                 >  1:4103195/1‑51 (MQ=255)
                                                       aTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCAcct            <  1:4642529/71‑2 (MQ=255)
                                                       aTCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCg            <  1:4636158/71‑1 (MQ=255)
                                                       aTCGGGCGAGTGCCGCACTTTTGATTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCg            <  1:2416813/71‑1 (MQ=255)
                                                          gggCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCa          <  1:524512/70‑1 (MQ=255)
                                                          gggCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCa          <  1:2862973/70‑1 (MQ=255)
                                                          gggCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCa          <  1:452873/70‑1 (MQ=255)
                                                                gTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCATACGcc    >  1:4343948/1‑70 (MQ=255)
                                                                gTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCATACGcc    >  1:4209712/1‑70 (MQ=255)
                                                                 tGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCATACGCCaa  <  1:3415817/71‑1 (MQ=255)
                                                                  |                                                                     
CGGCGATAGCCCCAGGGAACATCGGCGATCACTTCAGAGACATCGTGTTTCATTAATCGGGCGAGTGCCGCACTTTTGCTTCGCTGCTGTAAATCCACGCTGGCGTGTTGTTGCTGACAGCCACCGCATACGCCAA  >  Exported/1884097‑1884232

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