Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,359,867 G→A 38.0% noncoding (1228/2905 nt) rrlA ← 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,359,8670GA38.0% 73.1 / 78.8 71noncoding (1228/2905 nt)rrlA23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (19/8);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCC  >  Exported/2359815‑2359933
                                                    |                                                                  
cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                  <  1:2614026/71‑1 (MQ=31)
    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                              >  1:982960/1‑71 (MQ=34)
    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                              >  1:1724143/1‑71 (MQ=34)
    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                              <  1:1572449/71‑1 (MQ=34)
    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                              <  1:3142201/71‑1 (MQ=255)
    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                              <  1:346414/71‑1 (MQ=255)
     tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                             >  1:2479960/1‑71 (MQ=255)
      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                            >  1:33901/1‑71 (MQ=34)
          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                          >  1:1429680/1‑69 (MQ=255)
          cGTTACTTATGTCAGCATTCGCACTTCTGATACCTACAGCATGCCTCACAGCACACCTTCACAGGCTTa                                          >  1:1493397/1‑69 (MQ=25)
           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                       <  1:1044328/71‑1 (MQ=34)
           gTTACTTATGTCAGCATTAGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                          >  1:712985/1‑68 (MQ=11)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAAc                                    >  1:3155421/1‑71 (MQ=34)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     <  1:3279495/70‑1 (MQ=255)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     <  1:1831006/70‑1 (MQ=255)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     <  1:3275207/70‑1 (MQ=255)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                    >  1:682856/1‑71 (MQ=255)
                    gTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                              >  1:3267138/1‑71 (MQ=255)
                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:945720/70‑1 (MQ=35)
                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:949668/70‑1 (MQ=35)
                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                           >  1:2036231/1‑71 (MQ=34)
                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                           >  1:1046479/1‑71 (MQ=34)
                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                >  1:782118/1‑66 (MQ=11)
                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                >  1:2731854/1‑66 (MQ=11)
                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                      >  1:57780/1‑55 (MQ=11)
                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:609476/1‑71 (MQ=255)
                            tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:1340370/1‑71 (MQ=255)
                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:260759/1‑70 (MQ=35)
                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     >  1:1683142/1‑55 (MQ=11)
                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     >  1:1748464/1‑55 (MQ=11)
                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            >  1:1110881/1‑64 (MQ=255)
                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            >  1:284289/1‑64 (MQ=255)
                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            >  1:3079029/1‑64 (MQ=255)
                             agCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     >  1:961537/2‑55 (MQ=11)
                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     >  1:1238462/1‑70 (MQ=35)
                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     >  1:477067/1‑70 (MQ=35)
                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:22728/70‑1 (MQ=35)
                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:2558054/1‑67 (MQ=35)
                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                 >  1:2523233/1‑70 (MQ=255)
                                 cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                  >  1:912591/1‑70 (MQ=255)
                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                 >  1:2946382/1‑68 (MQ=37)
                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca               >  1:1851500/1‑70 (MQ=40)
                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca               >  1:1189157/1‑70 (MQ=40)
                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAAAGAACGCTCCCCTACCCAACAacac                >  1:1829806/1‑69 (MQ=255)
                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:2547174/54‑1 (MQ=35)
                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:1948075/54‑1 (MQ=35)
                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:2386666/54‑1 (MQ=35)
                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:1734987/54‑1 (MQ=35)
                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:1430701/54‑1 (MQ=35)
                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACATCACATAGTGTCGc      >  1:2137463/1‑70 (MQ=255)
                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt     >  1:2090222/1‑71 (MQ=255)
                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt     >  1:924656/1‑71 (MQ=255)
                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc        >  1:3188480/1‑67 (MQ=255)
                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc        >  1:676131/1‑67 (MQ=255)
                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg    <  1:2969483/71‑1 (MQ=255)
                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg    <  1:2607136/71‑1 (MQ=255)
                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                <  1:1440657/41‑1 (MQ=38)
                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc  >  1:2328205/1‑71 (MQ=255)
                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc  >  1:968100/1‑71 (MQ=255)
                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc   >  1:1392526/1‑70 (MQ=255)
                                                gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc   >  1:966382/1‑70 (MQ=255)
                                                gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                      <  1:1001653/35‑1 (MQ=14)
                                                gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                      <  1:2688902/35‑1 (MQ=14)
                                                gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                      <  1:1017131/35‑1 (MQ=14)
                                                 cATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                             >  1:2564914/1‑43 (MQ=14)
                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:3023979/49‑1 (MQ=38)
                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:612461/49‑1 (MQ=38)
                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:55083/49‑1 (MQ=38)
                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:3447536/49‑1 (MQ=38)
                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:3200085/49‑1 (MQ=38)
                                                   tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:2543887/49‑1 (MQ=38)
                                                    |                                                                  
CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCC  >  Exported/2359815‑2359933

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