Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,143,119 G→A 39.9% R521H (CGC→CAC)  katE → hydroperoxidase HPII(III)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,143,1190GA39.9% 54.7 / 74.8 111R521H (CGC→CAC) katEhydroperoxidase HPII(III)
Reads supporting (aligned to +/- strand):  ref base G (27/39);  new base A (19/25);  total (46/65)
Fisher's exact test for biased strand distribution p-value = 8.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

CCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAG  >  minE/1143052‑1143180
                                                                   |                                                             
ccATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAgcgc                                                              <  1:2537164/69‑1 (MQ=255)
ccATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAgcgc                                                              <  1:2341061/69‑1 (MQ=255)
ccATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGAAgcgc                                                              <  1:1155322/69‑1 (MQ=255)
     gTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCatat                                                         >  1:2700763/1‑69 (MQ=255)
     gTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCatat                                                         >  1:40254/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAgcg                                                               >  1:3064977/1‑53 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATgg                                                >  1:3217062/1‑68 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATgg                                                >  1:2813218/1‑68 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATgg                                                >  1:2783132/1‑68 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATgg                                                >  1:807149/1‑68 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATgg                                                >  1:2216727/1‑68 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATgg                                                >  1:2266351/1‑68 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATgg                                                >  1:163921/1‑68 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:635229/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:229569/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:1947842/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:253326/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:2565127/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:2581087/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:3298473/1‑69 (MQ=255)
               tattatTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGt                                               >  1:870498/1‑69 (MQ=255)
                   attCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCa                                           >  1:2950620/1‑69 (MQ=255)
                   attCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCa                                           >  1:740018/1‑69 (MQ=255)
                   attCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCa                                           >  1:2389738/1‑69 (MQ=255)
                     tCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:2951604/68‑1 (MQ=255)
                     tCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:1284745/68‑1 (MQ=255)
                     tCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:2683854/68‑1 (MQ=255)
                     tCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:2632795/68‑1 (MQ=255)
                     tCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:1370717/68‑1 (MQ=255)
                     tCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:2561276/68‑1 (MQ=255)
                     tCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:801526/68‑1 (MQ=255)
                      cccATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:2622827/67‑1 (MQ=255)
                      cccATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:2222745/67‑1 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAg                              >  1:873879/1‑67 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAg                              >  1:1174258/1‑67 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAg                              >  1:417626/1‑67 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCa                            <  1:2451664/69‑1 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCa                            <  1:431614/69‑1 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCa                            <  1:1250319/69‑1 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCa                            <  1:2902410/69‑1 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCa                            <  1:2589840/69‑1 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCa                            <  1:3416584/69‑1 (MQ=255)
                                  tGTTCTGGCTAAGTCAGACGCAATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCa                            <  1:793092/69‑1 (MQ=255)
                                      cTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAg                         <  1:311083/68‑1 (MQ=255)
                                         gCTAAGTCAGTCGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:17896/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCAGCATATTGTCGATGGTTTCAGtt                                        <  1:38409/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1894150/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:656403/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1235471/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:2885713/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:2934753/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:3157006/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1758702/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:3345524/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1929365/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1939018/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1684879/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:484206/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1424844/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:736368/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1604146/50‑1 (MQ=255)
                                         gCTAAGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:791830/50‑1 (MQ=255)
                                          cTACGTCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTGAGtt                                        <  1:381505/49‑1 (MQ=255)
                                             aGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAg                                          <  1:1901746/44‑1 (MQ=255)
                                               tCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTc                >  1:1842647/1‑68 (MQ=255)
                                               tCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTc                >  1:2970309/1‑68 (MQ=255)
                                               tCAGACGCCATTCGAGCAGCACCATATTGTCGATGGTTTCAGtt                                        <  1:1139671/44‑1 (MQ=38)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:929512/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:874494/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:1022247/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:1651229/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:3248791/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:1656830/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:337695/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:1693268/68‑1 (MQ=255)
                                                   aCGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGta            <  1:2430493/68‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:3011654/69‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:658373/69‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:535735/69‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:426144/69‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:2941427/69‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:2780022/69‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:2855660/69‑1 (MQ=255)
                                                    cGCCATTTGAGCAGCCCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGtata          <  1:2466645/69‑1 (MQ=255)
                                                           cgAGCAGGACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:3082105/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:2932065/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:892593/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:3117094/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:1270479/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:1562366/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:157573/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:1642074/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:1700764/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:542446/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:2580949/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:2110534/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:2508427/2‑69 (MQ=255)
                                                           cgAGCAGCACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGa   >  1:2507319/2‑69 (MQ=255)
                                                             agcagcACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGAg  <  1:2631239/68‑1 (MQ=255)
                                                             agcagcACCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTAGTTCGTCCGTATATTCGTGAg  <  1:1359999/68‑1 (MQ=255)
                                                                   |                                                             
CCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAG  >  minE/1143052‑1143180

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: