Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,975,684 A→G 10.2% intergenic (+628/‑516) PP_5727 → / → yifB hypothetical protein/ATP‑dependent chelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,975,6840AG10.2% 128.0 / 5.6 49intergenic (+628/‑516)PP_5727/yifBhypothetical protein/ATP‑dependent chelatase
Reads supporting (aligned to +/- strand):  ref base A (19/25);  new base G (3/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGAGCGC  >  NC_002947/5975540‑5975831
                                                                                                                                                |                                                                                                                                                   
aaaTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt                                                                                                                                                 <  2:236487/149‑1 (MQ=255)
        gCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTaa                                                                                                                                         <  2:225320/149‑1 (MQ=255)
         cGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAg                                                                                                                                        <  6:65608/149‑1 (MQ=255)
          ggaggaTTCTCGAGAGAGTAAAATAAATTCAGCATGCTCCATCAATTCAAAAGTGAAATTATACATATTGAACCTGGAAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGt                                                                                                                                       <  3:3907/149‑1 (MQ=255)
          ggaggaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGt                                                                                                                                       <  1:104670/149‑1 (MQ=255)
              gaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCt                                                                                                                                   >  4:301011/1‑149 (MQ=255)
                     gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGt                                                                                                                             <  7:141026/148‑1 (MQ=255)
                     gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTg                                                                                                                            >  2:198260/1‑149 (MQ=255)
                      agagagGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTgg                                                                                                                           <  2:120732/149‑1 (MQ=255)
                           ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAATATCTTATAGATGTTTTGCAgg                                                                                                                                            >  6:199581/1‑126 (MQ=255)
                           ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAATATCTTATAGATGTTTTGCAgg                                                                                                                                            <  5:199581/127‑2 (MQ=255)
                           ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa                                                                                                                      <  2:157339/149‑1 (MQ=255)
                           ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa                                                                                                                      >  3:73641/1‑149 (MQ=255)
                               aataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTAcaca                                                                                                                  <  1:198260/149‑1 (MQ=255)
                               aataaatTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTAcaca                                                                                                                  >  3:176782/1‑149 (MQ=255)
                                     ttCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt                                                                                                                                                 <  3:201276/112‑1 (MQ=255)
                                     ttCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt                                                                                                                                                 >  4:201276/1‑112 (MQ=255)
                                         gCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGcc                                                                                                        >  8:65769/1‑149 (MQ=255)
                                           aTGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCg                                                                                                       >  6:97750/1‑148 (MQ=255)
                                                     aaTGCACAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATCCGTATAAGCGGGCCATAAGAAGACCTGCAAAACATATATAAGATATTTTGAATTTTAGTGCCTCGTAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatc                                                                                             <  2:98434/148‑1 (MQ=255)
                                                       tGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc                                                                                          <  4:176782/149‑1 (MQ=255)
                                                        gCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcat                                                                                           >  2:250155/1‑147 (MQ=255)
                                                                         cATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTg                                                                         <  3:45394/148‑1 (MQ=255)
                                                                            atTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAACACGAGTGCCa                                                                      >  6:7190/1‑148 (MQ=255)
                                                                                    tGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc                                                                                          >  8:157260/1‑120 (MQ=255)
                                                                                    tGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc                                                                                          <  7:157260/120‑1 (MQ=255)
                                                                                      gCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagct                                                           <  1:250155/149‑1 (MQ=255)
                                                                                      gCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagc                                                            >  5:275016/1‑148 (MQ=255)
                                                                                        aaaGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagctag                                                         <  3:301011/149‑1 (MQ=255)
                                                                                                gATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTg                                                                                                                                                <  8:284870/54‑1 (MQ=255)
                                                                                                gATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTg                                                                                                                                                >  7:284870/1‑54 (MQ=255)
                                                                                                 aTGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACt                                                >  6:17155/1‑149 (MQ=255)
                                                                                                    cGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGAt                                             <  5:97750/149‑1 (MQ=255)
                                                                                                    cGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAACACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGAt                                             <  5:7190/149‑1 (MQ=255)
                                                                                                            cGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgt                                     <  6:275016/149‑1 (MQ=255)
                                                                                                               gCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgtg                                    >  3:258801/1‑147 (MQ=255)
                                                                                                               gCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgtg                                    <  4:258801/147‑1 (MQ=255)
                                                                                                                ccATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTc                                 >  5:234379/1‑149 (MQ=255)
                                                                                                                    aCGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCgg                              >  8:148530/1‑148 (MQ=255)
                                                                                                                        ttaacTGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta                                                          <  2:70729/112‑1 (MQ=38)
                                                                                                                        ttaacTGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta                                                          >  1:70729/5‑116 (MQ=38)
                                                                                                                             tGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACg                    >  5:280327/1‑149 (MQ=255)
                                                                                                                                  aaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTa               <  7:65769/149‑1 (MQ=255)
                                                                                                                                   aCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTaa              <  7:148530/149‑1 (MQ=255)
                                                                                                                                         ataAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt        >  7:267000/1‑149 (MQ=255)
                                                                                                                                           aaGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGAt        >  5:20813/1‑147 (MQ=255)
                                                                                                                                            aGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGa                                                  >  1:202466/1‑104 (MQ=255)
                                                                                                                                            aGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGa                                                  <  2:202466/104‑1 (MQ=255)
                                                                                                                                            aGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGa      <  6:20813/148‑1 (MQ=255)
                                                                                                                                              atatTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCAt                 >  4:47435/1‑135 (MQ=255)
                                                                                                                                              atatTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCAt                 <  3:47435/135‑1 (MQ=255)
                                                                                                                                               tatTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGAgcgc  <  2:201660/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGAAATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAACCAGATGAGCGC  >  NC_002947/5975540‑5975831

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