Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,763,636 A→G 14.6% noncoding (208/1542 nt) rrsG ← 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,763,6360AG14.6% 95.0 / 9.6 41noncoding (208/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (11/24);  new base G (2/4);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCC  >  Exported/1763569‑1763699
                                                                   |                                                               
cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCTAGGGt                                                              >  1:1028553/1‑71 (MQ=255)
cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGt                                                              >  1:4804307/1‑71 (MQ=255)
 ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                             >  1:1251940/1‑71 (MQ=255)
 ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                             >  1:4484325/1‑71 (MQ=255)
       aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                        <  1:304188/70‑1 (MQ=255)
       aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                       <  1:2904722/71‑1 (MQ=39)
       aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                       <  1:2615005/71‑1 (MQ=39)
       aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                       <  1:2200549/71‑1 (MQ=39)
            ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                  <  1:1494854/71‑1 (MQ=255)
               ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt                                               <  1:1016439/71‑1 (MQ=255)
               ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                               <  1:2987580/71‑1 (MQ=39)
               ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                               <  1:4917521/71‑1 (MQ=39)
                  caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTgg                                                  <  1:478146/65‑1 (MQ=255)
                  caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                            <  1:4913370/71‑1 (MQ=39)
                    ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                 >  1:3341024/1‑48 (MQ=255)
                     caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                          <  1:3641603/70‑1 (MQ=39)
                     caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:3287513/61‑1 (MQ=255)
                     caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:2732525/61‑1 (MQ=255)
                     caacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:4468679/61‑1 (MQ=255)
                      aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                          <  1:3915921/69‑1 (MQ=37)
                      aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAAGCCCGAAGGTCCCCCTCTTTg                                                   <  1:3690255/60‑1 (MQ=255)
                       acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                              >  1:419489/1‑48 (MQ=38)
                          aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:4272631/71‑1 (MQ=34)
                          aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:3532497/71‑1 (MQ=34)
                          aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:739830/71‑1 (MQ=34)
                          aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:1495861/71‑1 (MQ=34)
                          aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTAt                                    <  1:1768344/71‑1 (MQ=34)
                           gCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                   >  1:3521586/1‑71 (MQ=34)
                           gCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                   <  1:2947079/71‑1 (MQ=34)
                           gCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATg                                   >  1:1528437/1‑71 (MQ=34)
                                 cccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                        >  1:488861/1‑60 (MQ=31)
                                 cccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                        >  1:4744069/1‑60 (MQ=31)
                                         ggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                         <  1:2985546/67‑1 (MQ=35)
                                                  gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt            <  1:3109179/71‑1 (MQ=34)
                                                        aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt      <  1:4078482/71‑1 (MQ=34)
                                                           agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:3509991/1‑71 (MQ=34)
                                                           agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:1060731/1‑71 (MQ=34)
                                                           agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:1997506/1‑71 (MQ=34)
                                                           agGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc   >  1:3879270/1‑71 (MQ=34)
                                                            ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  <  1:4318533/71‑1 (MQ=34)
                                                            ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  <  1:463664/71‑1 (MQ=34)
                                                            ggCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc  <  1:1430355/71‑1 (MQ=34)
                                                                   aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:2229688/1‑36 (MQ=40)
                                                                   aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:1392609/1‑36 (MQ=40)
                                                                   aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:1178561/1‑36 (MQ=40)
                                                                   aGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                              >  1:3202234/1‑36 (MQ=40)
                                                                   aGGTCCCCCTCTTTGGTCTTGCGACATTATGCAGt                               >  1:161583/1‑35 (MQ=255)
                                                                   |                                                               
CTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCC  >  Exported/1763569‑1763699

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: