Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,619 A→C 9.4% L372F (TTA→TTC)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6190AC 9.4% 155.1 / 11.1 64L372F (TTA→TTC) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base C (4/2);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATG  >  minE/189551‑189688
                                                                    |                                                                      
aGCGGCGAGCTTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATc                                                                      >  1:31746/1‑71 (MQ=255)
 gcggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                                                                     <  1:1353903/71‑1 (MQ=255)
   ggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                                                                   <  1:3767439/71‑1 (MQ=255)
   ggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                                                                   <  1:168812/71‑1 (MQ=255)
   ggcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCa                                                                   <  1:3518502/71‑1 (MQ=255)
    gcgAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGCCGGCGTTATCGCa                                                                   <  1:2304605/70‑1 (MQ=255)
         cGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTc                                                             >  1:602715/1‑71 (MQ=255)
            gtTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCtt                                                              <  1:3125983/67‑1 (MQ=255)
             ttGCGGGAAAAGTGATCTTCGAAACGCAATCGACCTACAAAATGCTGGCGGCGTTATCGCAGGCt                                                               <  1:598359/65‑1 (MQ=255)
             ttGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCt                                                          >  1:2457588/1‑70 (MQ=255)
              tGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                         <  1:555725/70‑1 (MQ=255)
                 gggAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATc                                                      >  1:875236/1‑70 (MQ=255)
                  gggaaaGTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCca                                                    >  1:1732390/4‑71 (MQ=255)
                  gggaaaGTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCca                                                    >  1:3431554/4‑71 (MQ=255)
                   gAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcac                                                   <  1:1232237/71‑1 (MQ=255)
                    aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcaca                                                  <  1:2900179/71‑1 (MQ=255)
                     aaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                                                                     <  1:2430194/51‑1 (MQ=255)
                      aaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACAt                                                 >  1:3542458/1‑70 (MQ=255)
                       aGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                               <  1:1927645/71‑1 (MQ=255)
                        gTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                               >  1:2282388/1‑70 (MQ=255)
                           aTCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTaaa                                             <  1:3768447/69‑1 (MQ=255)
                           aTCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTa                                               <  1:238895/67‑1 (MQ=255)
                              ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                         <  1:2168222/54‑1 (MQ=255)
                              ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAgg                                           >  1:2955028/1‑68 (MQ=255)
                              ttCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                          >  1:754662/1‑69 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTTATCCACATTAAAGGc                                          <  1:3926652/67‑1 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                          <  1:1068284/67‑1 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGc                                          <  1:375758/67‑1 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                         <  1:1811161/68‑1 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                      <  1:2076571/70‑1 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                      >  1:1049595/1‑70 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑g                                      >  1:3454258/1‑70 (MQ=255)
                                cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTAGCTGATCCACATTAAAGGCg                                         <  1:2447914/68‑1 (MQ=255)
                                   aaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATg                                  <  1:2130320/71‑1 (MQ=255)
                                    aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                 >  1:2772430/1‑71 (MQ=255)
                                    aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGa                                 >  1:1912756/1‑71 (MQ=255)
                                     cGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                         <  1:3825564/63‑1 (MQ=255)
                                         aTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACg                               >  1:3144947/1‑68 (MQ=255)
                                          tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                             >  1:30969/1‑69 (MQ=255)
                                          tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                             >  1:1795071/1‑69 (MQ=255)
                                               ccACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGc                        >  1:577684/1‑69 (MQ=255)
                                                cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATCACGAAGAGGCCtt                     >  1:2614385/1‑71 (MQ=255)
                                                 acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                    <  1:3878029/71‑1 (MQ=255)
                                                  caAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTa                   <  1:106517/71‑1 (MQ=255)
                                                  caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                    >  1:342552/1‑70 (MQ=255)
                                                  caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                    >  1:1440890/1‑70 (MQ=255)
                                                  caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                    >  1:109014/1‑70 (MQ=255)
                                                   aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTaa                  >  1:3327390/1‑71 (MQ=255)
                                                   aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTaa                  >  1:3151622/1‑71 (MQ=255)
                                                   aaaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGa                              >  1:3602599/1‑59 (MQ=255)
                                                    aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAAc                 <  1:468205/71‑1 (MQ=255)
                                                     aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                <  1:1989331/71‑1 (MQ=255)
                                                     aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                <  1:2524304/71‑1 (MQ=255)
                                                      aTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTat                                              <  1:2838230/41‑2 (MQ=255)
                                                       tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                         >  1:489712/1‑45 (MQ=255)
                                                       tGCTGGCGGCGTTATCGCAGGCTTCGCTGATACACATTAAAGGCGA‑GTATGACGAag                            >  1:874904/1‑57 (MQ=255)
                                                         cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAg             <  1:1421299/70‑1 (MQ=255)
                                                         cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAg             <  1:1100066/70‑1 (MQ=255)
                                                         cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAg             <  1:949683/70‑1 (MQ=255)
                                                            gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCt          <  1:1810438/70‑1 (MQ=9)
                                                               gcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat      >  1:1164167/1‑71 (MQ=255)
                                                                  ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                             <  1:2776494/45‑1 (MQ=255)
                                                                  ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGaa                             <  1:1195601/45‑1 (MQ=255)
                                                                  ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgat   <  1:3050058/71‑1 (MQ=255)
                                                                    aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg  >  1:541136/1‑70 (MQ=255)
                                                                    |                                                                      
AGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATG  >  minE/189551‑189688

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: ‑