Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,676,371 C→T 7.6% A660A (GCG→GCA yhjG ← putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,676,3710CT 7.6% 160.0 / 2.9 54A660A (GCG→GCAyhjGputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base T (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTC  >  NC_000913/3676233‑3676517
                                                                                                                                          |                                                                                                                                                  
tttGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTa                                                                                                                                          >  3:181404/1‑149 (MQ=255)
  tGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTAc                                                                                                                                         <  3:143845/148‑1 (MQ=255)
    ccAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGTCTACCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGAGTTACCAg                                                                                                                                      >  8:35874/1‑149 (MQ=255)
     caacaaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGc                                                                                                                                     <  7:194455/149‑1 (MQ=255)
               aaTCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCa                                                                                                                              <  3:66209/146‑1 (MQ=255)
                   tGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCg                                                                                                                        >  6:116680/1‑148 (MQ=255)
                    gTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCg                                                                                                                        >  8:204995/1‑147 (MQ=255)
                     tttAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgca                                                                                                                      <  4:181404/148‑1 (MQ=255)
                      ttAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgcag                                                                                                                     >  1:23771/1‑148 (MQ=255)
                            gCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAg                                                                                                               <  5:116680/148‑1 (MQ=255)
                             ccGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGca                                                                                                             >  7:104497/1‑149 (MQ=255)
                             ccGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGGGTTACCAGCGTTGCCAGTGCCGCAGCAACAGc                                                                                                              <  5:190455/148‑1 (MQ=255)
                                      aGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGc                                                                                                                          <  6:241912/127‑1 (MQ=255)
                                      aGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGc                                                                                                                          >  5:241912/1‑127 (MQ=255)
                                             cGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCg                                                                                             >  2:248088/1‑149 (MQ=255)
                                             cGGCAATTATCATCACTTCTTAATCTGCGACAAAATCGTCAGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCAGCCGCCGGAGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCg                                                                                             >  6:150435/1‑149 (MQ=255)
                                                             ttcttcATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCAcccc                                                                             >  4:169844/1‑149 (MQ=255)
                                                             ttcttcATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCAccc                                                                              >  5:100731/1‑148 (MQ=255)
                                                                         gACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGAtt                                                                 >  3:117932/1‑149 (MQ=255)
                                                                            aaaaTCGTCCTACACTGATTAGGCTGCCCTTCGGATGGTGAGATCAGTGCCAGTATCGCCGCGGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACGACTGGGAATCAGAGGTCCGGCGTTGACCCCAGCCTGCGGAttttt                                                              <  5:155942/149‑1 (MQ=255)
                                                                            aaaaTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGAtttt                                                               >  1:34063/1‑148 (MQ=255)
                                                                             aaaTGGTCCTACACTGATTAGCCTCCACTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGGGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTa                                                             <  5:189186/149‑1 (MQ=255)
                                                                             aaaTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGAttttt                                                              >  8:87496/1‑148 (MQ=255)
                                                                               aTCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTaaa                                                           <  1:248088/149‑1 (MQ=255)
                                                                                tCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAAc                                                          <  7:204995/149‑1 (MQ=255)
                                                                                      gACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCTGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTcccc                                                    >  4:158475/1‑149 (MQ=255)
                                                                                      gACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCTGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTcccc                                                    >  2:5917/1‑149 (MQ=255)
                                                                                       acacTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTccccc                                                   <  8:104497/149‑1 (MQ=255)
                                                                                         acTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcac                                                                                                            >  5:40647/1‑90 (MQ=255)
                                                                                         acTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcac                                                                                                            <  6:40647/90‑1 (MQ=255)
                                                                                             aTTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGcaca                                              >  7:58684/1‑148 (MQ=255)
                                                                                                      cccTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAgcggcg                                    <  7:87496/149‑1 (MQ=255)
                                                                                                      cccTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAgcggc                                     >  7:60022/1‑148 (MQ=255)
                                                                                                       ccTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTcc                                                      >  8:177874/1‑130 (MQ=255)
                                                                                                       ccTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTcc                                                      <  7:177874/130‑1 (MQ=255)
                                                                                                              aaGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAgt                            >  8:140078/1‑149 (MQ=255)
                                                                                                                ggTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAgtgt                          <  4:182794/149‑1 (MQ=255)
                                                                                                                 gTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAgtgt                          >  7:73886/1‑148 (MQ=255)
                                                                                                                  tGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGa                        <  7:59322/149‑1 (MQ=255)
                                                                                                                  tGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGa                        <  8:60022/149‑1 (MQ=255)
                                                                                                                  tGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGa                        <  5:126881/149‑1 (MQ=255)
                                                                                                                   gagaTCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGa                        >  5:99684/1‑148 (MQ=255)
                                                                                                                         aGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCCCGGGCAATCAGCGGTCCGGCTGTCACCCCCGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATcc                  <  2:23771/148‑1 (MQ=255)
                                                                                                                         aGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATcc                  <  6:100731/148‑1 (MQ=255)
                                                                                                                          gTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCCCGGGCAATCAGCGGTCCGGCTTTCACCCCCGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCg                 >  1:168358/1‑148 (MQ=255)
                                                                                                                              cAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTc            <  2:201965/149‑1 (MQ=255)
                                                                                                                              cAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTc            >  1:201965/1‑149 (MQ=255)
                                                                                                                               aGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTcc           >  5:231172/1‑149 (MQ=255)
                                                                                                                               aGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCt          >  7:105807/1‑149 (MQ=255)
                                                                                                                                gTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACCTACAGCGGCGAACGCAGTGTGATAATCCGAATTCCt          >  6:195218/1‑149 (MQ=255)
                                                                                                                                gTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCt          >  3:129044/1‑149 (MQ=255)
                                                                                                                                  aaCGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCtt         >  1:104916/1‑148 (MQ=255)
                                                                                                                                   aCGCCGCTGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTa       <  1:5917/149‑1 (MQ=255)
                                                                                                                                   aCGCCGCTGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTa       <  3:158475/149‑1 (MQ=255)
                                                                                                                                      ccgccgccgGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTActc    >  4:101198/1‑149 (MQ=255)
                                                                                                                                        gccgccgGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGATTTCCTTTACtctc  >  3:123780/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
TTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCAATCAGCGGTCCGGCTTTCACCCCAGCCTGCGGATTTTTAAACGTCCCCCGCACATACAGCGGCGAACGCAGTGTGATAATCCGAATTCCTTTACTCTC  >  NC_000913/3676233‑3676517

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: