Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,672 C→T 16.7% N41N (AAC→AAT yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,6720CT16.7% 148.8 / 15.6 89N41N (AAC→AATyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (18/56);  new base T (7/8);  total (25/64)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

CGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGC  >  minE/1781612‑1781738
                                                            |                                                                  
cgTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCaa                                                            <  1:2739296/69‑1 (MQ=255)
   cGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:436996/69‑1 (MQ=255)
   cGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:2892819/69‑1 (MQ=255)
   cGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:2923340/69‑1 (MQ=255)
   cGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:1937653/69‑1 (MQ=255)
   cGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:1781817/69‑1 (MQ=255)
   cGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:1487119/69‑1 (MQ=255)
   cGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:1191337/69‑1 (MQ=255)
    gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                        <  1:1403515/69‑1 (MQ=255)
    gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                        <  1:1298126/69‑1 (MQ=255)
    gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                        <  1:1446769/69‑1 (MQ=255)
    gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                        <  1:328320/69‑1 (MQ=255)
    gAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATAtt                                                        <  1:2413440/69‑1 (MQ=255)
     aataaTTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAatat                                                         <  1:2364394/64‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg                                                                 >  1:2543062/1‑44 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg                                                                 >  1:1964919/1‑44 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg                                                                 >  1:1996059/1‑44 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg                                                                 >  1:418625/1‑44 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg                                                                 >  1:33349/1‑44 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTg                                                                 >  1:771960/1‑44 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:434376/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:229876/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:2284979/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:2705152/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:2779239/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:3065988/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:499907/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:12464/47‑1 (MQ=255)
                    aGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTc                                                              <  1:1866111/47‑1 (MQ=255)
                        aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:1926195/69‑1 (MQ=255)
                        aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:1447130/69‑1 (MQ=255)
                        aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:2527170/69‑1 (MQ=255)
                        aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:2515442/69‑1 (MQ=255)
                        aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:3330570/69‑1 (MQ=255)
                        aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:73585/69‑1 (MQ=255)
                        aGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:1233766/69‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:62715/68‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:790044/68‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:3415152/68‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCAc                                    <  1:1915364/68‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                   <  1:141617/69‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                   <  1:720225/69‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                   <  1:1568876/69‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                   <  1:1150781/69‑1 (MQ=255)
                         gACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACt                                   <  1:1425186/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGTCACTAttgtt                             <  1:2981319/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                             <  1:2844175/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                             <  1:1426249/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                             <  1:142868/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                             <  1:3250214/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                             <  1:2577819/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                             <  1:1908303/69‑1 (MQ=255)
                               cTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTAttgtt                             <  1:795239/69‑1 (MQ=255)
                                   tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGc                                      >  1:2427475/1‑56 (MQ=37)
                                   tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCggg                          >  1:2574209/1‑68 (MQ=255)
                                   tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                         >  1:1077078/1‑69 (MQ=255)
                                   tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                         >  1:562729/1‑69 (MQ=255)
                                   tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                         >  1:358299/1‑69 (MQ=255)
                                   tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                         >  1:1910621/1‑69 (MQ=255)
                                   tAAGCCTGGTGTTAATCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGa                         >  1:2081009/1‑69 (MQ=255)
                                    aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                        >  1:206445/1‑69 (MQ=255)
                                    aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                        >  1:1487992/1‑69 (MQ=255)
                                    aaGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGAt                        >  1:2650515/1‑69 (MQ=255)
                                              ttAATCGGCAATAACCTGGTCAATATTCTTGCCTc                                                >  1:615845/1‑35 (MQ=255)
                                              ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc              >  1:1378778/1‑69 (MQ=255)
                                              ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc              >  1:2387500/1‑69 (MQ=255)
                                                  tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:194329/68‑1 (MQ=255)
                                                  tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:1895089/68‑1 (MQ=255)
                                                  tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:635200/68‑1 (MQ=255)
                                                  tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:675052/68‑1 (MQ=255)
                                                  tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:2835426/68‑1 (MQ=255)
                                                  tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:3033260/68‑1 (MQ=255)
                                                  tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:964196/68‑1 (MQ=255)
                                                   cGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc           <  1:3029205/67‑1 (MQ=255)
                                                    ggCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGAt        >  1:2041821/1‑69 (MQ=255)
                                                          aaCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGc  >  1:1583131/1‑69 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:2883012/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   >  1:505960/1‑67 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:552721/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:419064/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:2233189/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:187804/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:2505300/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   >  1:2692847/1‑67 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:3083514/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:2846702/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   >  1:12834/1‑67 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   <  1:1066511/67‑1 (MQ=255)
                                                           aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg   >  1:2963043/1‑67 (MQ=255)
                                                            |                                                                  
CGTCGAAAAATTGCTGCGTAAGCCAGACCGCCTGATAAGCCTGGTGTTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGC  >  minE/1781612‑1781738

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: