Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,063 C→A 12.8% intergenic (+45/‑131) PP_5259 → / → ydcJ LysR family transcriptional regulator/metalloprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,0630CA12.8% 122.9 / 8.1 47intergenic (+45/‑131)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base A (4/2);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

CAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCC  >  NC_002947/6004919‑6005199
                                                                                                                                                |                                                                                                                                        
cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTGCGGGGGCCGGGTAGGCGCGGc                                                                                                                                     >  1:12937/1‑150 (MQ=255)
       ggAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGcc                                                                                                                              >  1:270134/1‑150 (MQ=255)
              aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGTTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGGAGGCGCGGCTTTAGCAGGGAAAg                                                                                                                       >  1:731310/1‑150 (MQ=255)
              aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                       >  1:770343/1‑150 (MQ=255)
              aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTCAGCAAGGCATTCGGTTTCTACGGAGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                       >  2:20559/1‑150 (MQ=255)
                        gTATGCCGAGCATGCGGCAGAGGTGCCAAAGGGGCGGGCGCTGGTGGAATGGATTGGCGCGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGCAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAg                                                                                                             <  2:600644/150‑1 (MQ=255)
                           tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                     <  1:224015/139‑1 (MQ=255)
                           tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                     >  2:224015/1‑139 (MQ=255)
                                     gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCaaaaaa                                                                                                >  1:89756/1‑150 (MQ=255)
                                        gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGCTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaaga                                                                                                                       >  2:334545/1‑122 (MQ=14)
                                        gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGCTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGaaga                                                                                                                       <  1:334545/124‑3 (MQ=14)
                                         cAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGaa                                                                                            <  2:481187/150‑1 (MQ=255)
                                                 ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCATACGCGCGGCGTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGacacac                                                                                    >  1:776729/1‑150 (MQ=255)
                                                       gTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGCCACACTCACtt                                                                              <  2:715106/150‑1 (MQ=255)
                                                          ggggCGCTGGAGGAATGGATTGGCGAGCAGTTGCAGGTTTGCCGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACttttt                                                                           <  2:770343/149‑1 (MQ=255)
                                                           gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCtt                                                                                                                                   <  1:274497/93‑1 (MQ=255)
                                                           gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCtt                                                                                                                                   >  2:274497/1‑93 (MQ=255)
                                                           gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtttttt                                                                          >  1:333085/1‑150 (MQ=255)
                                                                     ggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                <  2:553215/150‑1 (MQ=255)
                                                                      gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAg                                                               <  2:241155/150‑1 (MQ=255)
                                                                        aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCTCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCa                                                             >  2:495099/1‑150 (MQ=255)
                                                                          ggATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCaaa                                                           <  1:299367/150‑1 (MQ=255)
                                                                            aTTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                                     <  2:374116/74‑1 (MQ=255)
                                                                            aTTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                                     >  1:374116/1‑74 (MQ=255)
                                                                                 tgAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGATAAAGCGGCTCCTACGCAGACCGCGTAGGCGCAACTTTAGCCGGGAAAGGGCCGTCGAGGTATATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCAtt                                                    <  2:12937/149‑1 (MQ=255)
                                                                                 cGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAGCGCGTAGGCGCGGCGTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTATCTTGTGAGGAAAGTGCagg                                                    >  1:292762/1‑148 (MQ=255)
                                                                                     cAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTGGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCa                                                >  1:377516/1‑150 (MQ=255)
                                                                                            aGGTGTGAAGGCCTCTTGGCGGATAAAGCCGCTCTTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCGGTCGAGGTTTATCAAAAAAGGAATGACACACTCACTTTTTTTCGTTTGTGAGCAAAGTGCATTCCCAGCTTTCa                                         <  2:331303/150‑1 (MQ=255)
                                                                                              gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATg                                       >  1:611489/1‑150 (MQ=255)
                                                                                               tttGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTTCGGTTAAAGTGCATTCCCAGCTTTCATGt                                      <  1:20559/150‑1 (MQ=255)
                                                                                                   ccgatctcctcGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGTCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGc                                  <  2:420016/145‑1 (MQ=255)
                                                                                                   aCGGCCTCTTCGCGGGTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTTGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAgg                                  >  2:595017/1‑149 (MQ=255)
                                                                                                      gCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTc                               <  2:776729/150‑1 (MQ=255)
                                                                                                      gCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTc                               >  1:54881/1‑150 (MQ=255)
                                                                                                       ccTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                                            <  1:640051/72‑1 (MQ=255)
                                                                                                       ccTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                                            >  2:640051/1‑72 (MQ=255)
                                                                                                          ctTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGaa                           <  2:333085/150‑1 (MQ=255)
                                                                                                           ttCGCGGCTAAAGCCGCTCCTACGGGGTCCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCgagatcgg                          >  1:420016/1‑143 (MQ=255)
                                                                                                           ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAg                                                                                               <  1:478011/81‑1 (MQ=255)
                                                                                                           ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAg                                                                                               >  2:478011/1‑81 (MQ=255)
                                                                                                                        ccGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCTAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGc             >  2:311682/1‑150 (MQ=38)
                                                                                                                          gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa           >  1:422578/4‑150 (MQ=39)
                                                                                                                          gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                >  2:630743/4‑97 (MQ=17)
                                                                                                                          gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                <  1:630743/94‑1 (MQ=17)
                                                                                                                              cTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGa       >  1:426637/1‑150 (MQ=255)
                                                                                                                                aCGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATg     >  1:456562/1‑150 (MQ=255)
                                                                                                                                   gggCCCGCGTGCGCGCTGCTTCAGGCGGGAAAGGGACGTAGAGGTAGATCAAAAAAGTAATGACACACTCACTTTTTTTCGCTTGTGATCAAAGTGTATGCCCAGGTTTCATGTAAGCGTCCGACCCCAACCAGGAAGCTAGCGATGccc  <  2:292762/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
CAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCC  >  NC_002947/6004919‑6005199

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