Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 228,490 C→G 10.1% L71L (CTG→CTC phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE228,4900CG10.1% 284.5 / 16.9 99L71L (CTG→CTCphoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base C (50/39);  new base G (6/4);  total (56/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 113 total aligned reads displayed.

TTTTGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTG  >  minE/228422‑228559
                                                                    |                                                                     
ttttGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:3875087/1‑71 (MQ=255)
ttttGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:5551624/1‑71 (MQ=255)
ttttGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:2439556/1‑71 (MQ=255)
  ttGCTGTTCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:657787/1‑69 (MQ=255)
  ttGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:4049169/1‑69 (MQ=255)
  ttGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:3292966/1‑69 (MQ=255)
  ttGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGtt                                                                   >  1:4409006/1‑71 (MQ=255)
            gTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAAt                                                         <  1:2155505/71‑1 (MQ=255)
               tCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtc                                                                       <  1:842339/54‑1 (MQ=255)
               tCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtc                                                                       <  1:445216/54‑1 (MQ=255)
               tCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTtg                                                      >  1:4342868/1‑71 (MQ=255)
                 aCTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:3976261/1‑54 (MQ=255)
                 aCTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTtgtg                                                    >  1:1540829/1‑71 (MQ=255)
                 aCTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTtgt                                                     <  1:2981487/70‑1 (MQ=255)
                 aCTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTtgt                                                     <  1:3623989/70‑1 (MQ=255)
                 aCTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTtgt                                                     <  1:527229/70‑1 (MQ=255)
                  ctctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                   <  1:4934302/71‑1 (MQ=255)
                  ctctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                   <  1:4764004/71‑1 (MQ=255)
                  ctctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                   <  1:1072710/71‑1 (MQ=255)
                  ctctctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                   <  1:5472869/71‑1 (MQ=255)
                     tctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:795136/1‑50 (MQ=255)
                     tctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAAt                                                         >  1:5385966/1‑62 (MQ=255)
                     tctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGtt                                                  <  1:3227798/69‑1 (MQ=255)
                     tctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTc                                                >  1:1653989/1‑71 (MQ=255)
                     tctGCTTTATTACAGGCAAACTCTGCTTCCCAACGACCATAACcagtcag                                                                     >  1:4129171/1‑50 (MQ=255)
                      ctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCg                                               >  1:4096218/1‑71 (MQ=255)
                      ctGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCg                                               >  1:2094000/1‑71 (MQ=255)
                       tGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAAtt                                                        >  1:5108429/1‑61 (MQ=255)
                       tGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAAtt                                                        >  1:281934/1‑61 (MQ=255)
                       tGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAAtt                                                        >  1:1156845/1‑61 (MQ=255)
                       tGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTc                                                <  1:3628948/69‑1 (MQ=255)
                       tGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTc                                                <  1:3976426/69‑1 (MQ=255)
                           ttattaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                 >  1:1137187/1‑64 (MQ=255)
                             attaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                 >  1:5189293/1‑62 (MQ=255)
                             attaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                 >  1:5039094/1‑62 (MQ=255)
                             attaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                 >  1:1941681/1‑62 (MQ=255)
                              ttaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                        <  1:4666553/70‑1 (MQ=255)
                              ttaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                        <  1:2981255/70‑1 (MQ=255)
                              ttaCCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                        >  1:2897838/1‑70 (MQ=255)
                                  cGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCa                                   <  1:1633794/71‑1 (MQ=255)
                                   ggCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                  <  1:909865/71‑1 (MQ=255)
                                    gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGac                                                                >  1:3582552/1‑39 (MQ=255)
                                    gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAAt                                                         <  1:3207611/47‑1 (MQ=255)
                                    gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAAt                                                         <  1:1964778/47‑1 (MQ=255)
                                    gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                   >  1:2445531/1‑53 (MQ=255)
                                    gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                  >  1:4504005/1‑70 (MQ=255)
                                    gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                  >  1:1563322/1‑70 (MQ=255)
                                    gCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCa                                   >  1:3912527/1‑69 (MQ=255)
                                      aaaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGtt                                                  >  1:5351047/1‑52 (MQ=255)
                                      aaaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                        >  1:4154855/1‑62 (MQ=255)
                                      aaaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaa                                  <  1:1368592/68‑1 (MQ=255)
                                      aaaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAAc                               >  1:2632157/1‑71 (MQ=255)
                                      aaaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAAc                               >  1:2625588/1‑71 (MQ=255)
                                      aaaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGAGTTTCGCCTTTGa                                        >  1:1370057/1‑62 (MQ=255)
                                       aaCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACg                              <  1:4491906/71‑1 (MQ=255)
                                        aCTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACg                              >  1:4411069/1‑70 (MQ=255)
                                         ctctGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGa                            >  1:528606/1‑71 (MQ=255)
                                         ctctGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGa                            >  1:4827552/1‑71 (MQ=255)
                                          tctGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGa                            >  1:4073817/1‑70 (MQ=255)
                                           ctGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGata                          >  1:1964451/1‑71 (MQ=255)
                                            tGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGata                          >  1:5544841/1‑70 (MQ=255)
                                             gCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGatat                         >  1:5195967/1‑70 (MQ=255)
                                              cTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGatat                         <  1:3863308/69‑1 (MQ=255)
                                              cTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATaa                       >  1:4365819/1‑71 (MQ=255)
                                               ttCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCt                                            <  1:5085107/49‑1 (MQ=255)
                                               ttCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCt                                            <  1:4642093/49‑1 (MQ=255)
                                               ttCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCt                                            <  1:2792128/49‑1 (MQ=255)
                                               ttCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAAc                      >  1:4191204/1‑71 (MQ=255)
                                                tCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCttt                                          <  1:4294114/50‑1 (MQ=255)
                                                tCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGa                                        <  1:3324083/52‑1 (MQ=255)
                                                 cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                 >  1:3809329/1‑42 (MQ=255)
                                                 cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                 >  1:3966962/1‑42 (MQ=255)
                                                 cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGttt                                                 >  1:4288110/1‑42 (MQ=255)
                                                 cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaaaa                                >  1:5209169/1‑59 (MQ=255)
                                                 cccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCaaaa                                >  1:5140426/1‑59 (MQ=255)
                                                  ccAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAActc                    <  1:154052/70‑1 (MQ=255)
                                                   cAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTc                                                <  1:1548392/41‑1 (MQ=255)
                                                   cAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTc                                                <  1:2301464/41‑1 (MQ=255)
                                                   cAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCttt                                          <  1:188364/47‑1 (MQ=255)
                                                   cAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTg                  <  1:859528/71‑1 (MQ=255)
                                                    attgACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTgg                 >  1:2370112/4‑71 (MQ=12)
                                                    aaCGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCt                                            >  1:1296081/1‑44 (MQ=255)
                                                     aCGACCATAACCAGTCAGTTGATCGTTAATTTgtgt                                                   <  1:1497812/36‑1 (MQ=255)
                                                     aCGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTggg                <  1:3045975/71‑2 (MQ=255)
                                                      cgcccATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTc               >  1:2905969/4‑71 (MQ=255)
                                                       gACCGTAGCCAGTGAGTTGATCGTTAATTTGGGttt                                                 >  1:1636650/1‑36 (MQ=19)
                                                       gACCGTAGCCAGTGAGTTGATCGTTAATTTGGGttt                                                 >  1:2848781/1‑36 (MQ=19)
                                                       gACCGTAGCCAGTGAGTTGATCGTTAATTTGGGttt                                                 >  1:4597878/1‑36 (MQ=19)
                                                       gACCGTAGCCAGTGAGTTGATCGTTAATTTGGGttt                                                 >  1:2383779/1‑36 (MQ=19)
                                                        aCCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGc             <  1:703582/71‑1 (MQ=255)
                                                          cATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACgg                             <  1:1167005/53‑1 (MQ=255)
                                                          cATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCa           >  1:4700275/1‑71 (MQ=255)
                                                          cATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCa           >  1:5387717/1‑71 (MQ=255)
                                                           aTAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGt                >  1:1106622/1‑65 (MQ=255)
                                                           aTAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGt                >  1:4397428/1‑65 (MQ=255)
                                                            tAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATc         <  1:5620118/71‑1 (MQ=255)
                                                            tAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATc         <  1:2667654/71‑1 (MQ=255)
                                                              aCCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCt        <  1:3673667/70‑1 (MQ=255)
                                                               ccAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCttt      <  1:2156731/71‑1 (MQ=16)
                                                                  gtcagtTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACt   >  1:4747545/1‑71 (MQ=255)
                                                                    |                                                                     
TTTTGCTGTGCAGTATCACTCTCTGCTTTATTACCGGCAAACTCTGCTTCCCAACGACCATAACCAGTCAGTTGATCGTTAATTTGTGTTTCGCCTTTGAAACCAAAACGGATATAACTCTGGTCGCCATCTTTACTG  >  minE/228422‑228559

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: