Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,560,500 G→A 39.0% R169R (CGG→CGA evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,560,5000GA39.0% 42.1 / 52.9 82R169R (CGG→CGAevgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base A (11/21);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

AACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATC  >  minE/1560440‑1560560
                                                            |                                                            
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCaa                                                       <  1:982478/68‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCaa                                                       <  1:1871265/68‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCaa                                                       <  1:2615372/68‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCaa                                                       <  1:2868754/68‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCaa                                                       <  1:1303441/68‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:2044955/67‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:2405723/67‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:1685776/67‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:1570436/67‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:3095368/67‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:1367591/67‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:1088256/67‑1 (MQ=255)
aaCCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:96176/67‑1 (MQ=255)
 accaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:3108262/66‑1 (MQ=255)
 accaccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAAt                                                     >  1:2741257/1‑69 (MQ=255)
   caccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:523041/64‑1 (MQ=255)
   caccCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCa                                                        <  1:2165035/64‑1 (MQ=255)
    accCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAc                                                  >  1:2014905/1‑69 (MQ=255)
       cTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccc                                                >  1:648517/1‑68 (MQ=255)
       cTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccc                                                >  1:649373/1‑68 (MQ=255)
       cTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccc                                                >  1:720578/1‑68 (MQ=255)
       cTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccc                                                >  1:724204/1‑68 (MQ=255)
       cTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcc                                                 >  1:2740070/1‑67 (MQ=255)
       cTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGATCGGGTAGCAAATTAccc                                                >  1:358308/1‑68 (MQ=255)
          cACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCa                                            >  1:3052156/1‑69 (MQ=255)
          cACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCa                                            >  1:1070967/1‑69 (MQ=255)
          cACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCa                                            >  1:1270910/1‑69 (MQ=255)
                              tACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCatcaatc                        <  1:1097884/69‑1 (MQ=255)
                              tACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCatcaatc                        <  1:2733984/69‑1 (MQ=255)
                              tACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCatcaatc                        <  1:399491/69‑1 (MQ=255)
                               aCCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCatcaatca                       >  1:2299155/1‑69 (MQ=255)
                                ccTCATCTAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:459659/37‑1 (MQ=25)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:2848389/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:2065864/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:2977523/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:3074132/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:3081428/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:3183102/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:386051/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:1245629/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:1212875/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:666790/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:798337/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:956730/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:1073757/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:1039308/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:2046834/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:2092973/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:2145497/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:2521488/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:261142/37‑1 (MQ=38)
                                ccTCATCAAAACCAGTAAATATTGCTCGAGTAGCaaa                                                      <  1:1702333/37‑1 (MQ=38)
                                 cTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCAt                      <  1:3024513/68‑1 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCGCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:1405377/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:347378/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:1122477/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:682155/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:390309/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:351294/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:2736079/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:1544434/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:1525446/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa               >  1:299348/1‑68 (MQ=255)
                                        aaaCCAGTAAATATTGCTCGAGTAGCAAATTATCCCCCACACGAGGTAATTCATCAATCATTTCCaaa               >  1:1769137/1‑68 (MQ=255)
                                                aaaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAt      <  1:3272372/69‑1 (MQ=255)
                                                aaaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAt      <  1:3337918/69‑1 (MQ=255)
                                                aaaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAt      <  1:2844163/69‑1 (MQ=255)
                                                aaaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAt      <  1:1541983/69‑1 (MQ=255)
                                                aaaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAt      <  1:2149531/69‑1 (MQ=255)
                                                 aaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTTCAAAAGcaacaa       >  1:2190049/1‑67 (MQ=255)
                                                 aaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTTCAAAAGCAACAAt      >  1:2926413/1‑68 (MQ=255)
                                                 aaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAtt     >  1:1292572/1‑69 (MQ=255)
                                                 aaTATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAtt     >  1:2834593/1‑69 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTa    >  1:2032076/1‑68 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTa    >  1:2296698/1‑68 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTa    >  1:1921206/1‑68 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt   >  1:2931176/1‑69 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt   >  1:2440834/1‑69 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt   >  1:301037/1‑69 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt   >  1:305253/1‑69 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAt   >  1:1738539/1‑69 (MQ=255)
                                                   tatTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCACAATTAtc  >  1:3030000/1‑69 (MQ=255)
                                                            |                                                            
AACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATC  >  minE/1560440‑1560560

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: