Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,763,580 G→T 16.5% noncoding (264/1542 nt) rrsG ← 16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,763,5800GT16.5% 250.8 / 23.7 85noncoding (264/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (28/42);  new base T (4/10);  total (32/53)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTG  >  minE/1763510‑1763650
                                                                      |                                                                      
gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                                                                        >  1:1050636/1‑71 (MQ=34)
gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCg                                                                        >  1:2667952/1‑71 (MQ=34)
gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                                                                        >  1:268553/1‑71 (MQ=21)
gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCg                                                                        >  1:2768645/1‑71 (MQ=21)
 aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                                       >  1:2791192/1‑71 (MQ=21)
 aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                                       <  1:2657554/71‑1 (MQ=21)
 aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                                                       <  1:2984759/71‑1 (MQ=21)
  ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                                      <  1:2269088/71‑1 (MQ=34)
  ccGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                                      <  1:1534099/71‑1 (MQ=34)
   cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                                      >  1:2770499/1‑70 (MQ=35)
   cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGtt                                                                      <  1:1459344/70‑1 (MQ=35)
   cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTa                                                                     >  1:1172289/1‑71 (MQ=34)
       tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                                                  >  1:3082721/1‑70 (MQ=25)
         tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                                               >  1:2117497/1‑71 (MQ=34)
           aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGtt                                                                      <  1:2906679/62‑1 (MQ=25)
           aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                                                     >  1:2429005/1‑63 (MQ=25)
               ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                                          <  1:1784568/70‑1 (MQ=255)
               ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                                          <  1:2341536/70‑1 (MQ=255)
               ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                                          <  1:2534359/70‑1 (MQ=255)
               ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                                           >  1:143677/1‑69 (MQ=40)
                   tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                                     <  1:203697/71‑1 (MQ=255)
                   tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                                     <  1:2716738/71‑1 (MQ=255)
                   tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                                     <  1:3391515/71‑1 (MQ=255)
                   tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                                                     <  1:123096/71‑1 (MQ=255)
                   tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACAAGCTa                                                     <  1:1654433/71‑1 (MQ=255)
                     tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                                                    <  1:2333982/70‑1 (MQ=255)
                     tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAt                                                   <  1:2620077/71‑1 (MQ=255)
                       gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                                                    <  1:3071650/52‑1 (MQ=255)
                       gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                                                    <  1:1472666/52‑1 (MQ=255)
                       gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                                                  <  1:2296588/70‑1 (MQ=255)
                       gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                                                  <  1:736120/70‑1 (MQ=255)
                         tGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTAAt                                                   >  1:42192/1‑67 (MQ=11)
                           atCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATccga                                               <  1:1856091/68‑3 (MQ=11)
                            tCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                                            <  1:2828359/71‑1 (MQ=255)
                               tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                                         >  1:3225207/1‑71 (MQ=255)
                               tCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg                                         >  1:1322170/1‑71 (MQ=255)
                               tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTgg                                          <  1:3165672/70‑1 (MQ=255)
                                  tctcAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca                                       >  1:257865/1‑70 (MQ=255)
                                   ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                                                       >  1:2983573/1‑53 (MQ=255)
                                   ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                                                       >  1:2533351/1‑53 (MQ=255)
                                    tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                                                   >  1:1039970/1‑40 (MQ=40)
                                     cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                                            >  1:770004/1‑62 (MQ=255)
                                     cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                                            >  1:3144855/1‑62 (MQ=255)
                                      aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAt                                    <  1:3278459/69‑1 (MQ=255)
                                      aGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTaa                                                    <  1:2266679/53‑1 (MQ=11)
                                      aGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTaa                                                    <  1:2262613/53‑1 (MQ=11)
                                      aGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTaa                                                    <  1:359788/53‑1 (MQ=11)
                                      aGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTACCTCACCAACTAGCTaa                                                    <  1:662917/53‑1 (MQ=11)
                                        aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa                                <  1:2777904/71‑1 (MQ=255)
                                         ccAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt                               >  1:298506/1‑71 (MQ=255)
                                           aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                                             <  1:2272994/39‑1 (MQ=25)
                                           aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                                             <  1:2374791/39‑1 (MQ=25)
                                           aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                                             <  1:3302956/39‑1 (MQ=25)
                                           aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                                             <  1:1545529/39‑1 (MQ=25)
                                           aGCTACTGATCGTCGCCTTGGTGAGCCTTTACCtcacca                                                             <  1:1166168/39‑1 (MQ=25)
                                           aGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg                             <  1:1970495/71‑1 (MQ=255)
                                            gccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc                            <  1:873395/68‑1 (MQ=255)
                                            gccaGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc                            <  1:2416041/68‑1 (MQ=255)
                                             cTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCa                           <  1:2564884/71‑1 (MQ=255)
                                               acggATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCt                                                      >  1:1614449/3‑42 (MQ=14)
                                                gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                                                      <  1:1278351/41‑1 (MQ=255)
                                                gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa                          <  1:1683061/69‑1 (MQ=255)
                                                gggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCaa                          <  1:3391713/69‑1 (MQ=255)
                                                 ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:2654497/60‑1 (MQ=255)
                                                 ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:1444325/60‑1 (MQ=255)
                                                 ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:2549819/60‑1 (MQ=255)
                                                 ggATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAagag                       >  1:2379663/1‑71 (MQ=255)
                                                  gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                                             <  1:648730/48‑1 (MQ=255)
                                                  gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                                             <  1:1253892/48‑1 (MQ=255)
                                                  gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                                             <  1:1874019/48‑1 (MQ=255)
                                                  gATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAgg                      <  1:1173309/71‑1 (MQ=255)
                                                    tcgtcgCCTTGCTGAGCCGTTACCCCACCAACAAGCTAATCCCAt                                              >  1:1478389/1‑45 (MQ=255)
                                                     cgtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGccc                   >  1:620598/1‑71 (MQ=255)
                                                      gtcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGccc                   >  1:1554427/1‑70 (MQ=255)
                                                       tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                                                    >  1:816819/1‑36 (MQ=255)
                                                       tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGa                 <  1:1551600/71‑1 (MQ=255)
                                                       tcgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGa                 <  1:484548/71‑1 (MQ=255)
                                                        cgCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATCCGACCTAGGCTCAt                                    >  1:1053903/1‑51 (MQ=2)
                                                        cgCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATCCGACCTAGGCTCAt                                    >  1:2969043/1‑51 (MQ=2)
                                                         gCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGagg               >  1:2661987/1‑69 (MQ=255)
                                                            ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc            >  1:3403957/1‑71 (MQ=255)
                                                             tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc           >  1:391732/1‑71 (MQ=255)
                                                             tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc           >  1:1870370/1‑71 (MQ=255)
                                                             tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc           >  1:3468082/1‑71 (MQ=255)
                                                              ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccc          >  1:805235/1‑71 (MQ=255)
                                                                  aGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct       <  1:1641785/70‑1 (MQ=40)
                                                                   gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctct     >  1:339240/1‑71 (MQ=255)
                                                                    ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:2339024/41‑1 (MQ=255)
                                                                    ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc                                  <  1:159787/41‑1 (MQ=255)
                                                                      gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg  >  1:1700882/1‑71 (MQ=255)
                                                                      gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg  <  1:1753927/71‑1 (MQ=39)
                                                                      gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg  <  1:2076713/71‑1 (MQ=39)
                                                                      gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg  <  1:246218/71‑1 (MQ=39)
                                                                      gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTg  <  1:1127441/71‑1 (MQ=39)
                                                                      gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAACAGGCCCGAGGGTCCCCCTCTTTg  <  1:2876121/71‑1 (MQ=255)
                                                                      |                                                                      
GACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTG  >  minE/1763510‑1763650

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: ‑