Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,632 C→A 9.9% L377M (CTG→ATG)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,6320CA 9.9% 186.4 / 4.9 60L377M (CTG→ATG) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base A (4/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

TTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCT  >  minE/189564‑189701
                                                                    |                                                                      
ttGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCt                                                                       >  1:2457588/1‑70 (MQ=255)
 tGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg                                                                      <  1:555725/70‑1 (MQ=255)
    gggAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATc                                                                   >  1:875236/1‑70 (MQ=255)
     gggaaaGTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCca                                                                 >  1:3431554/4‑71 (MQ=255)
     gggaaaGTGATTTACGAAACCCAGTCCACTCACAAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCca                                                                 >  1:1732390/4‑71 (MQ=255)
      gAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcac                                                                <  1:1232237/71‑1 (MQ=255)
       aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCcaca                                                               <  1:2900179/71‑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)
                   cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAg                                                   <  1:2076571/70‑1 (MQ=255)
                   cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAg                                                   >  1:1049595/1‑70 (MQ=255)
                   cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAg                                                   >  1:3454258/1‑70 (MQ=255)
                   cGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTAGCTGATCCACATTAAAGGCg                                                      <  1:2447914/68‑1 (MQ=255)
                      aaCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATg                                               <  1:2130320/71‑1 (MQ=255)
                       aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGa                                              >  1:1912756/1‑71 (MQ=255)
                       aCGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGa                                              >  1:2772430/1‑71 (MQ=255)
                        cGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                                      <  1:3825564/63‑1 (MQ=255)
                            aTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACg                                            >  1:3144947/1‑68 (MQ=255)
                             tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                          >  1:1795071/1‑69 (MQ=255)
                             tCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                          >  1:30969/1‑69 (MQ=255)
                                  ccACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGc                                     >  1:577684/1‑69 (MQ=255)
                                   cacaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATCACGAAGAGGCCtt                                  >  1:2614385/1‑71 (MQ=255)
                                    acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                 <  1:3878029/71‑1 (MQ=255)
                                     caAACTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTa                                <  1:106517/71‑1 (MQ=255)
                                     caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                 >  1:1440890/1‑70 (MQ=255)
                                     caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                 >  1:109014/1‑70 (MQ=255)
                                     caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCttt                                 >  1:342552/1‑70 (MQ=255)
                                      aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTaa                               >  1:3151622/1‑71 (MQ=255)
                                      aaacTGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTaa                               >  1:3327390/1‑71 (MQ=255)
                                      aaaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGa                                           >  1:3602599/1‑59 (MQ=255)
                                       aaaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAAc                              <  1:468205/71‑1 (MQ=255)
                                        aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACg                             <  1:2524304/71‑1 (MQ=255)
                                        aaTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACg                             <  1:1989331/71‑1 (MQ=255)
                                         aTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTat                                                           <  1:2838230/41‑2 (MQ=255)
                                          tGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCg                                                      >  1:489712/1‑45 (MQ=255)
                                          tGCTGGCGGCGTTATCGCAGGCTTCGCTGATACACATTAAAGGCGAGTATGACGAag                                         >  1:874904/1‑57 (MQ=255)
                                            cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAg                          <  1:1100066/70‑1 (MQ=255)
                                            cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAg                          <  1:949683/70‑1 (MQ=255)
                                            cTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAg                          <  1:1421299/70‑1 (MQ=255)
                                               gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTAAACGAAGAAACCTTTAACGAAGCCt                       <  1:1810438/70‑1 (MQ=9)
                                                  gcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTat                   >  1:1164167/1‑71 (MQ=255)
                                                     ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                          <  1:2776494/45‑1 (MQ=255)
                                                     ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGaa                                          <  1:1195601/45‑1 (MQ=255)
                                                     ttATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgat                <  1:3050058/71‑1 (MQ=255)
                                                       aTCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg               >  1:541136/1‑70 (MQ=255)
                                                        tCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGc              <  1:1846554/70‑1 (MQ=255)
                                                                  cGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGc                                     >  1:1171744/1‑37 (MQ=255)
                                                                  cGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc   >  1:1936249/1‑71 (MQ=255)
                                                                  cGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc   >  1:2429243/1‑71 (MQ=255)
                                                                   gCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCt  >  1:1620737/1‑71 (MQ=255)
                                                                    |                                                                      
TTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGAGTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCT  >  minE/189564‑189701

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: