Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,969,214 G→A 8.2% S334S (AGC→AGT cheB ← fused chemotaxis regulator and protein‑glutamate methylesterase in two‑component regulatory system with CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,969,2140GA 8.2% 217.8 / 7.0 73S334S (AGC→AGTcheBfused chemotaxis regulator and protein‑glutamate methylesterase in two‑component regulatory system with CheA
Reads supporting (aligned to +/- strand):  ref base G (40/27);  new base A (3/3);  total (43/30)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACA  >  W3110S.gb/1969144‑1969284
                                                                      |                                                                      
tCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTgg                                                                        <  1:1241817/71‑1 (MQ=255)
tCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTgg                                                                        <  1:1347037/71‑1 (MQ=255)
  cGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCt                                                                      <  1:1203221/71‑1 (MQ=255)
  cGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCt                                                                      >  1:2828167/1‑71 (MQ=255)
   gCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCt                                                                      >  1:3810316/1‑70 (MQ=255)
       tttCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTAc                                                                   <  1:3410372/69‑1 (MQ=255)
       tttCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCt                                                                 <  1:354546/71‑1 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:3349455/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:3363604/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:2260081/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:4116927/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:1989009/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:1252787/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:1260862/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:584193/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTgg                                                               >  1:1391771/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTAACTgg                                                               >  1:1917115/1‑71 (MQ=255)
         tCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCAGGCTTACCTgg                                                               >  1:1387932/1‑71 (MQ=255)
           acacTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGc                                                              >  1:3621944/1‑70 (MQ=255)
           acacTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCt                                                             <  1:904098/71‑1 (MQ=255)
            cacTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTa                                                            <  1:1547633/71‑1 (MQ=255)
              cTCCTGATTTAAATACGTATCGCCTGTGCGGCACTAATTTTTGCCAACATTTGCTGGCTTAGCTGGCTAAg                                                          <  1:1767871/71‑1 (MQ=255)
              cTCCTGATTTAAATACGTATCGCCTGTGCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAg                                                          <  1:1762721/71‑1 (MQ=255)
              cTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAg                                                          <  1:607852/71‑1 (MQ=255)
              cTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAg                                                          <  1:2217153/71‑1 (MQ=255)
              cTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAg                                                          <  1:858955/71‑1 (MQ=255)
                 cTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGc                                                                       >  1:2013751/1‑55 (MQ=255)
                 cTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATc                                                       <  1:2652478/71‑1 (MQ=255)
                  tGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATc                                                       >  1:731523/1‑70 (MQ=255)
                    aTTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGa                                                     >  1:70374/1‑70 (MQ=255)
                    aTTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGAc                                                    >  1:748677/1‑71 (MQ=255)
                    aTTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGAc                                                    >  1:1404444/1‑71 (MQ=255)
                    aTTTAAATACGTATCGCCTGTCCGACACTAATTTTTGCCAACATTTGCTGGCtt                                                                     <  1:1757511/54‑1 (MQ=255)
                      ttAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCa                                                  >  1:3264436/1‑71 (MQ=255)
                       tAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCa                                                  >  1:2395862/1‑70 (MQ=255)
                       tAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCa                                                  >  1:301379/1‑70 (MQ=255)
                             cGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCa                                           <  1:1977846/71‑1 (MQ=255)
                             cGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCa                                           <  1:1013242/71‑1 (MQ=255)
                              gTATCGCCTGGCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAg                                          <  1:427468/71‑1 (MQ=255)
                                   gCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGacac                                      <  1:3159538/70‑1 (MQ=255)
                                     cTGTCCGTCACTTATTTTTGGCAACATTTGCTGGCTTACCTGGCTa                                                            <  1:223025/46‑1 (MQ=37)
                                            gCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATAtt                            <  1:193252/71‑1 (MQ=255)
                                            gCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATAtt                            <  1:1508648/71‑1 (MQ=255)
                                             cACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTg                           <  1:2864361/71‑1 (MQ=255)
                                               cTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACt                                                <  1:2326854/48‑1 (MQ=255)
                                                tAATTTTTGCCTACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAg                                          <  1:3168008/53‑1 (MQ=255)
                                                tAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCa                                           <  1:3545539/52‑1 (MQ=255)
                                                tAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATg                        >  1:4037742/1‑71 (MQ=255)
                                                 aaTTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGa                                         >  1:4256384/1‑53 (MQ=255)
                                                 aaTTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGa                                         >  1:1820258/1‑53 (MQ=255)
                                                  aTTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGc                      >  1:974217/1‑71 (MQ=255)
                                                  aTTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGc                      >  1:3860181/1‑71 (MQ=255)
                                                   tttttGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGcc                     <  1:1790508/71‑1 (MQ=255)
                                                    ttttGCCAACATTTGCTGGCTTACCTGGCTAAGATCGa                                                     >  1:2400176/1‑38 (MQ=255)
                                                     tttGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCt                    >  1:3554429/1‑70 (MQ=255)
                                                       tGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGAt                         >  1:4020352/1‑63 (MQ=255)
                                                         ccAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGTCACCACCCATAtt                            <  1:3230374/58‑1 (MQ=255)
                                                          cagcATTTGCTGACTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTcg                  <  1:2274981/64‑1 (MQ=255)
                                                          cagcATTTGCTGACTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTcg                  <  1:2024765/64‑1 (MQ=255)
                                                          cagcATTTGCTGACTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTcg                  <  1:1017841/64‑1 (MQ=255)
                                                           aaCATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCa                                           <  1:963314/41‑1 (MQ=255)
                                                           aaCATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCgg             >  1:2704712/1‑71 (MQ=255)
                                                           aaCATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCgg             >  1:2888089/1‑71 (MQ=255)
                                                           aaCATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCgg             >  1:1472111/1‑71 (MQ=255)
                                                           aaCATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCgg             <  1:1409376/71‑1 (MQ=255)
                                                            gcATTTGCTGACTTACCTGGCTAAGATCGACCACTTCGCAGAcacc                                     >  1:1757376/2‑46 (MQ=255)
                                                            gcATTTGCTGACTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGa                          >  1:2816792/2‑57 (MQ=255)
                                                            gcATTTGCTGACTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGa                          >  1:2611987/2‑57 (MQ=255)
                                                              aTTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCAt          >  1:3877915/1‑71 (MQ=255)
                                                              aTTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCAt          >  1:249660/1‑71 (MQ=255)
                                                              aTTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCAt          >  1:190485/1‑71 (MQ=255)
                                                                ttGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGa                          >  1:9968/1‑53 (MQ=255)
                                                                ttGCTGGCTTACCTGGCTAAGATCGACCACTTAGCAGCCACCACCCATATTGATGGCCTCGCGCGGCATGc        >  1:4010979/1‑71 (MQ=255)
                                                                 tGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCa                                >  1:2715238/1‑46 (MQ=255)
                                                                 tGCTGGCTTACCTGGCTAAGATCCACCACTTCGCa                                           >  1:4207017/1‑35 (MQ=255)
                                                                      gCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAaca  >  1:1856044/1‑71 (MQ=255)
                                                                      |                                                                      
TCCGCCATTTCACACTCCTGATTTAAATACGTATCGCCTGTCCGGCACTAATTTTTGCCAACATTTGCTGGCTTACCTGGCTAAGATCGACCACTTCGCAGACACCACCCATATTGATGGCCTCGCGCGGCATGCCGAACA  >  W3110S.gb/1969144‑1969284

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: