Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,137,295 G→A 47.2% G192G (GGC→GGT truB ← tRNA pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,137,2950GA47.2% 6.1 / 60.9 72G192G (GGC→GGTtruBtRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base G (15/23);  new base A (18/16);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCATCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCC  >  minE/2137228‑2137361
                                                                   |                                                                  
tCACCATCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAg                                                                    >  1:786468/1‑68 (MQ=255)
tCACCATCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAg                                                                    >  1:1805295/1‑68 (MQ=255)
tCACCATCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAg                                                                    >  1:1575008/1‑68 (MQ=255)
      tCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGtt                                                             >  1:2445072/1‑69 (MQ=255)
      tCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGtt                                                             >  1:256026/1‑69 (MQ=255)
      tCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGtt                                                             >  1:3070135/1‑69 (MQ=255)
      tCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGtt                                                             >  1:436142/1‑69 (MQ=255)
      tCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGtt                                                             >  1:836113/1‑69 (MQ=255)
      tCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCAAAGCCGAGtt                                                             >  1:1754110/1‑69 (MQ=255)
      tCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGAAGGTAAATAACATGCGCACCACAGCCGAGtt                                                             >  1:1984346/1‑69 (MQ=255)
        cGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGttt                                                            >  1:1304741/1‑68 (MQ=255)
        cGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGttt                                                            >  1:1794710/1‑68 (MQ=255)
        cGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGttt                                                            >  1:1690766/1‑68 (MQ=255)
          ttCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGttttt                                                          <  1:962993/68‑1 (MQ=255)
          ttCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGttttt                                                          <  1:1597108/68‑1 (MQ=255)
          ttCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGctttt                                                          <  1:276418/68‑1 (MQ=255)
               ccGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCAccc                                                     <  1:2800029/68‑1 (MQ=255)
               ccGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCAccc                                                     <  1:1918714/68‑1 (MQ=255)
               ccGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCAccc                                                     <  1:2911770/68‑1 (MQ=255)
                      tttACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGtcatc                                             >  1:1122222/1‑69 (MQ=255)
                      tttACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGtcat                                              >  1:727160/1‑68 (MQ=255)
                         acttacCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGtcatca                                            <  1:1596268/67‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:650198/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:3233151/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:3381719/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:709390/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:2561199/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:746069/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:839282/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:3135169/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:940054/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:1031119/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:1195247/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:1215347/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:1852467/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:1573370/63‑1 (MQ=255)
                                 ccAGACGACGCAAGTAAATAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:1660680/63‑1 (MQ=255)
                                    gACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGa                               <  1:968483/69‑1 (MQ=255)
                                    gACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGa                               <  1:2423797/69‑1 (MQ=255)
                                       ggCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                            <  1:1082734/69‑1 (MQ=255)
                                        gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                            <  1:2261311/68‑1 (MQ=255)
                                        gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                            <  1:1746818/68‑1 (MQ=255)
                                        gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                           <  1:2173170/69‑1 (MQ=255)
                                        gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                           <  1:166751/69‑1 (MQ=255)
                                        gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                           <  1:2297205/69‑1 (MQ=255)
                                         cgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                           <  1:1760907/68‑1 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:1292105/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:749617/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:1099491/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:658112/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:1886243/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:582229/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:1876516/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:355264/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:2788442/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:3308389/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:1338693/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:3160136/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:1363783/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:1453502/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:3134899/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                >  1:194312/1‑69 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTCCGGATATAAGTCCCtttt                 >  1:611783/1‑68 (MQ=255)
                                                   tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCCTCAATGATAGTGCGGATAAAAGTGCCTTTTg                >  1:3351491/1‑69 (MQ=255)
                                                       aTGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAatga                                        <  1:3383302/41‑1 (MQ=19)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         >  1:991624/1‑69 (MQ=255)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         <  1:882267/69‑1 (MQ=255)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         <  1:793578/69‑1 (MQ=255)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         <  1:315293/69‑1 (MQ=255)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         >  1:3088524/1‑69 (MQ=255)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         <  1:2206041/69‑1 (MQ=255)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         <  1:1095463/69‑1 (MQ=255)
                                                          cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg         >  1:1065793/1‑69 (MQ=255)
                                                                 caGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtcc  <  1:2040328/69‑1 (MQ=255)
                                                                 caGCCGAGTTTTTCACCCAGGGCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtcc  <  1:2767693/69‑1 (MQ=255)
                                                                   |                                                                  
TCACCATCCGTTCAACCGGATATTTACTTACCGCCAGACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCC  >  minE/2137228‑2137361

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: