Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,509,560 T→G 18.8% T486T (ACA→ACC yicH ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,509,5600TG18.8% 230.0 / 49.6 133T486T (ACA→ACCyicHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/66);  new base G (11/14);  total (53/80)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

AATTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGAATGCCTTTTTCAGTAAATGCA  >  minE/2509495‑2509626
                                                                 |                                                                  
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATgtgtgt                                                                  >  1:1263453/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:2244295/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:1267081/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:605041/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:1886344/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:917017/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:2222269/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:1881944/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:721272/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:2908459/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:1704060/1‑68 (MQ=255)
aaTTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:1514713/1‑68 (MQ=255)
aaTTAAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATgtgt                                                                  >  1:1003839/1‑68 (MQ=255)
    cAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGAc                                                             >  1:2924292/1‑69 (MQ=255)
    cAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGAc                                                             >  1:1935922/1‑69 (MQ=255)
    cAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGAc                                                             >  1:51919/1‑69 (MQ=255)
    cAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGAc                                                             >  1:1453734/1‑69 (MQ=255)
    cAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGAc                                                             >  1:84779/1‑69 (MQ=255)
    cAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGAc                                                             >  1:513542/1‑69 (MQ=255)
    cAGGCCATCCCCACTGTTGCAAAATATTCAACGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGAc                                                             >  1:3342012/1‑69 (MQ=255)
       gCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGt                                                           <  1:1560419/68‑1 (MQ=255)
        caaTCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGt                                                           <  1:3081879/65‑1 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:2369138/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:2022314/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:492756/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:228959/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:2493625/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:787794/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:3184010/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:1495602/1‑69 (MQ=255)
             cccACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCg                                                    >  1:2732277/1‑69 (MQ=255)
                 atGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGATGTTGTGGCGttt                                                 >  1:3255172/2‑68 (MQ=255)
                 atGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGttt                                                 >  1:864724/2‑68 (MQ=255)
                 atGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGttt                                                 >  1:1199575/2‑68 (MQ=255)
                 atGTTGAAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGttt                                                 >  1:4138/2‑68 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCTTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:2630070/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTtgt                                               >  1:2667773/1‑68 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:2949486/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:2962923/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:419121/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:2919634/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:1838177/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:1824811/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:3229247/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:3372142/69‑1 (MQ=255)
                   gTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:2689269/69‑1 (MQ=255)
                      gCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTTATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:593004/69‑1 (MQ=255)
                      gCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:1069892/69‑1 (MQ=255)
                      gCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:3379978/69‑1 (MQ=255)
                      gCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:2560994/69‑1 (MQ=255)
                      gCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:1504062/69‑1 (MQ=255)
                      gCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:2250094/69‑1 (MQ=255)
                      gCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:1346119/69‑1 (MQ=255)
                       cAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGaaa                                           <  1:2356370/68‑1 (MQ=255)
                        aaaaTATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTTTGAAac                                          <  1:2935758/68‑1 (MQ=255)
                        aaaaTATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAac                                          <  1:2296990/68‑1 (MQ=255)
                        aaaaTATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAac                                          <  1:3280654/68‑1 (MQ=255)
                            tatTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTgg                                     >  1:661710/1‑69 (MQ=255)
                                   cGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTgg                               <  1:3006862/68‑1 (MQ=255)
                                   cGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTgg                               <  1:3230539/68‑1 (MQ=255)
                                   cGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTgg                               <  1:1817129/68‑1 (MQ=255)
                                   cGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTgg                               <  1:363044/68‑1 (MQ=255)
                                   cGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTgg                               <  1:1344697/68‑1 (MQ=255)
                                   cGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTgg                               <  1:2428438/68‑1 (MQ=255)
                                   cGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTgg                               <  1:3397328/68‑1 (MQ=255)
                                     gCACACCGCGTCCATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:3224852/51‑1 (MQ=255)
                                         acCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:492887/69‑1 (MQ=255)
                                           cgcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:1473625/1‑67 (MQ=255)
                                           cgcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:2843911/1‑67 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGAGGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:3207567/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:758834/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:685025/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:624125/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:962371/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:1198181/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:3341114/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:2847987/66‑1 (MQ=255)
                                            gcgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:2546813/66‑1 (MQ=255)
                                             cgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:1387780/65‑1 (MQ=255)
                                             cgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:2647948/65‑1 (MQ=255)
                                             cgTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        <  1:1962845/65‑1 (MQ=255)
                                              gTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGTCGGTGGCTTCCAg                        <  1:2425307/64‑1 (MQ=255)
                                                 cATTCAGGCTGATATGGGTCTGACGTTGTGGCGTTtgtg                                              <  1:384791/39‑1 (MQ=255)
                                                    tCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:3075111/1‑58 (MQ=255)
                                                    tCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:1909148/1‑58 (MQ=255)
                                                    tCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:2639192/1‑58 (MQ=255)
                                                    tCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:3428918/1‑58 (MQ=255)
                                                    tCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:2054799/1‑58 (MQ=255)
                                                    tCAGGCTGATATGTGTCTGACGTTGTGGAGTTTGTGAAACACTGGCGGTGGCTTCCAg                        >  1:94034/1‑58 (MQ=255)
                                                    acAGGCTGATATGTGTCTGAAGTTGTGGCGTTTGTGAAac                                          >  1:2210016/2‑40 (MQ=39)
                                                     cAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGaa                      <  1:2108364/59‑1 (MQ=255)
                                                     cAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGaa                      <  1:2984832/59‑1 (MQ=255)
                                                     cAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGaa                      <  1:1817153/59‑1 (MQ=255)
                                                     cAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGaa                      <  1:871502/59‑1 (MQ=255)
                                                     cAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGaa                      <  1:876631/59‑1 (MQ=255)
                                                     cAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGaa                      <  1:348474/59‑1 (MQ=255)
                                                     cAGGCTGATATGTGTCTGACGTCGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGaa                      <  1:3113774/59‑1 (MQ=255)
                                                         cTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGAATGCCTTTTTCAGTa        <  1:376820/69‑1 (MQ=255)
                                                          tGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGAATGCCTTTTTCAGTa        <  1:1437120/68‑1 (MQ=255)
                                                                 tgtCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGAATGCCTTTTTCAGTAAATGCa  >  1:1443915/1‑67 (MQ=255)
                                                                 |                                                                  
AATTCAGGCCATCCCCACTGTTGCAAAATATTCACCGGCACACCGCGTCCATTCAGGCTGATATGTGTCTGACGTTGTGGCGTTTGTGAAACACTGGCGGTGGCTTCCAGAATGCCTTTTTCAGTAAATGCA  >  minE/2509495‑2509626

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: