Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,975,672 C→T 10.4% intergenic (+616/‑528) 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,6720CT10.4% 129.1 / 5.2 48intergenic (+616/‑528)PP_5727/yifBhypothetical protein/ATP‑dependent chelatase
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base T (3/2);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTTTTTGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAA  >  NC_002947/5975530‑5975819
                                                                                                                                                |                                                                                                                                                   
gTTTTTGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCta                                                                                                                                                >  2:104670/1‑148 (MQ=255)
   tttGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATaa                                                                                                                                             <  5:58577/148‑1 (MQ=255)
   tttGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATaa                                                                                                                                             <  6:159445/148‑1 (MQ=255)
         aaaaTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAGAAGTGAGGAATTATAAATATTGAACCTGGAAAAGATACGATGCGTATAAGCGGGCCATACGAATACCTGCAAAACATATATAatattt                                                                                                                                        >  1:98434/1‑147 (MQ=255)
          aaaTAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt                                                                                                                                     <  2:236487/149‑1 (MQ=255)
                  gCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTaa                                                                                                                             <  2:225320/149‑1 (MQ=255)
                   cGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAg                                                                                                                            <  6:65608/149‑1 (MQ=255)
                    ggaggaTTCTCGAGAGAGTAAAATAAATTCAGCATGCTCCATCAATTCAAAAGTGA‑‑AATTATACATATTGAACCTGGAAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGt                                                                                                                           <  3:3907/149‑1 (MQ=255)
                    ggaggaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGt                                                                                                                           <  1:104670/149‑1 (MQ=255)
                        gaTTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCt                                                                                                                       >  4:301011/1‑149 (MQ=255)
                               gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGt                                                                                                                 <  7:141026/148‑1 (MQ=255)
                               gagagaGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTg                                                                                                                >  2:198260/1‑149 (MQ=255)
                                agagagGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTgg                                                                                                               <  2:120732/149‑1 (MQ=255)
                                     ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAATATCTTATAGATGTTTTGCAgg                                                                                                                                >  6:199581/1‑126 (MQ=255)
                                     ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAATATCTTATAGATGTTTTGCAgg                                                                                                                                <  5:199581/127‑2 (MQ=255)
                                     ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa                                                                                                          <  2:157339/149‑1 (MQ=255)
                                     ggAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTa                                                                                                          >  3:73641/1‑149 (MQ=255)
                                         aataaatTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTAcaca                                                                                                      <  1:198260/149‑1 (MQ=255)
                                         aataaatTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTAcaca                                                                                                      >  3:176782/1‑149 (MQ=255)
                                               ttCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt                                                                                                                                     >  4:201276/1‑112 (MQ=255)
                                               ttCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATAtttt                                                                                                                                     <  3:201276/112‑1 (MQ=255)
                                                   gCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGcc                                                                                            >  8:65769/1‑149 (MQ=255)
                                                     aTGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCg                                                                                           >  6:97750/1‑148 (MQ=255)
                                                               aaTGCACAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATCCGTATAAGCGGGCCATAAGAAGACCTGCAAAACATATATAAGATATTTTGAATTTTAGTGCCTCGTAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatc                                                                                 <  2:98434/148‑1 (MQ=255)
                                                                 tGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcatc                                                                              <  4:176782/149‑1 (MQ=255)
                                                                  gCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGatcat                                                                               >  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                                                                                                                                    >  7:284870/1‑54 (MQ=255)
                                                                                                            gATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTg                                                                                                                                    <  8:284870/54‑1 (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                        <  4:258801/147‑1 (MQ=255)
                                                                                                                           gCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCtgtg                        >  3:258801/1‑147 (MQ=255)
                                                                                                                            ccATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTc                     >  5:234379/1‑149 (MQ=255)
                                                                                                                                aCGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCgg                  >  8:148530/1‑148 (MQ=255)
                                                                                                                                    ttaacTGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta                                              >  1:70729/5‑116 (MQ=38)
                                                                                                                                    ttaacTGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAagcta                                              <  2:70729/112‑1 (MQ=38)
                                                                                                                                         tGCAAAATATCTTATAGATGTTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACg        >  5:280327/1‑149 (MQ=255)
                                                                                                                                              aaCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTa   <  7:65769/149‑1 (MQ=255)
                                                                                                                                               aCATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTaa  <  7:148530/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTTTTTGTAAAAATAAGAGCGGAGGATTCTCGAGAGAGGAAAATAAATTCAGCATGCTCCATCAATGCAAAAGTGA‑‑AATTATACATATTGAACCTGGCAAAGATACGATGCGTATAAGCGGGCCATACGAAGACCTGCAAAACATCTATAAGATATTTTGCAGTTAAGTGCCTCGGAGTGGCTCTACACAGAAAAGGGCCGAATTTTGATCATCAGACTTGAGCACGAGTGCCATAGCAGAAGCTAGCTTGCGACTGATGGGGCTGTGCTCCGGTACCGAAACGCATTAA  >  NC_002947/5975530‑5975819

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: ‑