Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,592,188 C→T H178Y (CAT→TAT)  tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,1880CT86.7% 62.7 / 3.5 30H178Y (CAT→TAT) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (12/14);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

AACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGC  >  NC_000913/4591952‑4592425
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
accAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTTGCCGCGGGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTc                                                                                                                                                                                                                                  <  1:147019/249‑1 (MQ=255)
             aaaCTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATg                                                                                                                                                                                                                     >  1:102649/1‑251 (MQ=255)
             aaaCTGGGCAGATTACGAAGCGTTGCCGCGTAACCCGCGACAGAGCACCGCCGCAAAGGCAGAGATCAACCGAAATTACGATATTTATCACAATTCGCTGGCGGAGCTGATCCAACTATTTGGTGCAGGCAAAATCAACGCGTTATTTGATAAGCCGACCCAGGGATATCCGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGGACGGATCTATAATATCGCCGTCAGCGATAACAATg                                                                                                                                                                                                                     >  2:328884/1‑251 (MQ=255)
                      agATTACGAAGGTTTCCCGCGTGTCCCGTGTCAGAGATCCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGTGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACTATCGGCTCTATGATATCGCCGTCAGCGAGAACAATGCCTCCTACa                                                                                                                                                                                                            <  1:341447/250‑1 (MQ=255)
                                                                          gagaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGGCGTACTGGCGGTCAt                                                                                                                                                        >  2:75180/1‑251 (MQ=255)
                                                                           agaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTTTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATTAGGACTGTTTCGAGAAGCAGTATGTTGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACATCCAGGCGATTTGGAGTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATc                                                                                                                                                       <  1:113573/251‑1 (MQ=255)
                                                                            gaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCctct                                                                                                                                                      >  1:173040/1‑251 (MQ=255)
                                                                              tCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCAATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTAATCttt                                                                                                                                                    >  1:11613/1‑250 (MQ=255)
                                                                                                  tATCACAATTCGCGGGCGGATCTGATCCAACTATTAGGAGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCTAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTa                                                                                                                                 <  2:173040/250‑1 (MQ=255)
                                                                                                   agcacaATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaaa                                                                                                                               <  1:174174/249‑1 (MQ=255)
                                                                                                            cgcTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTCCCCGTTTGGGTTGGTGTTAAAGACTCGcct                                                                                                                      >  1:357417/1‑249 (MQ=255)
                                                                                                                           gccAACTATTAGGGGCTGGTAAAATCAACGAGTTCTTTGAGCCGCCGACCCAGGGATATCAGGACGGTTTCGATAAGCAGTATTTCGCTTACAGGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTTGTGGTCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAAt                                                                                                       <  1:217961/250‑1 (MQ=255)
                                                                                                                              aaCTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGCTCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTGTCCGTCTGGTTCGGGCTTTAAGCCTAGCTGGGAGCGCCAATGAATCGc                                                                                                    >  1:33635/1‑251 (MQ=255)
                                                                                                                                          cAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATTAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGATCAAAACGAGCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCTTGATGATCGGCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGAATGACAga                                                                                        <  2:76699/251‑2 (MQ=255)
                                                                                                                                          cAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGc                                                                                        <  1:240328/251‑1 (MQ=255)
                                                                                                                                               aaaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCg                                                                                   <  1:37362/251‑1 (MQ=255)
                                                                                                                                                  aTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATc                                                                                                      <  2:91065/229‑1 (MQ=255)
                                                                                                                                                  aTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCCCCCGCTGGTTCGGGATTAAAACCTCGCTGGGAGGGCCAATGaaacagaacgggaaagcacaactctg                                                                                >  1:91065/1‑227 (MQ=255)
                                                                                                                                                         aGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGACTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCATACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGCCTGGGTCGGGATTAAAGCGCCGCTGGTAGCGCCACAGAAACGCCTGGTTggccgccttagtcacctttc                                                                          >  1:379105/1‑231 (MQ=255)
                                                                                                                                                                     aGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTg                                                             <  1:49360/251‑1 (MQ=255)
                                                                                                                                                                          aCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATTACAATGCCTCCGACAGCCAGGCGATGTGGATTCTGGTGGGCGGGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGGAGCGCCAATGAAGCGCCTTATTGACAGCGTTCGTCATATTGCGGGCGTCGGTCTGGTGaa                                                       <  2:379105/251‑1 (MQ=255)
                                                                                                                                                                                 gATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATTGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTTACAGCATTCGTCATATTGCAGGCGGCGATCTGTTGAAACCGAtt                                                <  1:342409/251‑1 (MQ=255)
                                                                                                                                                                                    atCAGGCCGGTTGCGAGATGCAGTCTGTGGCTTACATGGAGCACAACGATCGGCGCTCTGATATCGCCGTCAGCGATAACAATGACTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAg                                             <  1:62878/251‑1 (MQ=255)
                                                                                                                                                                                         gACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTTATCTTCGCCGGCTGGTTCGGTATTAAAGCCCCGCTGGTAGCGCCAAAGAATCGCCTGATTGACACCAATCGTCATATTTCAGGCGGGGAACTGGGTGAAAccattgtgggggg                                         >  1:269130/1‑240 (MQ=255)
                                                                                                                                                                                            ggTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTc                                                                                    <  1:27078/205‑1 (MQ=255)
                                                                                                                                                                                                           tATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTCGGGATTAAAACCTCGCTGGTAGCGCCCATTAAACGCCCGATTGACAGCATTCGTCCTATTTCAGGGGGCGATTTTGGGAAACCCATTTGGGGGGAGGGGTCTAATTAGATggg                      >  1:119812/1‑251 (MQ=255)
                                                                                                                                                                                                                  cTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTTGTTCGGTATTTAAACCCCGCTGGTAGCGCCACTGAAACGCCCGATTTACAACCTTCGGCATATTGCAGGCGGCGCTATGGTGGAACCGCTTGAGGGGGAGGGGGCTAATGAGATGGGGCaca                 >  1:214246/1‑247 (MQ=255)
                                                                                                                                                                                                                    tACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGTGGCAACTgg              <  1:75180/250‑1 (MQ=255)
                                                                                                                                                                                                                               aaaCGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGTTCGGGATTAAAAACTCGCTGGGAGCGCCAAAGAAACGCCTGATTGACAGGATTCGTCATATTGTAAGGGGGGATCTGGTGAAAACGATTTAGGTGGATGGCTCTAAAGAGGTGGGGGCACTGGCAGAGAGTTTgt  >  1:291451/1‑250 (MQ=255)
                                                                                                                                                                                                                                aaCGACCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGaa                                                                                                        <  2:369227/149‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              
AACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGC  >  NC_000913/4591952‑4592425

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