Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 972,048 G→A 6.4% P285P (CCG→CCA)  lsrC → AI2 transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE972,0480GA 6.4% 253.2 / 5.5 78P285P (CCG→CCA) lsrCAI2 transporter
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base A (2/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGC  >  minE/971978‑972111
                                                                      |                                                               
tCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCg                                                                 >  1:479295/1‑71 (MQ=255)
 cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCgg                                                                >  1:1949145/1‑71 (MQ=255)
 cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCgg                                                                >  1:3391615/1‑71 (MQ=255)
 cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCgg                                                                >  1:3104327/1‑71 (MQ=255)
 cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCgg                                                                >  1:562453/1‑71 (MQ=255)
 cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCg                                                                 >  1:1794897/1‑70 (MQ=255)
 cATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCg                                                                 >  1:480017/1‑70 (MQ=255)
  aTTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCgg                                                                >  1:2151092/1‑70 (MQ=255)
   ttGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCa                                                              <  1:3216873/71‑1 (MQ=255)
    tGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCAt                                                             >  1:176187/1‑71 (MQ=255)
    tGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCAt                                                             >  1:773392/1‑71 (MQ=255)
      gTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCAtgg                                                           <  1:1274282/71‑1 (MQ=255)
      gTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCAtg                                                            <  1:3049517/70‑1 (MQ=255)
          ggTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGa                                                       <  1:3109060/71‑1 (MQ=255)
           gTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCACCATGCCGGaa                                                      <  1:1532652/71‑1 (MQ=255)
           gTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCACCATGCCGGaa                                                      <  1:2190330/71‑1 (MQ=255)
           gTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCACCATGCCGGaa                                                      <  1:2346588/71‑1 (MQ=255)
            tACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGaa                                                      <  1:1347220/70‑1 (MQ=255)
              cTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGa                                                   >  1:115932/1‑71 (MQ=255)
              cTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGa                                                   >  1:53865/1‑71 (MQ=255)
              cTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGa                                                   >  1:2922569/1‑71 (MQ=255)
              cTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGa                                                   <  1:2603140/71‑1 (MQ=255)
               tCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAt                                                  >  1:2351890/1‑71 (MQ=255)
               tCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAt                                                  >  1:2558054/1‑71 (MQ=255)
                cGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtt                                                 >  1:522407/1‑71 (MQ=255)
                 ggCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAttt                                                >  1:1711297/1‑71 (MQ=255)
                 ggCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtt                                                 <  1:1070939/70‑1 (MQ=255)
                 ggCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtt                                                 <  1:1042739/70‑1 (MQ=255)
                 ggCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtt                                                 <  1:2462510/70‑1 (MQ=255)
                   cgcATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTa                                              >  1:1421939/1‑71 (MQ=255)
                   cgcATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTa                                              >  1:1303960/1‑71 (MQ=255)
                    gcATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTa                                              >  1:1388405/1‑70 (MQ=255)
                    gcATGGATCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTAt                                             >  1:475189/1‑71 (MQ=255)
                      aTGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcg                                           >  1:2252203/1‑71 (MQ=255)
                      aTGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATc                                            >  1:854239/1‑70 (MQ=255)
                        ggTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCGATCATGGTGGa                                                       >  1:1344434/1‑57 (MQ=255)
                        ggTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCGATCATGGTGGa                                                       >  1:3141011/1‑57 (MQ=255)
                        ggTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcgcg                                         >  1:1783334/1‑71 (MQ=255)
                        ggTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcgcg                                         >  1:2772366/1‑71 (MQ=255)
                         gTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATcgcg                                         <  1:1722728/70‑1 (MQ=255)
                          ttCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCggg                                       >  1:791379/1‑71 (MQ=255)
                            ccTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGaa                                                      <  1:2625394/54‑1 (MQ=255)
                                aCGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGt                                 >  1:2418447/1‑71 (MQ=255)
                                aCGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGt                                 >  1:1375880/1‑71 (MQ=255)
                                aCGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGt                                 >  1:2921700/1‑71 (MQ=255)
                                aCGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGt                                 >  1:2931175/1‑71 (MQ=255)
                                 cGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGtt                                >  1:520307/1‑71 (MQ=255)
                                    aGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtttt                                               <  1:3350699/53‑1 (MQ=255)
                                    aGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtttt                                               <  1:2152831/53‑1 (MQ=255)
                                    aGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTctg                             >  1:1200757/1‑71 (MQ=255)
                                      aTCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtttt                                               <  1:3134441/51‑1 (MQ=255)
                                      aTCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtttt                                               <  1:2502896/51‑1 (MQ=255)
                                      aTCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtttt                                               <  1:3622045/51‑1 (MQ=255)
                                      aTCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTAt                                             >  1:2432318/1‑53 (MQ=255)
                                      aTCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTAt                                             >  1:2094076/1‑53 (MQ=255)
                                      aTCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTctgct                           <  1:1692736/71‑1 (MQ=255)
                                        cGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTctgctg                          >  1:3150752/1‑70 (MQ=255)
                                        cGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                         >  1:343339/1‑71 (MQ=255)
                                        cGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                         >  1:3418689/1‑71 (MQ=255)
                                        cGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                         >  1:3478225/1‑71 (MQ=255)
                                        cGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                         >  1:2521575/1‑71 (MQ=255)
                                        cGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                         >  1:324693/1‑71 (MQ=255)
                                        cGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTgg                         >  1:2209732/1‑71 (MQ=255)
                                         gATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGc                        <  1:3630133/71‑1 (MQ=255)
                                          aTAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCg                       <  1:1220915/71‑1 (MQ=255)
                                            aGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGt                     >  1:51672/1‑71 (MQ=255)
                                             gCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGt                     <  1:1369369/70‑1 (MQ=255)
                                              cGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGc                   >  1:2731814/1‑71 (MQ=255)
                                                  cTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGgt               >  1:3489280/1‑71 (MQ=255)
                                                  cTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGgt               >  1:821924/1‑71 (MQ=255)
                                                   tGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtt                                                 <  1:1077170/36‑1 (MQ=255)
                                                   tGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAtt                                                 <  1:304934/36‑1 (MQ=255)
                                                   tGGTGCTGTTGCGCATTCCGGCATGGTGGAATGAt                                                  <  1:947453/35‑1 (MQ=255)
                                                     gTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGgtg              <  1:2491698/69‑1 (MQ=255)
                                                     gTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGgtg              <  1:535337/69‑1 (MQ=255)
                                                      tGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTctgc                            <  1:2231473/54‑1 (MQ=255)
                                                      tGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGttt           <  1:2039131/71‑1 (MQ=255)
                                                       gCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGttt           >  1:1739532/1‑70 (MQ=255)
                                                       gCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGttt           >  1:681068/1‑70 (MQ=255)
                                                        cTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTg          <  1:1586071/70‑1 (MQ=255)
                                                           ttGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATgg      <  1:3458378/71‑1 (MQ=255)
                                                           ttGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATg       <  1:657438/70‑1 (MQ=255)
                                                           ttGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATg       >  1:3160779/1‑70 (MQ=255)
                                                           ttGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATg       <  1:962186/70‑1 (MQ=255)
                                                               gcATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGc  <  1:2944718/71‑1 (MQ=255)
                                                                   tCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTctgctg                          <  1:1332084/43‑1 (MQ=255)
                                                                      |                                                               
TCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCCGGCATGGTGGAATGATTTTATCGCGGGTCTGGTTCTGCTGGCGGTGCTGGTGTTTGATGGACGC  >  minE/971978‑972111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑