Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 838,311 +G coding (146/267 nt) ybiI ← DksA‑type zinc finger protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913838,3081.G97.4% 132.7 / ‑4.0 39coding (149/267 nt)ybiIDksA‑type zinc finger protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (18/20);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGA  >  NC_000913/838161‑838450
                                                                                                                                                    |                                                                                                                                              
ataCTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACggg                                                                                                                                             <  1:612607/151‑1 (MQ=255)
  aCTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACggggc                                                                                                                                           >  1:651076/1‑149 (MQ=255)
   cTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACggggcc                                                                                                                                          >  2:608798/1‑148 (MQ=255)
   cTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACggggcc                                                                                                                                          >  2:111436/1‑148 (MQ=255)
    tCGCTATGCATCTCGCCAGATGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACggggcct                                                                                                                                         >  2:535136/1‑147 (MQ=255)
        tATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCgg                                                                                                                                      >  1:671902/1‑151 (MQ=255)
         aTGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCg                                                                                                                                       <  2:87188/149‑1 (MQ=255)
           gCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGAt                                                                                                                                   <  1:208309/151‑1 (MQ=255)
                 cGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGc                                                                                                                             >  2:30446/1‑151 (MQ=255)
                           aGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCActct                                                                                                                   <  1:92147/151‑1 (MQ=255)
                                  gTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACAtt                                                                                                            >  2:82991/1‑151 (MQ=255)
                                  gTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACAtt                                                                                                            >  2:23597/1‑151 (MQ=255)
                                         gCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAg                                                                                                     <  1:185957/151‑1 (MQ=255)
                                           tGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGc                                                                                                   >  2:211254/1‑151 (MQ=255)
                                              cTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCttt                                                                                                >  1:671678/1‑151 (MQ=255)
                                                  cGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTcgcc                                                                                            >  1:660928/1‑151 (MQ=255)
                                                  cGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTcgc                                                                                             <  2:612921/150‑1 (MQ=255)
                                                           gCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAAt                                                                                   <  2:398488/151‑1 (MQ=255)
                                                             gATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAAttt                                                                                 <  1:240668/151‑1 (MQ=255)
                                                             gATTATATCCTATATAAGCTGGTTTTTGTAAAACTTTCTCCTGCTGACAATGAATACATAAGCGCACTCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGa                                                                                     <  2:386080/147‑1 (MQ=255)
                                                                 atatCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCAcc                                                                             >  2:566364/1‑151 (MQ=255)
                                                                 atatCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCAc                                                                              <  1:246075/150‑1 (MQ=255)
                                                                     ccTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCg                                                                          >  2:581477/1‑150 (MQ=255)
                                                                        gTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAgc                                                                       <  1:188637/150‑1 (MQ=255)
                                                                            aaGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCggg                                                                  <  1:34043/151‑1 (MQ=255)
                                                                                       tGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCt                                                       >  1:288640/1‑151 (MQ=255)
                                                                                         tAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTc                                                     >  2:262209/1‑151 (MQ=255)
                                                                                          aaaTCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTc                                                     <  1:19643/150‑1 (MQ=255)
                                                                                                ttCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGt                                               <  1:328622/150‑1 (MQ=255)
                                                                                                     ctgctgACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTGCTGtt                                         >  1:500519/1‑151 (MQ=255)
                                                                                                        ctgACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCTCTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGa                                       <  2:249824/150‑1 (MQ=255)
                                                                                                        ctgACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGAt                                      >  2:596571/1‑151 (MQ=255)
                                                                                                                    aCATAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGa                           <  1:111438/150‑1 (MQ=255)
                                                                                                                       tAAGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGc                       >  1:563391/1‑151 (MQ=255)
                                                                                                                        aaGCGCACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCg                      <  1:597818/151‑1 (MQ=255)
                                                                                                                          gcgcACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTc                    >  1:663277/1‑151 (MQ=255)
                                                                                                                           cgcACGCCAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCa                   <  2:234970/151‑1 (MQ=255)
                                                                                                                                 ccAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTa             <  1:632419/151‑1 (MQ=255)
                                                                                                                                  cAGGAATGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGtt              >  2:564219/1‑149 (MQ=255)
                                                                                                                                  cAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGc           >  2:211096/1‑151 (MQ=255)
                                                                                                                                       aTGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCaa       <  1:580898/151‑1 (MQ=255)
                                                                                                                                        tGGCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAAc      >  2:91106/1‑151 (MQ=255)
                                                                                                                                         ggCTTCCCGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAAcc     >  2:247896/1‑151 (MQ=255)
                                                                                                                                               ccGACGGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCACGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGa  <  1:127225/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGA  >  NC_000913/838161‑838450

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: