Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,676,271 A→T 8.7% intergenic (+163/+19) yhjE → / ← yhjG putative MFS transporter; membrane protein/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,2710AT 8.7% 154.9 / 5.7 46intergenic (+163/+19)yhjE/yhjGputative MFS transporter; membrane protein/putative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base T (2/2);  total (24/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 = 9.96e-01

CTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCA  >  NC_000913/3676123‑3676419
                                                                                                                                                    |                                                                                                                                                    
cTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCa                                                                                                                                                      <  1:478821/149‑1 (MQ=255)
     gcgcTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACAt                                                                                                                                                 <  2:561811/149‑1 (MQ=255)
        cTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCg                                                                                                                                              <  3:311461/149‑1 (MQ=255)
         tACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCgg                                                                                                                                             >  2:156390/1‑149 (MQ=255)
            gCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCaa                                                                                                                                          <  3:65553/149‑1 (MQ=255)
                   aGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcat                                                                                                                                   <  2:335435/149‑1 (MQ=255)
                        tGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCActt                                                                                                                              <  1:156390/149‑1 (MQ=255)
                                 aTACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTg                                                                                                                     >  3:611817/1‑149 (MQ=255)
                                  tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGc                                                                                                                    >  1:155109/1‑149 (MQ=255)
                                   aCAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCg                                                                                                                   >  4:87598/1‑149 (MQ=255)
                                     aaTATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGAc                                                                                                                 <  4:611817/149‑1 (MQ=255)
                                       tatatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCaa                                                                                                                                          <  3:596876/122‑1 (MQ=255)
                                       tatatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCaa                                                                                                                                          >  4:596876/1‑122 (MQ=255)
                                              aaTTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTc                                                                                                        >  4:204062/1‑149 (MQ=255)
                                                    cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGacac                                                                                                  >  1:228739/1‑149 (MQ=255)
                                                      aGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCa                                                                                                                                                      >  2:256061/1‑95 (MQ=255)
                                                      aGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCa                                                                                                                                                      <  1:256061/95‑1 (MQ=255)
                                                      aGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTg                                                                                                >  2:353564/1‑149 (MQ=255)
                                                      aGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTg                                                                                                >  4:14263/1‑149 (MQ=255)
                                                          tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTa                                                                                            <  3:14263/149‑1 (MQ=255)
                                                               aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaaaa                                                                                                             <  2:332316/127‑1 (MQ=255)
                                                               aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaaaa                                                                                                             >  1:332316/1‑127 (MQ=255)
                                                                 gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACttcttc                                                                                                                          >  1:549795/1‑112 (MQ=38)
                                                                 gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACttcttc                                                                                                                          <  2:549795/112‑1 (MQ=37)
                                                                 gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTccc                                                                                     >  2:625136/1‑149 (MQ=255)
                                                                      aTAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCg                                                                                <  1:353564/149‑1 (MQ=255)
                                                                        aaGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGa                                                                              >  2:488361/1‑149 (MQ=255)
                                                                             gTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTg                                                                         >  1:343998/1‑149 (MQ=255)
                                                                              ttCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTga                                                                        <  3:87598/149‑1 (MQ=255)
                                                                              ttCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTga                                                                        <  1:488361/149‑1 (MQ=255)
                                                                               tCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTgag                                                                       >  4:193398/1‑149 (MQ=255)
                                                                                 aCGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGAt                                                                     >  3:598366/1‑149 (MQ=255)
                                                                                   gccgcATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCa                                                                   >  1:7869/1‑149 (MQ=255)
                                                                                     cgcATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGt                                                                 <  2:155109/149‑1 (MQ=255)
                                                                                     cgcATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGt                                                                 <  2:343998/149‑1 (MQ=255)
                                                                                             gCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTaa                                                         <  4:598366/149‑1 (MQ=255)
                                                                                                         cacaCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCgg                                             >  3:135919/1‑149 (MQ=255)
                                                                                                            acTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACt                                                                                                 >  3:248488/1‑94 (MQ=255)
                                                                                                            acTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACt                                                                                                 <  4:248488/94‑1 (MQ=255)
                                                                                                                  ccAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAg                                    >  3:284469/1‑149 (MQ=255)
                                                                                                                   caacaaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGc                                   <  3:204062/149‑1 (MQ=255)
                                                                                                                     acaaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGt                                 >  1:3528/1‑149 (MQ=255)
                                                                                                                     acaaTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGt                                 <  2:7869/149‑1 (MQ=255)
                                                                                                                        aTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGc                              <  4:284469/149‑1 (MQ=255)
                                                                                                                                   tttAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgcag                   <  2:3528/149‑1 (MQ=255)
                                                                                                                                      aaTTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCaa                >  2:195671/1‑149 (MQ=255)
                                                                                                                                               aTGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcaccac       >  3:162004/1‑149 (MQ=255)
                                                                                                                                                 gTCTGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCg                                                                                >  2:602567/1‑74 (MQ=255)
                                                                                                                                                 gTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACgg     <  4:135919/149‑1 (MQ=255)
                                                                                                                                                    aGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACAGTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCTACAGCACCACGGGCa  >  2:181580/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCACGGGCA  >  NC_000913/3676123‑3676419

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑