Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,633,534 T→C 45.6% G512G (GGT→GGC narQ → sensory histidine kinase in two‑component regulatory system with NarP (NarL)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,633,5340TC45.6% 11.0 / 75.7 90G512G (GGT→GGCnarQsensory histidine kinase in two‑component regulatory system with NarP (NarL)
Reads supporting (aligned to +/- strand):  ref base T (33/16);  new base C (23/18);  total (56/34)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGC  >  minE/1633479‑1633583
                                                       |                                                 
gCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAAcc                                        >  1:3435933/1‑67 (MQ=255)
gCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAAcc                                        >  1:3323721/1‑67 (MQ=255)
    cACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                  <  1:326766/69‑1 (MQ=255)
    cACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                  <  1:2516715/69‑1 (MQ=255)
    cACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                  <  1:1438784/69‑1 (MQ=255)
    cACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                  <  1:2520848/69‑1 (MQ=255)
     aCCGCGCCGTACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:3241011/1‑69 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGaaa                                         <  1:556631/61‑2 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:2568942/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:2644226/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:2974562/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:3033088/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:353217/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:37100/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:2080436/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:2058506/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:2023913/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:2001778/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:1963468/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:1013953/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:752040/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:892916/61‑1 (MQ=255)
     aCCGCGCCGGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:1462525/61‑1 (MQ=255)
     aCCGCGCCGGATGACAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:166545/61‑1 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:2016668/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:2649093/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:1666209/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:1650083/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:3126119/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:944571/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:3093622/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:2487172/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGaa                                 >  1:2913046/1‑69 (MQ=255)
     aCCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                  <  1:97783/68‑1 (MQ=255)
     aCCGCGCCGATGGCAATCACACAGTCTATATTCGCGACAACGGGATTGGCATCGGTGAAc                                         <  1:907889/60‑1 (MQ=255)
        tcgcCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAcc                               >  1:2831842/2‑68 (MQ=255)
        gcgcCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAcc                               >  1:768402/1‑68 (MQ=255)
        gcgcCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAcc                               >  1:3408336/1‑68 (MQ=255)
        gcgcCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAcc                               >  1:270501/1‑68 (MQ=255)
            cGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAg                          <  1:348808/69‑1 (MQ=255)
                        acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt               >  1:2768475/1‑68 (MQ=255)
                        acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt               >  1:1406361/1‑68 (MQ=255)
                        acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt               >  1:3014736/1‑68 (MQ=255)
                        acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt               >  1:645201/1‑68 (MQ=255)
                               ttATATTCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:441647/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGaaa                                    >  1:2301360/1‑40 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:1417266/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:3433130/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:2274134/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:1157487/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:742873/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:68106/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:3150190/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:3139942/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:1816704/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAg                                   >  1:195205/1‑41 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat       <  1:1063339/69‑1 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat       <  1:252171/69‑1 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat       <  1:1837728/69‑1 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat       <  1:1787130/69‑1 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat       <  1:1306956/69‑1 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat       <  1:1088279/69‑1 (MQ=255)
                               ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat       <  1:3069985/69‑1 (MQ=255)
                                tataTCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTg            <  1:914314/63‑1 (MQ=255)
                                tataTCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTg            <  1:239041/63‑1 (MQ=255)
                                    tCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgc  >  1:1239349/1‑69 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:946550/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:1173032/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:90068/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:994286/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:159683/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:767058/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:1621864/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:3095366/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:1808971/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:1857189/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:492753/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:404004/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:2165866/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:2253647/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:2347877/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:246209/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:335600/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:2985394/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:3006792/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:3244194/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:3016622/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata        >  1:3034686/1‑53 (MQ=255)
                                              cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGGCTGAat         >  1:1904761/1‑52 (MQ=255)
                                                 gATTGGTATAGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTg            <  1:3358390/46‑1 (MQ=255)
                                                       |                                                 
GCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGC  >  minE/1633479‑1633583

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: