Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,822,231 A→T 10.4% intergenic (+31/+27) PP_mr22 → / ← fdxA PrrB_RsmZ ncRNA/ferredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,822,2310AT10.4% 127.6 / 6.0 47intergenic (+31/+27)PP_mr22/fdxAPrrB_RsmZ ncRNA/ferredoxin 1
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base T (3/2);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGA  >  NC_002947/1822084‑1822376
                                                                                                                                                   |                                                                                                                                                 
aTGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAg                                                                                                                                                 <  1:958452/150‑1 (MQ=255)
        gCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCtttt                                                                                                                                         >  1:1139223/1‑150 (MQ=255)
           ttCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAATGGCCTTTTTcg                                                                                                                                      >  1:1486447/1‑150 (MQ=255)
             aaGGATGAAGGCACGGGACAGCGCAAGAAGCGATTTCATCAGGATGATGTTTGCGACAAGCAGGGACTACGGAAAAACAGGTGGGCGGGTCAAACCGCCTCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGTAGTGCCTTTTTcgcg                                                                                                                                    <  2:1497593/149‑1 (MQ=255)
               ggATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCACGCAGAGCTttttacccgac                                                                                                                                  >  2:349014/1‑143 (MQ=255)
                gATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCt                                                                                                                                 >  1:337763/1‑150 (MQ=255)
                         cAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTAt                                                                                                                        >  1:985120/1‑150 (MQ=255)
                             gACATCACAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAgg                                                                                                                    >  1:622606/1‑149 (MQ=255)
                                 tCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttttgtt                                                                                                                                    >  1:1159699/1‑126 (MQ=255)
                                 tCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttttgtt                                                                                                                                    <  2:1159699/130‑5 (MQ=255)
                                   gCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAACAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCa                                                                                                              <  2:770272/150‑1 (MQ=255)
                                        aagCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAACGGCCCGAAATGCATCTTGGGCAGCTGGTCGCTATCATCGCTCCAGGTCg                                                                                                         >  2:835482/1‑150 (MQ=255)
                                         agCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCTCGTAGGAAAACAAAAAAGGCCCGaa                                                                                                                                                  <  2:515360/108‑1 (MQ=255)
                                         agCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCTAGGAGGAAAACAAAAAAGGCCCGaa                                                                                                                                                  >  1:515360/1‑108 (MQ=255)
                                                tcatcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCtt                                                                                                 <  2:1486447/150‑1 (MQ=255)
                                                   tcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGACAATCTTGcc                                                                                              >  1:1535245/1‑150 (MQ=255)
                                                    caGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCa                                                                                             <  2:622606/150‑1 (MQ=255)
                                                       gatgatGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAgcgc                                                                                                                  >  2:789574/1‑126 (MQ=255)
                                                       gatgatGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAgcgc                                                                                                                  <  1:789574/126‑1 (MQ=255)
                                                              tttGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCAt                                                                                   >  2:1342163/1‑150 (MQ=255)
                                                                   gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAt                                                                              >  1:113171/1‑150 (MQ=255)
                                                                            ggACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTg                                                                                                                                <  1:1072161/91‑1 (MQ=255)
                                                                            ggACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTg                                                                                                                                >  2:1072161/1‑91 (MQ=255)
                                                                                    gAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAg                                                             <  2:113171/150‑1 (MQ=255)
                                                                                     aaaaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAgg                                                            <  2:696363/150‑1 (MQ=255)
                                                                                           aTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGc                                                                                                                                              >  2:1034009/1‑62 (MQ=255)
                                                                                           aTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGc                                                                                                                                              <  1:1034009/62‑1 (MQ=255)
                                                                                            tgtgGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCt                                                     <  2:337763/150‑1 (MQ=255)
                                                                                              tgGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATc                                                                                  >  2:1272566/1‑119 (MQ=255)
                                                                                              tgGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATc                                                                                  <  1:1272566/119‑1 (MQ=255)
                                                                                                  cGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTcc                                                                                                               >  2:1560084/1‑86 (MQ=255)
                                                                                                  cGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTcc                                                                                                               <  1:1560084/86‑1 (MQ=255)
                                                                                                      tCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                           >  1:81936/1‑150 (MQ=255)
                                                                                                      tCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                           >  1:629441/1‑150 (MQ=255)
                                                                                                       cAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCTTTTTGCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTg                                          <  1:1635132/150‑1 (MQ=255)
                                                                                                        aaaCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGa                                         >  1:718363/1‑150 (MQ=255)
                                                                                                            cGCCCCTTTTTTTGCTCGTTTGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGTTAtt                                     >  2:1343253/1‑150 (MQ=255)
                                                                                                              ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAg                                   <  1:1290472/150‑1 (MQ=255)
                                                                                                              ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAg                                   <  2:1139223/150‑1 (MQ=255)
                                                                                                                           ccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATc                                                                                  >  2:967302/1‑90 (MQ=255)
                                                                                                                           ccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATc                                                                                  <  1:967302/90‑1 (MQ=255)
                                                                                                                           ccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGc                      <  2:691856/150‑1 (MQ=255)
                                                                                                                           ccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGCCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGc                      >  1:1520006/1‑150 (MQ=255)
                                                                                                                            cGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTc                                                       <  2:1034171/116‑1 (MQ=255)
                                                                                                                            cGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTc                                                       >  1:1034171/1‑116 (MQ=255)
                                                                                                                                        aaaaaGGCCCTTCGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTc         >  2:550488/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGCCCTTCGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTc         >  2:1553415/1‑150 (MQ=255)
                                                                                                                                         aaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCa        <  2:629441/150‑1 (MQ=255)
                                                                                                                                               cccGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGa  >  2:524729/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGA  >  NC_002947/1822084‑1822376

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