Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,902 G→C 9.3% A1129G (GCA→GGA)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,648,9020GC 9.3% 132.7 / 4.3 42A1129G (GCA→GGA) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base C (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGT  >  NC_000913/2648757‑2649047
                                                                                                                                                 |                                                                                                                                                 
gATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCg                                                                                                                                                >  5:187550/1‑149 (MQ=255)
       cGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGc                                                                                                                                         <  8:15966/149‑1 (MQ=255)
            gATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCcggac                                                                                                                                    >  6:314734/1‑149 (MQ=255)
                                      acaacaaTTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGTAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGc                                                                                                          >  7:184594/1‑149 (MQ=255)
                                      acaacaaTTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGc                                                                                                          <  1:202465/149‑1 (MQ=255)
                                       caacaaTTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGCCGGACGCCGATTTTCCACTGCTTATGCt                                                                                                         >  5:211858/1‑149 (MQ=255)
                                        aacaaTTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCcggacgga                                                                                                                                 >  8:34871/1‑124 (MQ=255)
                                        aacaaTTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCcggacgga                                                                                                                                 <  7:34871/124‑1 (MQ=255)
                                        aacaaTTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCt                                                                                                         >  3:281618/1‑148 (MQ=255)
                                           aaTTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCcggac                                                                                                                                    <  2:31922/118‑1 (MQ=255)
                                                  cTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACg                                                                                              <  7:294849/149‑1 (MQ=255)
                                                      cAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGttt                                                                                          >  5:357799/1‑149 (MQ=255)
                                                        gCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCa                                                                                        >  6:216020/1‑149 (MQ=255)
                                                            aaCAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGttt                                                                                          <  5:117240/143‑1 (MQ=255)
                                                            aaCAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGttt                                                                                          >  6:117240/1‑143 (MQ=255)
                                                                ggCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGATGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACTCCGGACGGACTCCGATTTTCCACTGCTTATGCGGATCGGCAACGGTTTCACCCGGTaa                                                                                <  6:103744/149‑1 (MQ=255)
                                                                 gCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAAc                                                                               <  6:130234/149‑1 (MQ=255)
                                                                         ttttGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAAccc                                                                       >  3:354693/1‑149 (MQ=255)
                                                                         ttttGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCAGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAAccc                                                                       >  5:327951/1‑149 (MQ=255)
                                                                              cAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTaa                                                                  <  5:314734/149‑1 (MQ=255)
                                                                                atccatccGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAAtg                                                                <  6:187550/149‑1 (MQ=255)
                                                                                 tccatccGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAAtgg                                                               >  2:54626/1‑149 (MQ=255)
                                                                                   catccGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAAtgg                                                               >  8:158758/1‑147 (MQ=255)
                                                                                   catccGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAAtgg                                                               <  7:158758/147‑1 (MQ=255)
                                                                                           ggAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAAt                                                     >  8:277377/1‑149 (MQ=255)
                                                                                                 gCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTtctc                                                >  2:50035/1‑148 (MQ=255)
                                                                                                 gCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTcc                                               <  4:354693/149‑1 (MQ=255)
                                                                                                 gCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTcc                                               <  8:184594/149‑1 (MQ=255)
                                                                                                               aGGCTGTAACGCCGTCCCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCACCATAACCCg                                 >  1:98828/1‑149 (MQ=255)
                                                                                                                ggCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAGCGCTaa                                                                          >  4:114656/1‑107 (MQ=255)
                                                                                                                ggCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAGCGCTaa                                                                          <  3:114656/107‑1 (MQ=255)
                                                                                                                ggCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCgg                                <  6:211858/149‑1 (MQ=255)
                                                                                                                        cGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTc                        <  5:216020/149‑1 (MQ=255)
                                                                                                                         gCCGTACCGTAATTAACCGTTTGTCCCGGGAACGCCGGACGGCCGCCGATTTTCCACTGCTTATGCTGCTCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCg                       >  5:229117/1‑149 (MQ=255)
                                                                                                                         gCCGTACCGTAATTAACCGTTTGTCCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCg                       >  5:234448/1‑149 (MQ=255)
                                                                                                                                        aCCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAAccc                                                                       <  8:266787/86‑1 (MQ=255)
                                                                                                                                        aCCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCGCACCGGGATAAACa         <  6:357799/148‑1 (MQ=255)
                                                                                                                                        aCCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCGCACCGGGATAAACa         <  4:281618/148‑1 (MQ=255)
                                                                                                                                          cGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTATCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCACCATAACCCGGCAATGCTCGCACCGGGATAAACAGCg      <  2:98828/149‑1 (MQ=255)
                                                                                                                                            tttGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCGCACCGGGATAAACAGCg      >  5:207520/1‑147 (MQ=255)
                                                                                                                                            tttGTCCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCGCACCGGGATAAACAGCgta    <  6:234448/149‑1 (MQ=255)
                                                                                                                                            tttGTCCCGGGAACGCCGGACGGACGCAGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCAGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCAATAACCCGGCAATGCTCGAACCGGGATAAACAGCgta    <  6:229117/149‑1 (MQ=255)
                                                                                                                                               gtgCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCGCACCGGGATAAACAGCgtagt  <  6:327951/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCTCACCAGGCTGTAACGCCGTACCGTAATTAACCGTTTGTGCCGGGAACGCCGGACGGACGCCGATTTTCCACTGCTTATGCTGATCGGCAACGGTTTCACCCGGTAACGCTAACCCGCTAATGGTGGCCTGAATTTCTCCATCGCCATAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGT  >  NC_000913/2648757‑2649047

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