Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,455,945 T→C 10.4% intergenic (‑301/+380) yfcV ← / ← sixA putative fimbrial‑like adhesin protein/phosphohistidine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,455,9450TC10.4% 98.4 / 3.0 40intergenic (‑301/+380)yfcV/sixAputative fimbrial‑like adhesin protein/phosphohistidine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (20/15);  new base C (2/2);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TGATGTTTTCACTCCCCTGATTCAATTTAATTTTACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTT  >  NC_000913/2455797‑2456088
                                                                                                                                                    |                                                                                                                                               
tGATGTTTTCACTCCCCTGATTCAATTTAATTTTACTAATGATTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATAGCAACaat                                                                                                                                                 >  2:597164/1‑149 (MQ=255)
         ccacccccTGATTCAATTTAATTTTACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCaa                                                                                                                                        <  2:460356/145‑1 (MQ=255)
          aaacacatGATGCAATTTAATTTTACTAATGAGTAAATGTAGATTTAAATAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATAAAAAAAAAAAAAAAACAGGAAGCATATCAACAATGAATGCCAAt                                                                                                                                       <  2:473018/142‑1 (MQ=255)
                cTGATTCAATTTAATTTTACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAGAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTa                                                                                                                                 <  1:443710/149‑1 (MQ=255)
                cTGATTCAATTTAATTTTAATAAGGAGTAAAGGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTa                                                                                                                                 <  2:133133/149‑1 (MQ=255)
                 tGATTCAATTTAATTTTACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGACAATAGCTTaa                                                                                                                                >  1:102915/1‑149 (MQ=255)
                    ttCAATTTAATTTTACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTaaaaa                                                                                                                             >  2:799387/1‑149 (MQ=255)
                       aatttaatttTACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAaaaca                                                                                                                          >  1:388137/1‑149 (MQ=255)
                       aatttaatttTACTAATGAGTAAAAGTAGATTTAATTAATATATTAAAAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAagaca                                                                                                                          >  1:57818/1‑149 (MQ=255)
                       aatttaatttTACCAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATAAAAAAATAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAaaaca                                                                                                                          <  2:614674/149‑1 (MQ=255)
                              ttttACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATg                                                                                                                   >  1:623704/1‑149 (MQ=255)
                              ttttACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATg                                                                                                                   >  2:514883/1‑149 (MQ=255)
                              ttttACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAAGAAGCATATCAACAATGAAAGACAAAAGATTAAAAAACAAACAAta                                                                                                                   >  2:763822/1‑148 (MQ=255)
                                  aCTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACa                                                                                                               <  1:799387/149‑1 (MQ=255)
                                           tAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATGATGCg                                                                                                      >  1:681231/1‑149 (MQ=255)
                                           tAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAACGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCg                                                                                                      >  1:106060/1‑149 (MQ=255)
                                           tAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAACGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCg                                                                                                      >  1:106052/1‑149 (MQ=255)
                                                          ttaatATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCg                                                                                       <  2:152328/149‑1 (MQ=255)
                                                             atatatTGATAGGGGAAAGATTATCTTAAGGATCTACTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCtt                                                                                    >  1:389171/1‑149 (MQ=255)
                                                               atatTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCtttt                                                                                  >  1:323350/1‑149 (MQ=255)
                                                               atatTGATAGGGGAAAGATTATCTTAAGGATATTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATACAGCGATTGGCATTAAACCGCtttt                                                                                  >  2:9136/1‑149 (MQ=255)
                                                                tatTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGt                                                                                <  1:365985/149‑1 (MQ=255)
                                                                  ttGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCATAAAAAACAAACAATGAACATATAATGCGACTAGCATTAACCCGCTTCTGTc                                                                               >  2:148469/1‑149 (MQ=255)
                                                                     aaaGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGGTTTTGTCTat                                                                            <  2:210311/147‑1 (MQ=255)
                                                                         gggAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGatataa                                                                        >  1:699619/1‑149 (MQ=255)
                                                                         gggAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGatataa                                                                        >  1:643053/1‑149 (MQ=255)
                                                                          ggAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGatataat                                                                       >  1:245881/1‑149 (MQ=255)
                                                                              aGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGAATGGCATTAACCAGCTTTTGTCGATATAatatat                                                                   <  1:366898/149‑1 (MQ=255)
                                                                               gATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACACGCTTTTGTCGATATAAtatata                                                                  >  1:29836/1‑149 (MQ=255)
                                                                                          ggATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATCGc                                                       >  1:292033/1‑149 (MQ=255)
                                                                                                 ctAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAAtt                                                <  1:549008/149‑1 (MQ=255)
                                                                                                 ctAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAAtt                                                <  1:401520/149‑1 (MQ=255)
                                                                                                   aaaaaaaCACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTa                                              >  1:707625/1‑149 (MQ=255)
                                                                                                                 tatatatTTAGAAAAAACAGGAATCATATCAACAGTGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATCGAATTGCATTATAGCCAGAATTTAATCAAGTTTCAAtt                                <  1:793741/149‑1 (MQ=255)
                                                                                                                 tatatacATAGAAAAAACAGAAAGCAAATCAACAAAGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAAtt                                <  1:707145/149‑1 (MQ=255)
                                                                                                                               aaaCAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCGAGAATTTAATCAAGTTTCAATTACGTcctccctgct                  >  2:259734/1‑147 (MQ=255)
                                                                                                                                aaCAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATCAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAATTACGTCCTGCCTGatt                 <  2:509137/149‑1 (MQ=255)
                                                                                                                                             tCAACAACGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTAtt    <  2:106060/149‑1 (MQ=255)
                                                                                                                                             tCAACAACGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTAtt    <  2:106052/149‑1 (MQ=255)
                                                                                                                                             acaacaaTGAATGCCAATAGCTTAAAAAACAAACCATGCACATGTAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTAtt    <  2:681231/148‑1 (MQ=255)
                                                                                                                                               aacaaTGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAATTACGTCGTGCCTGATTATTCAACTCTAtttt  >  1:243958/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGATGTTTTCACTCCCCTGATTCAATTTAATTTTACTAATGAGTAAATGTAGATTTAATTAATATATTGATAGGGGAAAGATTATCTTAAGGATCTTCTAAAAAAACACGAAATATATATTTAGAAAAAACAGGAAGCATATCAACAATGAATGCCAATAGCTTAAAAAACAAACAATGAACATATAATGCGATTGGCATTAACCCGCTTTTGTCGATATAATATATAAACACTATAGCCAGAATTTAATCAAGTTTCAATTACGTCCTGCCTGATTATTCAACTCTATTTT  >  NC_000913/2455797‑2456088

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