Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,084,237 A→G 10.2% K288E (AAA→GAA)  efeB → deferrrochelatase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,084,2370AG10.2% 107.2 / 3.9 40K288E (AAA→GAA) efeBdeferrrochelatase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base G (2/2);  total (19/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 = 9.77e-01

ATCCCGATAGCCAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCC  >  NC_000913/1084091‑1084381
                                                                                                                                                  |                                                                                                                                                
aTCCCGATAGCCATAATGATAAGTTGATGCAAAAAGTGTTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTTCGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGaat                                                                                                                                                <  2:352360/149‑2 (MQ=255)
aTCCCGATAGCCAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGaaa                                                                                                                                                <  1:145344/149‑1 (MQ=255)
aTCCCGATAGCCAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGaaa                                                                                                                                                <  1:150687/149‑1 (MQ=255)
aTCCCGAGAGCCACAATCATAAGAGGATGCAAAAAGAGGTGGGGGTAACGGCAGATCAGCACGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGaaa                                                                                                                                                <  2:125774/149‑1 (MQ=255)
 tCCCGATAGCCAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGa                                                                                                                                                  >  5:5908/1‑146 (MQ=255)
     gATAGCCAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAAca                                                                                                                                           <  5:32812/149‑1 (MQ=255)
      aTAGCCAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGGTTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAAcag                                                                                                                                          <  7:101240/149‑1 (MQ=255)
          ccAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAAca                                                                                                                                           >  7:162837/1‑144 (MQ=255)
          ccAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAAca                                                                                                                                           <  8:162837/144‑1 (MQ=255)
                 gATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGGAAGAACAGCAGACGAtttt                                                                                                                               >  1:340871/1‑149 (MQ=255)
                 gATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGACTCTGGAAGCACAGCAGACGAtttt                                                                                                                               >  1:346332/1‑149 (MQ=255)
                  aTAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGAttttt                                                                                                                              <  6:5908/149‑1 (MQ=255)
                      gTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGcc                                                                                                                          >  8:395051/1‑149 (MQ=255)
                             cAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATaa                                                                                                                   <  7:395051/149‑1 (MQ=255)
                             cAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATaa                                                                                                                   <  3:361141/149‑1 (MQ=255)
                              aaaaaGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAg                                                                                                                  <  2:117836/149‑1 (MQ=255)
                                    tggtgTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAAcc                                                                                                            <  3:128991/149‑1 (MQ=255)
                                           ggTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTgcgc                                                                                                     >  8:322206/1‑149 (MQ=255)
                                            gTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGcgcc                                                                                                    >  1:106250/1‑149 (MQ=255)
                                                    aGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGCAAGAACAGCAGACGATTTTTGGACGTGATAATCAAACCGGTGCGCCGCTGGGaa                                                                                            >  5:20858/1‑149 (MQ=255)
                                                                     cGTGGACAATCTGTGTCAGCTATCAGGAAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAAGGCCGCTTAAAGATCAGAAGACGATTTTTGGCCGTTATAAGCAAACCGGTGCGCCTCTGTGAATGCAGCATGAGCatgat                                                                           <  6:20858/149‑1 (MQ=255)
                                                                              tCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGTCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTgcg                                                                                                      <  6:190076/113‑1 (MQ=255)
                                                                              tCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGTCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTgcg                                                                                                      >  5:190076/1‑113 (MQ=255)
                                                                                ggTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGAtt                                                                 <  8:58351/148‑1 (MQ=255)
                                                                                                   tACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCggaagg                                              >  3:341788/1‑148 (MQ=255)
                                                                                                    aCGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGcgccgc                                                                                                  <  1:256687/95‑1 (MQ=255)
                                                                                                    aCGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGcgccgc                                                                                                  >  2:256687/1‑95 (MQ=255)
                                                                                                     cGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCggaagggaa                                           <  2:106250/149‑1 (MQ=255)
                                                                                                      gCTTGATTCAGTTTCGAGTTGAATTTTGGGACAGAACGCCGCTGGAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCggaagggaag                                          <  2:346332/149‑1 (MQ=255)
                                                                                                      gCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGGAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCggaagggaag                                          <  2:340871/149‑1 (MQ=255)
                                                                                                      gCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCggaagggaag                                          >  3:181543/1‑149 (MQ=255)
                                                                                                       cTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCTGGATGTGCCTGATTACGCCAGCGACCCGgaagggaagg                                         >  6:121912/1‑149 (MQ=255)
                                                                                                          gATTCAGTTTCGAGTGGAATTTTGGGACGGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCggaa                                                <  1:195534/139‑1 (MQ=255)
                                                                                                          gATTCAGTTTCGAGTGGAATTTTGGGACGGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCggaa                                                >  2:195534/1‑139 (MQ=255)
                                                                                                           aTTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGAt                                     >  2:73949/1‑149 (MQ=255)
                                                                                                            ttCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGAGTACGCCAGCGACCCGGAAGGGAAGGTGATc                                    >  4:16633/1‑149 (MQ=255)
                                                                                                                      aGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCTTGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACa                          <  7:175203/149‑1 (MQ=255)
                                                                                                                        tGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGc                        >  8:268052/1‑149 (MQ=255)
                                                                                                                        tGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGc                        <  1:73949/149‑1 (MQ=255)
                                                                                                                        tGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGc                        >  8:54958/1‑149 (MQ=255)
                                                                                                                                  ggACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGc              >  4:357454/1‑149 (MQ=255)
                                                                                                                                      aGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGGGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGc          <  7:322206/149‑1 (MQ=255)
                                                                                                                                              gcTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATcccc  <  7:268052/149‑1 (MQ=255)
                                                                                                                                               cTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTg                              <  3:63779/120‑1 (MQ=255)
                                                                                                                                               cTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTg                              >  4:63779/1‑120 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCCCGATAGCCAGAATGATAAGTTGATGCAAAAAGTGGTGTGGGTAACGGCAGATCAGCAGGAGCCTGCGTGGACAATCGGTGGCAGCTATCAGGCAGTACGCTTGATTCAGTTTCGAGTGGAATTTTGGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCC  >  NC_000913/1084091‑1084381

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