Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 623,375 C→T 20.0% 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,3750CT20.0% 65.1 / 10.9 30T19I (ACC→ATC) msbAfused lipid transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (11/13);  new base T (2/4);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGCAGCCAGCGATA  >  minE/623306‑623443
                                                                     |                                                                    
tGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAAc                                                                      >  1:339247/1‑70 (MQ=255)
   tGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCa                                                                    <  1:551646/69‑1 (MQ=255)
        ttttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCATTGCGcc                                                             >  1:89500/1‑71 (MQ=255)
          ttGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCt                                                            <  1:1962092/70‑1 (MQ=255)
                 cATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAg                                                     <  1:1459479/70‑1 (MQ=255)
                  aTAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCg                                                   >  1:94520/1‑71 (MQ=255)
                  aTAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGAGGCCAACCATTGCGCCTTTCAAAGCg                                                   >  1:2014065/1‑71 (MQ=255)
                     aCGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAATCa                                                                    >  1:551903/1‑51 (MQ=255)
                     aCGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGt                                                >  1:1068045/1‑71 (MQ=255)
                         cAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTg                                             <  1:614151/70‑1 (MQ=255)
                              aTCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCaa                                                       <  1:1269579/55‑1 (MQ=255)
                              aTCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCaa                                                       <  1:281611/55‑1 (MQ=255)
                                ctctctACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggc                                     >  1:1602604/1‑71 (MQ=255)
                                ctctctACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTg                                       >  1:1983111/1‑69 (MQ=255)
                                   tctACGTGGCAGACATTCCGCCGACTGTGGCCAACCa                                                                    >  1:1552309/1‑37 (MQ=255)
                                   tctACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTg                                       <  1:728598/66‑1 (MQ=255)
                                   tctACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTg                                       <  1:1893367/66‑1 (MQ=255)
                                     tACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcgggcg                                >  1:1445186/1‑71 (MQ=255)
                                     tACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcggg                                  <  1:1751334/69‑1 (MQ=255)
                                     tACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTggcggg                                  <  1:802822/69‑1 (MQ=255)
                                           gCAGACATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGt                          <  1:2131686/71‑1 (MQ=255)
                                           gCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCa                                                        <  1:2029314/41‑1 (MQ=255)
                                           gCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGt                          >  1:1655211/1‑71 (MQ=255)
                                                cATTCCGCCGACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATc                     <  1:86898/71‑1 (MQ=255)
                                                cATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTaa                       >  1:39191/1‑69 (MQ=255)
                                                   tCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCt                   >  1:58579/1‑70 (MQ=255)
                                                    ccgccgACTGTGGCCAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTc                  <  1:1493457/70‑1 (MQ=255)
                                                              tgGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGcagcc         <  1:1321516/69‑1 (MQ=255)
                                                                  cAATCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGcagccagc      <  1:172619/68‑1 (MQ=255)
                                                                   aaCCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGCAGCCAGCGATa  >  1:1637229/1‑71 (MQ=255)
                                                                    aCCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGcagccagc      <  1:660555/66‑1 (MQ=255)
                                                                     |                                                                    
TGCTGGTTTTTTGAATGCATAACGACAAAGATCTCTCTACGTGGCAGACATTCCGCCGACTGTGGCCAACCATTGCGCCTTTCAAAGCGGGTCTGATCGTGGCGGGCGTAGCGTTAATCCTCAACGCAGCCAGCGATA  >  minE/623306‑623443

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: