Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,973,098 G→A 12.0% A436A (GCC→GCT tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,973,0980GA12.0% 130.7 / 7.7 50A436A (GCC→GCTtarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base G (19/25);  new base A (4/2);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

ATCTCGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATG  >  W3110S.gb/1973036‑1973152
                                                              |                                                      
aTCTCGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTc                                                <  1:479855/71‑1 (MQ=255)
aTCTCGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTc                                                <  1:1006218/71‑1 (MQ=255)
 tctcGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTc                                                >  1:2332390/1‑70 (MQ=255)
 tctcGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTc                                                >  1:1992675/1‑70 (MQ=255)
 tctcGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTc                                                >  1:3192184/1‑70 (MQ=255)
 tctcGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTc                                                >  1:2692842/1‑70 (MQ=255)
     gCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCg                                               <  1:2712765/67‑1 (MQ=255)
     gCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGAc                                             <  1:1131026/69‑1 (MQ=255)
          tAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGgcgc                                                    <  1:65025/57‑1 (MQ=255)
            aTGTCAGTCACGCTAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCAGCGCTTTAGACCAGCACCGa                                    <  1:2630569/71‑1 (MQ=38)
            aTGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCa                                        >  1:670023/1‑67 (MQ=255)
            aTGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCa                                        >  1:1715214/1‑67 (MQ=255)
            aTGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCa                                        >  1:3456527/1‑67 (MQ=255)
            aTGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCg                                     <  1:1875607/70‑1 (MQ=255)
            aTGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCAGCGCTTTCGACCAGCACCGa                                    <  1:2927454/71‑1 (MQ=255)
                 agtcaCGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCg                                <  1:1080397/70‑1 (MQ=255)
                     aCGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGAACAGCACCGAACCGGTAt                            >  1:3686000/1‑70 (MQ=255)
                        cgAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCaa                         >  1:1978770/1‑70 (MQ=255)
                              aCAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGt                    >  1:105753/1‑69 (MQ=255)
                                  cATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACg              <  1:2056653/71‑1 (MQ=255)
                                  cATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACg              <  1:1409974/71‑1 (MQ=255)
                                   aTTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGa                                              >  1:3548859/1‑38 (MQ=255)
                                    ttGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGt                    <  1:3574753/63‑1 (MQ=255)
                                    ttGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACg              >  1:1818890/1‑69 (MQ=255)
                                    ttGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGa            <  1:2378877/71‑1 (MQ=255)
                                    ttGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGa            <  1:3820402/71‑1 (MQ=255)
                                    ttGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGa            <  1:3012786/71‑1 (MQ=255)
                                      gACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTga                  <  1:2855934/63‑1 (MQ=255)
                                      gACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGt          >  1:1672640/1‑71 (MQ=255)
                                      gACGATATTGTTCATTGTTTCCCCAGCGCTTTCGa                                              >  1:370620/1‑35 (MQ=255)
                                      gACGATATTGTTCATTGTTTCCCCAGCGCTTTCGa                                              >  1:2056584/1‑35 (MQ=255)
                                      gACGATATTGTTCATTGTTTCCCCAGCGCTTTCGa                                              >  1:290090/1‑35 (MQ=255)
                                      gACGATATTGTTCATTGTTTCCCCAGCGCTTTCGa                                              >  1:3361731/1‑35 (MQ=255)
                                         gATATTGTTCATTGTTTCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACg              <  1:3184473/63‑1 (MQ=255)
                                          atatTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTc      <  1:1539787/71‑1 (MQ=255)
                                             ttGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATc                           <  1:2625920/47‑1 (MQ=255)
                                             ttGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAAt   <  1:1751005/71‑1 (MQ=255)
                                             ttGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAAt   <  1:2516448/71‑1 (MQ=255)
                                              tGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAAt   <  1:2181186/70‑1 (MQ=255)
                                               gTTCATTGTTTCCCCGGCGCTTTCGACCAGCCCCg                                     >  1:2960769/1‑35 (MQ=255)
                                               gTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAAc                                  >  1:1302716/1‑38 (MQ=255)
                                               gTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATg  <  1:956308/70‑1 (MQ=255)
                                               gTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATg  <  1:268854/70‑1 (MQ=255)
                                                ttCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTc         <  1:3579484/62‑1 (MQ=255)
                                                ttCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTc         <  1:2188657/62‑1 (MQ=255)
                                                ttCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTc         <  1:1338052/62‑1 (MQ=255)
                                                 tCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTa                             >  1:1569414/1‑41 (MQ=255)
                                                     tGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTgag                 >  1:3607683/1‑49 (MQ=255)
                                                     tGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTgag                 >  1:374465/1‑49 (MQ=255)
                                                     tGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTgag                 >  1:3542378/1‑49 (MQ=255)
                                                              |                                                      
ATCTCGCCCATAATGTCAGTCACGCGAGTGACAGCATTGACGATATTGTTCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATG  >  W3110S.gb/1973036‑1973152

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: