Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,676,271 A→T 15.0% intergenic (+163/+19) yhjE → / ← yhjG putative transporter YhjE/AsmA family protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,676,2710AT15.0% 83.7 / 11.9 41intergenic (+163/+19)yhjE/yhjGputative transporter YhjE/AsmA family protein YhjG
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base T (2/4);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-02

TATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCA  >  NC_000913/3676138‑3676413
                                                                                                                                     |                                                                                                                                              
tATCAGGCCTGAATTTCCATACAATAGATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCgg                                                                                                                                       <  1:399652/143‑1 (MQ=255)
  tCAGGCGTGAATTCCCATACAATCTATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAAGCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCa                                                                                                                                     <  2:325063/143‑1 (MQ=37)
  ccAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCGGCa                                                                                                                                     <  1:698418/142‑1 (MQ=255)
      gCCTGAATTTCCATACAATATATTGAATTCGCCAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACCTCCGGccacta                                                                                                                                 >  1:790252/1‑138 (MQ=255)
          gAATTTTCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcat                                                                                                                             <  2:829014/143‑1 (MQ=255)
                 cacaCAATACATGGAAGTCTGCAGGATTGGCCGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCGAAATCTGTTTAATTGCCGGATGTCAGACATTCGGCAATTATCATCActtct                                                                                                                      <  1:849850/140‑1 (MQ=25)
                 cATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATCATCATAActtca                                                                                                                      >  2:817395/1‑142 (MQ=255)
                        tatagTTAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACCCAGCACACTTTGCCAACAATCTAAATCGGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTg                                                                                                               <  2:210536/143‑1 (MQ=38)
                         atatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGc                                                                                                              <  1:353445/143‑1 (MQ=255)
                          tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAAAAAGCACACTTTGCCAACAATCTAAATCTGTTAAATAGCCGGATGTCAGACACCCGGAAACAATCATAACTTTTTCATcagag                                                                                                             >  2:151115/1‑139 (MQ=25)
                             tGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTGGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAATATCAGCACTGCTGCCTCTGCGACa                                                                                                          >  2:823153/1‑143 (MQ=37)
                                     cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAGACCGTc                                                                                                  >  1:701627/1‑143 (MQ=255)
                                     cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTc                                                                                                  >  1:461584/1‑143 (MQ=255)
                                     cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTc                                                                                                  >  1:104689/1‑143 (MQ=255)
                                      aaGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGATTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTcc                                                                                                 >  2:736051/1‑143 (MQ=255)
                                       aGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCg                                                                                                <  1:853997/143‑1 (MQ=255)
                                        ggATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGa                                                                                               <  2:59718/143‑1 (MQ=255)
                                          aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCgg                                                                                                                                       <  1:776384/101‑1 (MQ=255)
                                          aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCTGACATCCgg                                                                                                                                       >  2:776384/1‑101 (MQ=255)
                                                atgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCACACTTCGCCAACACTCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAGTATAATCACTTCTTCATCTGCGCCACAATCGTCCGAAACAGagc                                                                                       >  2:706896/3‑141 (MQ=25)
                                                   ccGGAGAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGc                                                                                    <  2:767778/143‑1 (MQ=255)
                                                           ggggTTCACTCCGCATCCGGCATGAACAAAGCACACTTTGTCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCtt                                                                            <  2:566731/140‑1 (MQ=255)
                                                                   cgccgccgcCGGCATGAACAAAGCACACGTTGACAACAATCTTAATCTGTTTAATTGCCGGACGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGt                                                                    <  2:735382/143‑1 (MQ=255)
                                                                     ccgcATCCGGCATGAACAGAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTga                                                                  <  2:897939/143‑1 (MQ=255)
                                                                           ccGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAg                                                            <  2:578274/143‑1 (MQ=255)
                                                                            cGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGt                                                           <  1:50497/143‑1 (MQ=255)
                                                                              gCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGc                                                         <  2:396562/143‑1 (MQ=255)
                                                                                   aaCAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGtg                                                    >  1:317928/1‑142 (MQ=255)
                                                                                               tttGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCCCTTCTTCATCTGCGACAAACTCGTCCGCCACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGccgccgccg                                        <  1:823153/143‑1 (MQ=255)
                                                                                                 tGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTGCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGgt                                      <  2:104689/143‑1 (MQ=255)
                                                                                                   ccAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGgtgt                                    >  2:252760/1‑143 (MQ=255)
                                                                                                                 cTGTTTCATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCa                      <  2:88042/143‑1 (MQ=255)
                                                                                                                 cTGTTTAATTGCCGGCGTGCACTCAGCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCa                      <  2:143424/143‑1 (MQ=255)
                                                                                                                  tGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCGTCCCCTTCGGAAGGTGAGATCAGTGCGAGTAACGCGGCCGCCGGTGTTACCAGCGTTGCCAg                     >  1:459887/1‑143 (MQ=255)
                                                                                                                        aTTGCCGGATGTGAGACATCCGGCAATTACCATCACTTCTTCATCGGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgc               <  2:537684/143‑1 (MQ=255)
                                                                                                                        aTTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgc               <  1:736051/143‑1 (MQ=255)
                                                                                                                          tGCCGGATGTCTGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCgcag             >  1:171355/1‑143 (MQ=255)
                                                                                                                            ccGGATGTCAGACATCCGGCAATTATCATCACTTCGTCATATGCGACAAAAACGGCCGACACTGAGTAGCCTCCCCTTCGGAAGGTCAGATCAGTGCCAGTACCGACGCCGCCGGTGATACCAGCGTTTCCCGTGCCgcagca           >  2:155291/1‑143 (MQ=255)
                                                                                                                             aggATGTCAGACAGCCGGCAATTATCATCACTTCTGCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCaa          <  1:125192/142‑1 (MQ=255)
                                                                                                                                ctGTCCGACATCCGGCAATTATCATCACTTCTTCATCGGCGACAAAATCGTCCGACACTGATTAGCCTCCTCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAg       <  2:529584/142‑1 (MQ=255)
                                                                                                                                   tCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcac    >  2:346436/1‑143 (MQ=255)
                                                                                                                                     aGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACCCTGATTAGCCTCCCCTTCGGAAGGCGCGATCAGTGCCAGTAACGCCGCCACCGGTGTTACCAGCTCTGCCCGACCAGCACCAGCAGcacca  >  2:284163/1‑143 (MQ=255)
                                                                                                                                     aGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGcacca  >  2:806074/1‑143 (MQ=255)
                                                                                                                                     aGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAACGTGAGACCCGTGCCAGTAACGCCGCCGCCGCTGTTACCAGCGTTGCCCGTGCCCGCGCAACAGcacca  >  2:843297/1‑143 (MQ=255)
                                                                                                                                     |                                                                                                                                              
TATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGACACTGATTAGCCTCCCCTTCGGAAGGTGAGATCAGTGCCAGTAACGCCGCCGCCGGTGTTACCAGCGTTGCCAGTGCCGCAGCAACAGCACCA  >  NC_000913/3676138‑3676413

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