Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,763,558 G→C 7.6% noncoding (286/1542 nt) rrsG ← 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,763,5580GC 7.6% 242.4 / 8.6 79noncoding (286/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base C (3/3);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAA  >  Exported/1763502‑1763626
                                                        |                                                                    
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                        <  1:153424/71‑1 (MQ=11)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                        >  1:1287326/1‑71 (MQ=32)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                         <  1:2053499/70‑1 (MQ=14)
   gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                      >  1:245794/1‑70 (MQ=35)
      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTgg                                                       >  1:2854884/1‑66 (MQ=35)
      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                  >  1:2822562/1‑71 (MQ=34)
      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                  >  1:1542170/1‑71 (MQ=21)
        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 <  1:36093/70‑1 (MQ=35)
        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 <  1:1533134/70‑1 (MQ=35)
        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                   <  1:54745/68‑1 (MQ=25)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               <  1:92239/71‑1 (MQ=34)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               <  1:2291173/71‑1 (MQ=34)
          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                              <  1:322810/71‑1 (MQ=34)
          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               >  1:725313/1‑70 (MQ=35)
          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               >  1:34137/1‑70 (MQ=35)
           cGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                             >  1:2234487/1‑71 (MQ=11)
           cGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                             >  1:2064553/1‑71 (MQ=11)
            gtgtCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAc                                            >  1:2812637/1‑71 (MQ=11)
            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGAACAGCCAGGGATCGTCGCCtt                                                         >  1:612039/1‑58 (MQ=11)
              gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCATACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                            <  1:629313/69‑1 (MQ=255)
              gtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                          >  1:117713/1‑71 (MQ=34)
               tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcctc                                         <  1:1445443/71‑3 (MQ=11)
               tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcctc                                         <  1:2114044/71‑3 (MQ=11)
               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                               <  1:1206147/65‑1 (MQ=35)
               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                          <  1:2317573/70‑1 (MQ=35)
                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:1441911/1‑61 (MQ=35)
                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:2419411/1‑61 (MQ=35)
                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                       >  1:1439129/1‑71 (MQ=21)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                     <  1:2998771/71‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                     <  1:805551/71‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                      >  1:2588559/1‑70 (MQ=35)
                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   <  1:1069538/71‑1 (MQ=39)
                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   >  1:553835/1‑70 (MQ=255)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:567127/1‑55 (MQ=35)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 >  1:1399471/1‑55 (MQ=35)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                  <  1:1426133/70‑1 (MQ=255)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                  <  1:2970047/70‑1 (MQ=255)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   >  1:708595/1‑69 (MQ=40)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                   >  1:2829712/1‑69 (MQ=40)
                       ccAGTGTGGCGGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGcc                                                 <  1:2338003/55‑1 (MQ=255)
                        cAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                <  1:2579144/71‑1 (MQ=255)
                        cAGTGTGGCTGGTCATCCTCCCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                                 <  1:1201307/70‑1 (MQ=255)
                         aGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                                <  1:1998922/70‑1 (MQ=255)
                         aGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                                 <  1:2043415/69‑1 (MQ=255)
                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                            >  1:2376979/1‑71 (MQ=255)
                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGAGAGCCGTTACCCCACCAACAAGCTaa                            >  1:857008/1‑71 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                            <  1:972724/70‑1 (MQ=255)
                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                         >  1:1797169/1‑71 (MQ=255)
                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                      <  1:996444/41‑1 (MQ=38)
                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                      <  1:2720448/41‑1 (MQ=38)
                                cTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                         >  1:2771500/1‑70 (MQ=255)
                                  ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                            <  1:1307396/65‑1 (MQ=255)
                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                         <  1:1603974/35‑1 (MQ=33)
                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                         <  1:517172/35‑1 (MQ=33)
                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                         <  1:2014764/35‑1 (MQ=33)
                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                     <  1:112475/71‑1 (MQ=255)
                                   gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                      <  1:1798062/70‑1 (MQ=255)
                                   gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                      <  1:766141/70‑1 (MQ=255)
                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                  >  1:236996/1‑71 (MQ=255)
                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca               >  1:2567421/1‑71 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                             >  1:1407705/1‑41 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                             >  1:2310381/1‑41 (MQ=255)
                                          tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                 >  1:148330/1‑68 (MQ=255)
                                          tctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca               >  1:1789183/1‑70 (MQ=255)
                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                             >  1:1160361/1‑55 (MQ=255)
                                           ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                             >  1:1535652/1‑55 (MQ=255)
                                            tcAGACCAGCTAGGGATCGACGCCTTGGTGAGCCGTTa                                             >  1:2522269/1‑38 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          <  1:677553/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          <  1:1231157/71‑1 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt            <  1:828174/69‑1 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          >  1:2513018/1‑71 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc          <  1:1356616/71‑1 (MQ=255)
                                               gACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg         <  1:521077/71‑1 (MQ=255)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa        <  1:915099/71‑1 (MQ=255)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa        <  1:928956/71‑1 (MQ=255)
                                                  cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                       >  1:1707202/1‑54 (MQ=255)
                                                  cAGCTAGGGATCGTCCCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                       >  1:303607/1‑54 (MQ=255)
                                                    gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc    <  1:388795/71‑1 (MQ=255)
                                                      caGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa  <  1:1840125/70‑1 (MQ=255)
                                                        |                                                                    
GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAA  >  Exported/1763502‑1763626

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