Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,966,986 G→C intergenic (‑38/+107) tar ← / ← cheW methyl‑accepting chemotaxis protein II/chemotaxis protein CheW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,966,9860GC97.6% 275.3 / NA 83intergenic (‑38/+107)tar/cheWmethyl‑accepting chemotaxis protein II/chemotaxis protein CheW
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base C (49/31);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAT  >  NZ_CP009273/1966899‑1967069
                                                                                       |                                                                                     
aCCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCa                                                                                     <  2:2457775/90‑1 (MQ=255)
  ccccAGCACCATTACCAACAGGGGGACTAGCCGGTTACGGTTAATCATAAGGCCCCTTCCGGAAAAAAAGTTGTTCTCGTTATCGCCagg                                                                                   >  2:2887668/1‑88 (MQ=255)
  ccccAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAag                                                                                   >  2:1464276/1‑90 (MQ=255)
  ccccAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAag                                                                                   >  1:200554/1‑90 (MQ=255)
    ccAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAagga                                                                                 >  1:2409580/1‑90 (MQ=255)
     cAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAaggag                                                                                <  1:2151549/90‑1 (MQ=255)
       gCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAggg                                                                              >  2:2218954/1‑90 (MQ=255)
          ccATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAgggggg                                                                           >  1:2960530/1‑90 (MQ=255)
           cATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAaggag                                                                                <  2:2690159/84‑1 (MQ=255)
           cATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAaggag                                                                                >  1:2690159/1‑84 (MQ=255)
           cATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGa                                                                          >  1:1242831/1‑90 (MQ=255)
           cATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGa                                                                          >  2:2085447/1‑90 (MQ=255)
           cATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGa                                                                          >  2:2542283/1‑90 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  1:2575918/90‑1 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  1:354042/90‑1 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         >  1:3018020/1‑90 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  1:479521/90‑1 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  2:2576254/90‑1 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  2:69535/90‑1 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  1:1935304/90‑1 (MQ=255)
            aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         >  2:61936/1‑90 (MQ=255)
                 cAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACttt                                                                   >  2:259806/1‑90 (MQ=255)
                 cAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACttt                                                                   >  2:544247/1‑90 (MQ=255)
                  aaCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTa                                                                  >  2:143450/1‑90 (MQ=255)
                  aaCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTa                                                                  >  1:1150722/1‑90 (MQ=255)
                   aCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTAt                                                                 <  1:2733839/90‑1 (MQ=255)
                   aCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTAt                                                                 <  2:1072990/90‑1 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAaggag                                                                                <  2:1923096/75‑1 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAaggag                                                                                >  1:1923096/1‑75 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         >  2:1818924/1‑82 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         >  2:1640001/1‑82 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  1:1818924/82‑1 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         >  1:571822/1‑82 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  1:1640001/82‑1 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGaa                                                                         <  2:571822/82‑1 (MQ=255)
                    cAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTAtt                                                                >  1:2300130/1‑90 (MQ=255)
                      gCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGc                                                             >  2:2057987/1‑90 (MQ=255)
                      gCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGc                                                             >  2:524939/1‑90 (MQ=255)
                       cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCt                                                            <  2:1749456/90‑1 (MQ=255)
                       cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCt                                                            <  2:1186487/90‑1 (MQ=255)
                       cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCt                                                            <  1:982045/90‑1 (MQ=255)
                         tGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGa                                                          >  1:46068/1‑90 (MQ=255)
                           aCTACGCGGATACGGTTAATCAAAAGCCCCCTTCCTAAAAAAAAGTTGATCTCGTTATCGCAAAGGAGGGGGGAAACTTTATTGCTGATg                                                        <  1:2160583/90‑1 (MQ=255)
                              aCGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCa                                                     >  1:1058073/1‑90 (MQ=255)
                              aCGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGAAACTTTATTGCTGATGCCAc                                                    >  1:1838678/1‑90 (MQ=255)
                                gcgGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCAcc                                                   >  2:634770/1‑90 (MQ=255)
                                gcgGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCAcc                                                   >  1:765725/1‑90 (MQ=255)
                                gcgGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCAcc                                                   >  2:1479955/1‑90 (MQ=255)
                                    aTACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACccgcc                                               >  1:2429644/1‑90 (MQ=255)
                                         gTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGaa                                          >  1:2917586/1‑90 (MQ=255)
                                            aaTCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAAtt                                       >  2:2002404/1‑90 (MQ=255)
                                            aaTCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAAtt                                       >  2:2495827/1‑90 (MQ=255)
                                                 tAAGGCACCTTCCTGAAAACTAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAAt                                  >  2:309368/1‑90 (MQ=255)
                                                 tAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAAt                                  >  2:393355/1‑90 (MQ=255)
                                                 tAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAAt                                  >  1:263002/1‑90 (MQ=255)
                                                   aGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATaa                                >  1:562486/1‑90 (MQ=255)
                                                   aGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATaa                                >  1:2844221/1‑90 (MQ=255)
                                                    ggCACCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATaaa                               <  2:2300130/90‑1 (MQ=255)
                                                       aCCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATaaaaaa                            >  2:712622/1‑90 (MQ=255)
                                                       aCCTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATaaaaaa                            >  1:2677190/1‑90 (MQ=255)
                                                        ccTTCCTGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAAc                           >  2:732451/1‑90 (MQ=255)
                                                          ttCCTAAAAACAAGTTGATCTCTTTATCGCAAAGGGGGGGGAAACTTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAAccc                         <  2:2967371/90‑1 (MQ=255)
                                                          ttCCTAAAAACAAGTTGATCTCGTTACCGCAAAGGGGGGGGAAAACTTTTTTGCTGATCCCCCCCGCCGCGAAATTGAAATAAAAAAccc                         <  1:61936/90‑1 (MQ=255)
                                                            ccTAAAAACAAGTTGATCTCTTTACCGCCAGGGGGGGGGGAAATTTTTTTGCTGATGCCCCCCGCCGCGAAATTGAAATAAAAAACCCGa                       <  2:2773018/90‑1 (MQ=255)
                                                            ccTAAAAAAAAGTTGCTCTCTTTACCGCGAGGGGGGGGGAAAACTTTATTCCTGATGCCCCCCGCCGCGAAATTGAAATAAAAAACCCGa                       <  1:2896393/90‑1 (MQ=255)
                                                              taaaaaCAAGTTAATCTCGTTACCGCCAAGGAGGGGGAAAACTTTATTGCTGATCCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATg                     <  2:2429644/88‑1 (MQ=255)
                                                               gAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATgc                    >  1:1733508/1‑90 (MQ=255)
                                                               gAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATgc                    >  1:1769325/1‑90 (MQ=255)
                                                               gAAAAAAAGTTGATCCCTTACGCGCCAAGGGGGGGGAAAATTTTATTGCGGAACCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATgc                    <  1:1464276/90‑1 (MQ=255)
                                                                 aaaCAAGTTCATCCCGTTACCGCAAAGGGGGGGGAAATTTTTTTGGCTGATCCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATgcgc                  <  1:1741409/90‑1 (MQ=255)
                                                                  aaCAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCa                 >  2:1400567/1‑90 (MQ=255)
                                                                  aaCAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCa                 >  1:1820343/1‑90 (MQ=255)
                                                                  aaCAAGTTGATCTCGTTATCGCAAAGGAGGGGGAAAATTTTATTGCTGATCCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCa                 <  1:524939/90‑1 (MQ=255)
                                                                  aaCAAGTGGATTTCGTTACCGCCAGGGAGGGGGAAAATTTTATTGCGGATCCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCa                 <  1:1659610/90‑1 (MQ=255)
                                                                   aaaaGTTAATCTCGTTATCGCCAGGGGGGGGGAAAATTTTTTTGCTGATCCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAg                <  1:934514/88‑1 (MQ=255)
                                                                           aTCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCgg        >  1:2812486/1‑90 (MQ=255)
                                                                             ctctTTATCGCCAGGGGGGGGGGAATTTTTATTGCTGATGCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGt      <  2:1500411/90‑1 (MQ=255)
                                                                             ctcGTTATCGCCAGGGGGGGGGAAATCTTTTTTGCTGATGCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGt      <  2:2443735/90‑1 (MQ=255)
                                                                             ctcGTTATCGCAAAGGGGGGGGGAAACTTTATTGCTGATGCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGt      <  2:434376/90‑1 (MQ=255)
                                                                             ctcGTTATCGAAAGGGAGGGGGAAATCTTTTTTGCTGATCCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGt      <  1:2972672/90‑1 (MQ=255)
                                                                              tcGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGtt     >  2:1817459/1‑90 (MQ=255)
                                                                                 ttATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAt  >  2:47377/1‑90 (MQ=255)
                                                                                 ttATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAt  >  1:1913245/1‑90 (MQ=255)
                                                                                       |                                                                                     
ACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAT  >  NZ_CP009273/1966899‑1967069

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: