Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,277,633 G→A 20.1% T372T (ACC→ACT btuB ← vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,277,6330GA20.1% 93.8 / 14.0 64T372T (ACC→ACTbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (35/16);  new base A (7/6);  total (42/22)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATC  >  minE/2277581‑2277694
                                                    |                                                             
ttCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACtt                                               <  1:835398/69‑1 (MQ=255)
 tCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACttg                                              >  1:1716949/1‑69 (MQ=255)
 tCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACtt                                               <  1:361020/68‑1 (MQ=255)
       acgacCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAg                                        >  1:2220327/1‑69 (MQ=255)
        cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                        >  1:1898742/1‑68 (MQ=255)
        cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                        >  1:442789/1‑68 (MQ=255)
        cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                        >  1:358266/1‑68 (MQ=255)
        cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                        >  1:2652026/1‑68 (MQ=255)
        cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                        >  1:1140641/1‑68 (MQ=255)
        cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTCGCTGCAg                                        >  1:387305/1‑68 (MQ=255)
        cgacCAAACTGAGAATTATCGTCGCAGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                                        >  1:861381/1‑68 (MQ=255)
               aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                <  1:1543739/69‑1 (MQ=255)
               aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                <  1:416026/69‑1 (MQ=255)
               aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                <  1:323795/69‑1 (MQ=255)
                 tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                              <  1:3300946/69‑1 (MQ=255)
                 tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                              <  1:1669890/69‑1 (MQ=255)
                 tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                              <  1:1835946/69‑1 (MQ=255)
                     aGTTATCGTCACTGCTTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagata                            <  1:2651043/67‑1 (MQ=255)
                     aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga                          >  1:1108553/1‑69 (MQ=255)
                     aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga                          >  1:2376560/1‑69 (MQ=255)
                     aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagata                            <  1:2777580/67‑1 (MQ=255)
                      attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                         <  1:1994254/68‑1 (MQ=255)
                      attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                         <  1:1826605/68‑1 (MQ=255)
                      attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                         <  1:2314032/68‑1 (MQ=255)
                      attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                         <  1:1743977/68‑1 (MQ=255)
                      attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                         <  1:448117/68‑1 (MQ=255)
                      attATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAGCCCAGTCagatagat                         <  1:1990297/68‑1 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTTAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:1224645/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                >  1:3329377/1‑59 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                >  1:2820658/1‑59 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                >  1:2971345/1‑59 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                >  1:331614/1‑59 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                >  1:1307019/1‑59 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa                                >  1:1497401/1‑59 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATGGATGcc                      >  1:1805626/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:447381/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:963159/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:385943/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:2651802/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:2764923/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:1087534/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:1287081/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:1364270/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:1765728/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:2409963/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:2566580/1‑69 (MQ=255)
                         aTCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:2783055/1‑69 (MQ=255)
                         aTCGTCACTGCGTACTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGcc                      >  1:28792/1‑69 (MQ=255)
                            gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCgg                    <  1:1088309/68‑1 (MQ=255)
                            gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                   >  1:1008764/1‑69 (MQ=255)
                            gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                   >  1:59919/1‑69 (MQ=255)
                            gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                   >  1:2892193/1‑69 (MQ=255)
                            gTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGt                   >  1:2039455/1‑69 (MQ=255)
                                   cGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACg             >  1:501596/1‑68 (MQ=255)
                                   cGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACg             >  1:671523/1‑68 (MQ=255)
                                     tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCTGTCagataga                          >  1:1042312/1‑53 (MQ=255)
                                     tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga                          >  1:2663400/1‑53 (MQ=255)
                                     tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagatag                           >  1:2695397/1‑52 (MQ=255)
                                     tgctgcGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga                              >  1:1720942/1‑49 (MQ=255)
                                          cgcCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTAt                 <  1:2522180/57‑1 (MQ=255)
                                             cTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc  <  1:598140/69‑1 (MQ=255)
                                             cTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc  <  1:2550816/69‑1 (MQ=255)
                                              ttCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc  <  1:2905848/68‑1 (MQ=255)
                                               tCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATatc  <  1:2676456/67‑1 (MQ=255)
                                                    |                                                             
TTCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGATGCCGGTATTACGTTGATCATATC  >  minE/2277581‑2277694

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: