Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,277,603 G→A 42.6% N382N (AAC→AAT 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,6030GA42.6% 22.7 / 48.8 68N382N (AAC→AATbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (11/28);  new base A (7/22);  total (18/50)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGA  >  minE/2277537‑2277670
                                                                  |                                                                   
gCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGtt                                                                   <  1:1156755/69‑1 (MQ=255)
gCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGtt                                                                   <  1:2910390/69‑1 (MQ=255)
   aTAACCTTCGATGAATTCCCAACCGGCTCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATc                                                                >  1:1156907/1‑69 (MQ=255)
   aTAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGtt                                                                   <  1:1718846/66‑1 (MQ=255)
   aTAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATc                                                                >  1:2518156/1‑69 (MQ=255)
   aTAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATc                                                                <  1:2015337/69‑1 (MQ=255)
     taCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:175132/68‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:2322560/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              >  1:1208041/1‑69 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:96220/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:189982/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:3136614/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:2272650/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:1822995/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:176609/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:376700/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:1692704/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:163970/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:49696/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:515162/69‑1 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              >  1:636213/1‑69 (MQ=255)
     aaCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGt                                                              <  1:747021/69‑1 (MQ=255)
      aCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTc                                                             >  1:1841791/1‑69 (MQ=255)
      aCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTc                                                             >  1:807695/1‑69 (MQ=255)
      aCCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTc                                                             >  1:2017643/1‑69 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:2677232/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:994301/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:2675729/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:2805733/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:2830014/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:315295/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:3189382/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:2464005/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:835301/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:407345/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:621038/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:1224927/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:1234133/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:220802/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:2136951/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:1899679/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:1302995/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:1394548/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:142466/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:1471330/47‑1 (MQ=255)
                                  gTTTGCCAGGTACCATGACGACCAAACTGAGAATTATCGTCGCTGCg                                                       <  1:1518640/47‑1 (MQ=255)
                                            ttCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACtt                       <  1:835398/69‑1 (MQ=255)
                                             tCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACttg                      >  1:1716949/1‑69 (MQ=255)
                                             tCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACtt                       <  1:361020/68‑1 (MQ=255)
                                              ccATGACGACCAAACTGAGAATTATCTTCGCTGCg                                                       <  1:510142/35‑1 (MQ=19)
                                                   acgacCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAg                >  1:2220327/1‑69 (MQ=255)
                                                    cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                >  1:2652026/1‑68 (MQ=255)
                                                    cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                >  1:442789/1‑68 (MQ=255)
                                                    cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                >  1:1140641/1‑68 (MQ=255)
                                                    cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                >  1:358266/1‑68 (MQ=255)
                                                    cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                >  1:1898742/1‑68 (MQ=255)
                                                    cgacCAAACTGAGAATTATCGTCGCTGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTCGCTGCAg                >  1:387305/1‑68 (MQ=255)
                                                    cgacCAAACTGAGAATTATCGTCGCAGCGTCCTGCGCCTTCAAAAGTAAAATCGCCGACTTGCTGCAg                >  1:861381/1‑68 (MQ=255)
                                                           aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa        <  1:416026/69‑1 (MQ=255)
                                                           aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa        <  1:1543739/69‑1 (MQ=255)
                                                           aCTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCa        <  1:323795/69‑1 (MQ=255)
                                                             tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga      <  1:1669890/69‑1 (MQ=255)
                                                             tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga      <  1:3300946/69‑1 (MQ=255)
                                                             tgtgAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCaga      <  1:1835946/69‑1 (MQ=255)
                                                                 aGTTATCGTCACTGCTTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagata    <  1:2651043/67‑1 (MQ=255)
                                                                 aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga  >  1:2376560/1‑69 (MQ=255)
                                                                 aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagataga  >  1:1108553/1‑69 (MQ=255)
                                                                 aGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCagata    <  1:2777580/67‑1 (MQ=255)
                                                                  |                                                                   
GCGATAACCTTCGATGAATTCCCAACCGGCGCTGGTTTGCCAGGTTCCATGACGACCAAACTGTGAGTTATCGTCACTGCGTGCTGCGCCTTCAAAGGTAAAATCGCCGACTTGTTGCAGCCCGGTCAGATAGA  >  minE/2277537‑2277670

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: