Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,773,862 G→A 8.5% A189T (GCC→ACC)  PP_5743 → TonB‑dependent receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,773,8620GA 8.5% 131.4 / 3.2 47A189T (GCC→ACC) PP_5743TonB‑dependent receptor protein
Reads supporting (aligned to +/- strand):  ref base G (18/25);  new base A (2/2);  total (20/27)
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.65e-01

CAGCGTGTACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCA  >  NC_002947/3773717‑3774001
                                                                                                                                                 |                                                                                                                                           
cAGCGTGTACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCACCTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAagc                                                                                                                                            >  8:170800/1‑147 (MQ=255)
 aGCGTGTACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAagccag                                                                                                                                         <  6:204661/149‑1 (MQ=255)
  gCGTGTACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGCCCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCGTCCGCACGGCCACCGAACAagccagc                                                                                                                                        >  8:119341/1‑149 (MQ=255)
   cGTGTACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTACGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAagccagcc                                                                                                                                       >  2:224080/1‑149 (MQ=255)
     tgtACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCt                                                                                                                                      >  1:102902/1‑148 (MQ=255)
       tACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACTTCCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGAt                                                                                                                                   <  5:193524/149‑1 (MQ=255)
       tACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGAt                                                                                                                                   <  1:147457/149‑1 (MQ=255)
        aCAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGAtt                                                                                                                                  >  5:50736/1‑149 (MQ=255)
          aaaGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTgg                                                                                                                                >  7:91932/1‑149 (MQ=255)
          aaaGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTgg                                                                                                                                >  7:88144/1‑149 (MQ=255)
          aaaGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTgg                                                                                                                                >  5:267571/1‑149 (MQ=255)
                      ccgcagAAGTTCGGCAACACCTACGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGa                                                                                                                    <  1:224080/149‑1 (MQ=255)
                         cagAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGAttt                                                                                                                 <  8:166228/149‑1 (MQ=255)
                         cagAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGAttt                                                                                                                 <  2:102902/149‑1 (MQ=255)
                             aGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGAc                                                                                                             <  4:262118/149‑1 (MQ=255)
                               ttCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTa                                                                                                           >  5:115064/1‑149 (MQ=255)
                                          cTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACAcc                                                                                                  >  2:62032/1‑147 (MQ=255)
                                          cTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCa                                                                                                <  1:62032/149‑1 (MQ=255)
                                            tCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAgg                                                                                              >  4:229319/1‑149 (MQ=255)
                                                              cGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATgg                                                                            <  2:97821/149‑1 (MQ=255)
                                                                cccAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGc                                                                          <  6:50736/149‑1 (MQ=255)
                                                                       ttGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAACCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCAt                                                                                   <  3:249727/133‑1 (MQ=255)
                                                                       ttGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAACCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCAt                                                                                   >  4:249727/1‑133 (MQ=255)
                                                                       ttGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAACCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCAt                                                                                   <  3:145263/133‑1 (MQ=255)
                                                                       ttGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAACCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCAt                                                                                   >  4:145263/1‑133 (MQ=255)
                                                                          ccGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATAc                                                                <  6:267571/149‑1 (MQ=255)
                                                                           cGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCa                                                                                                <  3:173862/116‑1 (MQ=255)
                                                                           cGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCa                                                                                                >  4:173862/1‑116 (MQ=255)
                                                                                        aaCCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTg                                                  <  5:96510/149‑1 (MQ=255)
                                                                                          ccaccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                >  1:150403/1‑149 (MQ=255)
                                                                                          ccaccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCACGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                <  8:88144/149‑1 (MQ=255)
                                                                                          ccaccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCACGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                <  8:91932/149‑1 (MQ=255)
                                                                                           caccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCa                                                           <  2:161014/137‑1 (MQ=255)
                                                                                                    gAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACTGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCt                                      >  7:104157/1‑149 (MQ=255)
                                                                                                       cTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCt                                                   <  8:96658/133‑1 (MQ=255)
                                                                                                       cTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCt                                                   >  7:96658/1‑133 (MQ=255)
                                                                                                          cAACTCACCGGCGGCAGCTTCCGCACGGCCACTGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCa                                <  8:104157/149‑1 (MQ=255)
                                                                                                          cAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCa                                <  6:115064/149‑1 (MQ=255)
                                                                                                                aCCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                <  6:271373/127‑1 (MQ=255)
                                                                                                                aCCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                >  5:271373/1‑127 (MQ=255)
                                                                                                                aCCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTg                           <  2:150403/148‑1 (MQ=255)
                                                                                                                           cTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATg               <  7:170800/149‑1 (MQ=255)
                                                                                                                            ttCCGCACGGCCACCGAACAAGCCAGCCTGATT‑GCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGc             >  6:54000/1‑149 (MQ=255)
                                                                                                                              ccGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGc             <  5:72126/148‑1 (MQ=255)
                                                                                                                                gCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTa          <  5:227289/149‑1 (MQ=255)
                                                                                                                                    ggCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACa                                                                               <  3:218483/76‑1 (MQ=255)
                                                                                                                                    ggCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACa                                                                               >  4:218483/1‑76 (MQ=255)
                                                                                                                                       cACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACGCCCAGGCACGTGCCCAATACATGGGCCTATGCATACTGGCACCGTTTATGGTGTAACGGTAAGTAGCAGCCTGCGCGGGGTAATTGCTTAGCAGTg    >  1:202763/1‑148 (MQ=255)
                                                                                                                                          cGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTg    >  3:97977/1‑145 (MQ=255)
                                                                                                                                           gAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCa  >  3:247028/1‑146 (MQ=255)
                                                                                                                                           gAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCa  <  4:247028/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CAGCGTGTACAAAGGCCCGCAGCCGCAGAAGTTCGGCAACACCTTCGCGGCCATCGACCTGACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCA  >  NC_002947/3773717‑3774001

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: ‑