Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,844,797 G→A 59.8% D96D (GAC→GAT cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,844,7970GA59.8% 31.7 / 57.7 82D96D (GAC→GATcysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base G (16/17);  new base A (35/14);  total (51/31)
Fisher's exact test for biased strand distribution p-value = 4.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCC  >  minE/1844732‑1844851
                                                                 |                                                      
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2291777/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2037072/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:285476/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2033244/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2012579/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2974800/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2570858/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2480871/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:666738/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2971518/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:1180208/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:726758/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2335718/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2377717/1‑68 (MQ=255)
tCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:878152/1‑68 (MQ=255)
 cGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    >  1:2865178/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATc                                                      >  1:2308181/1‑63 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAAtag                                                 >  1:1041629/1‑66 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:2887897/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:2471966/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:27478/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:2118874/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:2169801/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:2671011/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:2302887/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:242076/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:2922765/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:1006267/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:1113933/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:851247/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:381834/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:1555895/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:1547381/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:1194362/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:1464000/1‑69 (MQ=255)
     cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                >  1:690198/1‑69 (MQ=255)
                 ggCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTc                                    >  1:165683/1‑69 (MQ=255)
                 ggCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTc                                    >  1:531122/1‑69 (MQ=255)
                       aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATTGt                                                >  1:2238042/1‑51 (MQ=255)
                       aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                >  1:3156448/1‑51 (MQ=255)
                       aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                >  1:410587/1‑51 (MQ=255)
                       aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                >  1:353037/1‑51 (MQ=255)
                       aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                >  1:3239518/1‑51 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:72675/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:3194934/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:1019482/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:1661965/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:2210678/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:1712495/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:1922536/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:2061794/44‑1 (MQ=255)
                          aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCaa                                                    <  1:1984771/44‑1 (MQ=255)
                             tttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                        <  1:2901993/69‑1 (MQ=255)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                        <  1:989999/68‑1 (MQ=255)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                        <  1:1070175/68‑1 (MQ=255)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                        <  1:1524828/68‑1 (MQ=255)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                        <  1:2109447/68‑1 (MQ=255)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:1571744/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:698926/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:1321044/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:739054/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:65939/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:2510057/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:1164686/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:3365525/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:2549087/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:1865437/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:3123852/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:3062076/48‑1 (MQ=39)
                              ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:2025240/48‑1 (MQ=39)
                              ttACGCTTCTAGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                            <  1:1560726/48‑1 (MQ=255)
                                aCGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                        <  1:398088/66‑1 (MQ=255)
                                                   tAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGcc  <  1:1221712/69‑1 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:2250083/1‑65 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:2482173/1‑65 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:950700/1‑65 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:1178773/1‑65 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:1369663/1‑65 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:2194153/1‑65 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:1874058/1‑65 (MQ=255)
                                                    aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga     >  1:1874136/1‑65 (MQ=255)
                                                      cGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                        <  1:1537350/44‑1 (MQ=255)
                                                                 |                                                      
TCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCC  >  minE/1844732‑1844851

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: