Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,763,556 T→C 14.0% noncoding (288/1542 nt) rrsG ← 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,763,5560TC14.0% 212.6 / 24.3 93noncoding (288/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (35/45);  new base C (8/5);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  Exported/1763502‑1763624
                                                      |                                                                    
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      >  1:5382931/1‑71 (MQ=32)
 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:3474438/70‑1 (MQ=33)
 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     >  1:147033/1‑71 (MQ=21)
 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTgg                                                     >  1:338465/1‑71 (MQ=21)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:4613090/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    <  1:2342091/71‑1 (MQ=34)
    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:1920382/66‑1 (MQ=31)
    tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:5486318/66‑1 (MQ=31)
      tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                >  1:1243845/1‑71 (MQ=34)
      tGGACCGTGTCTCAGTTCCAGTGCGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       >  1:1500206/1‑64 (MQ=14)
        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                 <  1:2845927/68‑1 (MQ=35)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:4500763/71‑1 (MQ=34)
          ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:1419050/71‑1 (MQ=34)
            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                 <  1:3250513/64‑1 (MQ=25)
             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:1540128/1‑71 (MQ=34)
             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:4931068/1‑71 (MQ=34)
             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:3250443/1‑71 (MQ=21)
             tgtCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTAcc                                         >  1:2313179/1‑71 (MQ=11)
               tctcATTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                                              <  1:5240251/64‑1 (MQ=25)
               tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTAcct                                        <  1:1611581/70‑2 (MQ=11)
               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:2737638/56‑1 (MQ=33)
               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                                              <  1:5208726/64‑1 (MQ=35)
               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                                              <  1:4462332/64‑1 (MQ=35)
               tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                        <  1:2445349/70‑1 (MQ=25)
                ctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGcc                                               <  1:1115589/62‑1 (MQ=25)
                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            >  1:1606252/1‑64 (MQ=35)
                 tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            >  1:5653006/1‑64 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         <  1:1221911/65‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:1119337/71‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:3574374/71‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacca                                   <  1:5240096/71‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     <  1:5592627/69‑1 (MQ=35)
                    gTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                  <  1:290378/71‑1 (MQ=38)
                    gTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                  <  1:3154924/71‑1 (MQ=38)
                     ttCCAGTGTGGCTGGTCATCCTCTCAGCCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcccc                                       >  1:1766172/1‑65 (MQ=255)
                     ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaa                                  <  1:5250760/70‑1 (MQ=39)
                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:2096116/1‑70 (MQ=40)
                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:693361/1‑70 (MQ=40)
                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:1804816/1‑70 (MQ=40)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:5462507/48‑1 (MQ=35)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 >  1:2582220/1‑69 (MQ=40)
                        cAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                               <  1:4756359/70‑1 (MQ=255)
                         aGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGc                             <  1:2889981/71‑1 (MQ=11)
                         aGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                             <  1:3019785/71‑1 (MQ=255)
                          gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:3962393/71‑1 (MQ=255)
                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:4508462/1‑56 (MQ=35)
                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          >  1:3124287/1‑71 (MQ=255)
                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          >  1:4387212/1‑71 (MQ=255)
                            gtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           >  1:1710521/1‑70 (MQ=255)
                             tgGCTGGTCATCCTCTCAGCCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                         >  1:4295387/1‑71 (MQ=255)
                               gCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAATca                       <  1:127246/71‑2 (MQ=11)
                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:1531678/66‑1 (MQ=255)
                               gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                        >  1:2422744/1‑70 (MQ=255)
                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                          >  1:306370/1‑66 (MQ=255)
                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:4770613/1‑71 (MQ=255)
                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:1197977/1‑71 (MQ=255)
                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:5263956/1‑71 (MQ=255)
                                 tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:3369130/1‑71 (MQ=255)
                                  ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:4941676/1‑46 (MQ=38)
                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:529456/69‑1 (MQ=255)
                                   gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   <  1:4141065/71‑1 (MQ=255)
                                    tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   <  1:3600084/70‑1 (MQ=255)
                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                               >  1:2067189/1‑57 (MQ=255)
                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                 <  1:2528231/70‑1 (MQ=255)
                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:1280253/44‑1 (MQ=38)
                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:988604/44‑1 (MQ=38)
                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:1946398/71‑1 (MQ=255)
                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:4996728/1‑71 (MQ=255)
                                       tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                >  1:1616559/1‑70 (MQ=255)
                                        ccTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                            <  1:4674745/41‑1 (MQ=38)
                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:206717/1‑70 (MQ=255)
                                        ccTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:2850317/1‑70 (MQ=255)
                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           >  1:2559796/1‑57 (MQ=255)
                                         ctctCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGc              >  1:4820539/1‑70 (MQ=255)
                                           ctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca           >  1:4347826/1‑71 (MQ=255)
                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                       <  1:1780049/58‑1 (MQ=255)
                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                       <  1:974581/58‑1 (MQ=255)
                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                       <  1:267960/58‑1 (MQ=255)
                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac            <  1:1028295/69‑1 (MQ=255)
                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt          >  1:4302521/1‑71 (MQ=255)
                                             cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc         >  1:2704789/1‑71 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:1747788/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:444188/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        >  1:1514139/1‑71 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:143287/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:2203260/71‑1 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:2951737/71‑1 (MQ=255)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                             <  1:5031461/48‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:3106954/49‑1 (MQ=255)
                                                  cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg    >  1:1766011/1‑71 (MQ=255)
                                                  cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg    >  1:2722573/1‑71 (MQ=255)
                                                  cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg    >  1:3862394/1‑71 (MQ=255)
                                                  cAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATg    >  1:5199526/1‑71 (MQ=255)
                                                      tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                               >  1:3672310/1‑40 (MQ=255)
                                                      tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaacaa                               >  1:3912247/1‑40 (MQ=255)
                                                      tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            >  1:1195907/1‑43 (MQ=255)
                                                      tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           >  1:1705776/1‑44 (MQ=255)
                                                      tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                    >  1:1145283/1‑51 (MQ=255)
                                                      tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   >  1:1446288/1‑52 (MQ=255)
                                                      tAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  >  1:5686103/1‑69 (MQ=255)
                                                      |                                                                    
GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  Exported/1763502‑1763624

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: