Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 623,375 C→T 67.6% T19I (ACC→ATC)  msbA → fused lipid transporter subunits and membrane component and ATP‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE623,3750CT67.6% 73.7 / 55.9 63T19I (ACC→ATC) msbAfused lipid transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (14/6);  new base T (30/12);  total (45/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGCAGCCAGC  >  minE/623305‑623439
                                                                      |                                                                
aTGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  >  1:1076211/1‑71 (MQ=255)
aTGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  >  1:1234952/1‑71 (MQ=255)
aTGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  >  1:4187498/1‑71 (MQ=255)
aTGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  >  1:1312927/1‑71 (MQ=255)
aTGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  <  1:1799587/71‑1 (MQ=255)
aTGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  >  1:893468/1‑71 (MQ=255)
 tGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAcc                                                                 <  1:1687833/71‑1 (MQ=255)
   cTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCAt                                                               >  1:29049/1‑71 (MQ=255)
   cTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCAt                                                               >  1:2842070/1‑71 (MQ=255)
   cTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCAt                                                               >  1:2216433/1‑71 (MQ=255)
          tttGAATGGATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCt                                                        >  1:2303158/1‑71 (MQ=255)
          tttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCt                                                        >  1:1578915/1‑71 (MQ=255)
          tttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCt                                                        >  1:700961/1‑71 (MQ=255)
          tttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCt                                                        >  1:564524/1‑71 (MQ=255)
          tttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCt                                                        >  1:346974/1‑71 (MQ=255)
          tttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCt                                                        >  1:3658607/1‑71 (MQ=255)
           ttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCt                                                        <  1:2527960/70‑1 (MQ=255)
                  cATAACGACAAAGATCTCTCTACGTGGGAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGc                                                >  1:3701355/1‑71 (MQ=255)
                  cATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAg                                                 >  1:3009296/1‑70 (MQ=255)
                  cATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGc                                                >  1:3587444/1‑71 (MQ=255)
                  cATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGc                                                >  1:3906961/1‑71 (MQ=255)
                  cATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGc                                                <  1:2674779/71‑1 (MQ=255)
                  cATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGc                                                >  1:808810/1‑71 (MQ=255)
                  aaTAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGc                                                >  1:1257852/2‑71 (MQ=255)
                    tAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCa                                                                <  1:1615868/53‑1 (MQ=255)
                      aCGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  >  1:3084136/1‑49 (MQ=255)
                      aCGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                  >  1:3078786/1‑49 (MQ=255)
                      aCGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTgcgc                                                          >  1:2676048/1‑57 (MQ=255)
                      aCGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTg                                                             >  1:3221269/1‑54 (MQ=255)
                      aCGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTg                                                             >  1:3640869/1‑54 (MQ=255)
                        gACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGc                                                >  1:4193730/1‑65 (MQ=255)
                        gACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCt                                          >  1:1719779/1‑71 (MQ=255)
                        gACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCATTTCAAAGCGGGTCt                                          >  1:3990834/1‑71 (MQ=255)
                            aaGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATc                                      >  1:98805/1‑71 (MQ=255)
                            aaGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATc                                      >  1:4074037/1‑71 (MQ=255)
                            aaGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATc                                      >  1:3458174/1‑71 (MQ=255)
                              gATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGt                                    <  1:2464907/71‑1 (MQ=255)
                                     ctACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTggc                                 <  1:954459/67‑1 (MQ=255)
                                      tACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTc                                                     >  1:4088482/1‑46 (MQ=255)
                                      tACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcgggcg                            >  1:2860147/1‑71 (MQ=255)
                                      tACGTGGCAGACATTCCGCCGACTGTGGACAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcgggcg                            >  1:316824/1‑71 (MQ=255)
                                         gTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAg                         <  1:3479223/71‑1 (MQ=255)
                                         gTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAg                         <  1:1159562/71‑1 (MQ=255)
                                         gTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAg                         >  1:4130883/1‑71 (MQ=255)
                                             cAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGtt                     <  1:1613261/71‑1 (MQ=255)
                                              aGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGtt                     <  1:2406333/70‑1 (MQ=255)
                                              aGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTa                    <  1:3965867/71‑1 (MQ=255)
                                              aGACATTCCGCCGACTGTGGCCAAACATTGCGCCTTTCAAAGCGGGTCTGATCGTgg                                  >  1:1737833/1‑57 (MQ=255)
                                                aCATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAAt                  >  1:3884939/1‑71 (MQ=255)
                                                aCATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAAt                  >  1:1517420/1‑71 (MQ=255)
                                                aCATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAAt                  >  1:1846560/1‑71 (MQ=255)
                                                aCATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAAt                  >  1:1847693/1‑71 (MQ=255)
                                                aCATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAAt                  >  1:785582/1‑71 (MQ=255)
                                                 cATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATc                 <  1:297253/71‑1 (MQ=255)
                                                  aTTCCGGCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATcc                >  1:540772/1‑71 (MQ=255)
                                                  aTTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATGGTGGCGGGCGTAGCGTTAATcc                >  1:2159995/1‑71 (MQ=255)
                                                  aTTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATcc                >  1:3010453/1‑71 (MQ=255)
                                                  aTTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAAt                  >  1:1688822/1‑69 (MQ=255)
                                                  aTTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAAt                  >  1:1693549/1‑69 (MQ=255)
                                                    tCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAg                                                 >  1:1677659/1‑36 (MQ=255)
                                                    tCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAg                                                 >  1:3263746/1‑36 (MQ=255)
                                                     ccgccgACTGTGGCCAATCATTGCGCCTTTCAAAg                                                 >  1:2061519/1‑35 (MQ=255)
                                                     ccgccgACTGTGGCCAATCATTGCGCCTTTCAAAg                                                 >  1:3741719/1‑35 (MQ=255)
                                                      cgccgACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCa             <  1:3518362/70‑1 (MQ=255)
                                                      cgccgACTGTGGCCAATCATTGCGCCTTTCAAAGCGAGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCaa            <  1:197374/71‑1 (MQ=255)
                                                       gccgACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATACTCaa            <  1:4197124/70‑1 (MQ=255)
                                                         cgACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGc         <  1:4221783/71‑1 (MQ=255)
                                                         cgACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGc         <  1:2152337/71‑1 (MQ=255)
                                                             tgtgGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGcagc      >  1:1628822/1‑70 (MQ=255)
                                                             tgtgGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGcagc      >  1:3399556/1‑70 (MQ=255)
                                                                tggcAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGcagccagc  <  1:3842937/68‑1 (MQ=255)
                                                                 gCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGcagccagc  <  1:2939662/70‑1 (MQ=255)
                                                                      |                                                                
ATGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGCAGCCAGC  >  minE/623305‑623439

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: