Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,972,741 (G)6→7 intergenic (‑50/+95) tar ← / ← cheW methyl‑accepting chemotaxis protein Tar/chemotaxis protein CheW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,972,7351.G87.1% 92.6 / 7.4 31intergenic (‑44/+101)tar/cheWmethyl‑accepting chemotaxis protein Tar/chemotaxis protein CheW
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (22/5);  total (24/7)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AAGAAAAAAACAGACTGCCGGAAATAAGCTGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAG  >  NC_000913/1972598‑1972868
                                                                                                                                            |                                                                                                                                        
aagaaAAAAACAGACTGCCGGAAATAAGCTGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAaggagg                                                                                                                                        >  1:257604/1‑140 (MQ=255)
  gaaAAAAACAGACTGCCGGAAATAAGCTGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAgggg                                                                                                                                      <  2:257604/140‑1 (MQ=255)
      aaaaCAGACTGCCGGAAATAAGCTGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAggggggga                                                                                                                                 <  1:87198/140‑3 (MQ=255)
                               tAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACccgcc                                                                                                       >  2:30830/1‑140 (MQ=255)
                                     tGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGaaa                                                                                                 >  1:236284/1‑140 (MQ=255)
                                       cGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAAtt                                                                                               >  1:99754/1‑140 (MQ=255)
                                       cGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAAtt                                                                                               >  1:710502/1‑140 (MQ=255)
                                          acagccccAGCACCATTACCAACAGCGTGACTACGGGGATACGGTTAATCAAAAGGCCCCTTCCTAAAAACAAGTTGATCTCTTTATCGGAAAGGGGGGGGGGAAACTTTATTGCTGATGCCCCCCGCCGCGAAATTGaa                                                                                            <  2:495444/136‑1 (MQ=255)
                                                    cACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAAcc                                                                                  <  1:30830/140‑1 (MQ=255)
                                                        aTTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGAt                                                                              >  1:124204/1‑140 (MQ=255)
                                                         ttACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAAcccctgt                                                                             >  1:495444/1‑136 (MQ=255)
                                                              aaCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAg                                                                        >  1:688216/1‑140 (MQ=255)
                                                               aCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGa                                                                       >  1:77022/1‑140 (MQ=255)
                                                                 aGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGTGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGatc                                                                     >  1:833470/1‑140 (MQ=255)
                                                                   cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAGAAACCCGATGCGCAGatcat                                                                   >  1:712070/1‑140 (MQ=255)
                                                                   cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGatcct                                                                   >  1:417667/1‑138 (MQ=255)
                                                                   cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGatcat                                                                   >  1:388336/1‑140 (MQ=255)
                                                                       aCTCCGCGGATACGGTTAATCATAAGGCTCCTTCTGAAAAACAAGTTGATCTCGTTATCGGCAAGGGGGGGGGGAACCTTTTTTGTTGTCGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCACATCATCggg                                                               <  2:833470/140‑1 (MQ=255)
                                                                         agcgcgGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGtt                                                             >  2:555356/3‑140 (MQ=255)
                                                                            gcgGATACGTTAAATCATAAGCCCCCTTCCTAAAAACAAGTTGATCTCGTTACCGGAAAGGGGGGGGGAAAATTTTTTTGCTGATGCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAt                                                          <  2:124204/140‑1 (MQ=255)
                                                                                aTACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCa                                                      >  2:376843/1‑140 (MQ=255)
                                                                                 tACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCaa                                                     >  1:896431/1‑140 (MQ=255)
                                                                                   cGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAAtt                                                   >  1:323964/1‑140 (MQ=255)
                                                                                     gTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGa                                                 >  1:69174/1‑140 (MQ=255)
                                                                                          tCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGaaa                                            >  2:500443/1‑140 (MQ=255)
                                                                                               agcccccTCCCTAAAAACAAGTTGATCTCTTTATGCGCAGGGGGGGGGGAAAATTTTATTGCTGAGGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCggg                                       <  2:663519/135‑1 (MQ=255)
                                                                                                   aCCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGaa                                   <  1:555356/140‑1 (MQ=255)
                                                                                                    ccTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAAt                                  >  2:930224/1‑140 (MQ=255)
                                                                                                    ccTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAAt                                  >  1:447932/1‑140 (MQ=255)
                                                                                                    ccTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAAt                                  >  2:645685/1‑140 (MQ=255)
                                                                                                    ccTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAAt                                  <  1:500443/140‑1 (MQ=255)
                                                                                                       tCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTAc                               >  1:267407/1‑140 (MQ=255)
                                                                                                            aaaaCAAGTTGATCTCGTTATCGGCAAGGAGGGGAGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCAc                          >  1:768091/1‑140 (MQ=255)
                                                                                                                   tttGATCTCGTTATGGGCAAGGAGGGGAGAAAATTTTATTGCTGATCCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAc                   <  2:768091/139‑1 (MQ=255)
                                                                                                                       ttctcTTTACCGGAAAGGGGGGGGGAAAATTTTTTTGCTGATGCCCCCCGCCGCGAAATTGAAATAAGAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAcgccg               <  2:712070/139‑1 (MQ=255)
                                                                                                                       aTCTCGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAcgccg               <  2:388336/140‑1 (MQ=255)
                                                                                                                       aTCTCGTTACCGGAAGGGGGGGGGGAAAATTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAcgccg               <  2:447932/140‑1 (MQ=255)
                                                                                                                        tctcGTTAGGGGCAGGGGGGGGGGAAATTTTTATTGTTGATCCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAcgccgc              <  2:236284/140‑1 (MQ=255)
                                                                                                                         ctcGTTATCGGCAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCcgcg             <  1:645685/140‑1 (MQ=255)
                                                                                                                         agacagaTCGGAAAGGGGGGGGGAAAATTTTATTGCTGATGCCCCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCcgcg             <  2:28896/134‑1 (MQ=255)
                                                                                                                                    cAAGGAGGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAg  >  2:64790/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                        
AAGAAAAAAACAGACTGCCGGAAATAAGCTGTAACAGTGCGAATACCCCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAG  >  NC_000913/1972598‑1972868

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: