Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 863,898 T→A 12.2% intergenic (+15/+21) hemH → / ← uraA ferrochelatase/uracil permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947863,8980TA12.2% 113.2 / 6.4 48intergenic (+15/+21)hemH/uraAferrochelatase/uracil permease
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base A (4/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCC  >  NC_002947/863755‑864044
                                                                                                                                               |                                                                                                                                                  
gAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCg                                                                                                                                               <  7:162491/149‑1 (MQ=255)
gAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCg                                                                                                                                               >  5:8413/1‑149 (MQ=255)
 aaGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCgg                                                                                                                                              >  5:200835/1‑149 (MQ=255)
   gagaTAGGCCTGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCGGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGc                                                                                                                                             >  7:188954/1‑148 (MQ=255)
           cATGCGTGGCCGGGAGCTGTTTGTCGAGGGGGGCGGCCTGTCGCTGGGGCTGGTGCCGTGCCTGAATGATCACCCGTAGTGGGTTAGGGTGTTGGCTTATATGTGTGAAACGGCCTGAATCCATGGGGCCGCTTTGCGGCCCCTTCCGa                                                                                                                                    <  6:200835/149‑1 (MQ=255)
             tGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCTTc                                                                                                                                       >  2:190927/1‑144 (MQ=255)
             tGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCTTc                                                                                                                                       <  1:190927/144‑1 (MQ=255)
                                 gTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTCGCGGCCCCTTCCGACCGGTCATACATCCAGcttctt                                                                                                              >  8:151324/1‑149 (MQ=255)
                                    gAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtgtg                                                                                                           >  3:256816/1‑149 (MQ=255)
                                         gggcggCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTcc                                                                                                      <  2:102526/149‑1 (MQ=255)
                                         gggcggCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTCGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTcc                                                                                                      <  7:151324/149‑1 (MQ=255)
                                               ccAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGt                                                                                                <  3:87771/149‑1 (MQ=255)
                                                       tggtgctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTAccc                                                                                        >  7:229032/1‑149 (MQ=255)
                                                           gctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCgg                                                                                      >  5:199613/1‑147 (MQ=255)
                                                             tggtgcGGTGCCTTAATGAGCACCCGGAGTGGGTGAGGTGGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGGTTCCAGCTGTCATTACCCGGCAGa                                                                                  <  6:8413/149‑1 (MQ=255)
                                                                  ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCa                                                                              <  8:229032/148‑1 (MQ=255)
                                                                  ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAg                                                                             >  4:201061/1‑149 (MQ=255)
                                                                    gtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCa                                                                                    <  5:132170/140‑1 (MQ=255)
                                                                    gtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCa                                                                                    >  6:132170/1‑140 (MQ=255)
                                                                                tCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATAGAAGGGGCCGCTTTGCGGCCCCTTCCGAACGGTCATACATACAGCTTCTTGTGCTTCCAGATGTAATTACACGGCAGAATCAGGTTCAGGGCAATAg                                                               >  5:227983/1‑149 (MQ=255)
                                                                                  aCCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGc                                                              >  3:214265/1‑148 (MQ=255)
                                                                                        agtttctGAGTGTGCTGGCTGATATGTTTGAACAGTCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATg                                                       <  6:227983/143‑1 (MQ=255)
                                                                                             gTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGca                                                  >  8:147746/1‑149 (MQ=255)
                                                                                                   gtgtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGc                                            >  7:296352/1‑149 (MQ=255)
                                                                                                   gtgtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGATTCTTGTGATTCCAGCTGTCATTACCCGGCAGAATAAGTTTCAGGGCAATAGCCACAATGCCGCACAggtc                                            >  8:303199/1‑147 (MQ=255)
                                                                                                         gCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTGATACATACAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGACACAATGGCGCAAAGGGCGATGca                                      >  3:14787/1‑148 (MQ=255)
                                                                                                          cTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGccc                                     >  4:285040/1‑149 (MQ=255)
                                                                                                           tGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCt                                    >  1:265769/1‑149 (MQ=255)
                                                                                                           tGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCATTGCGGCCCCTTACGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGACACAATGGCGCACAGGGCGATGCCCt                                    >  5:172940/1‑149 (MQ=255)
                                                                                                            gATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt                                   >  6:140184/1‑149 (MQ=255)
                                                                                                                tgtgttAATAGGTCTGTATCGAAGGGGCCGCTTTGCGGCCCCTTACGACCGGTGATAAATACAGCTTCTTGTGCTTCCAGCTGGAATTACCCGGCAGAATCAGGTTCAGTTCAATATGCACAATGGCGAACAGGGCGATGCCCTTCAg                                <  8:188954/148‑1 (MQ=255)
                                                                                                                    tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGcccc                           <  3:201061/149‑1 (MQ=255)
                                                                                                                      aaaaGGCCTGAATCGAAGGTGCCGCTTTGCGGCCCCTTCCGACCGGTCATACTTCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                          <  6:172940/148‑1 (MQ=255)
                                                                                                                      aaaaGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGTTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                          <  2:265769/148‑1 (MQ=255)
                                                                                                                      aaaaGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                          <  6:116184/148‑1 (MQ=255)
                                                                                                                       aaaGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGGCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAg                         >  6:73256/1‑148 (MQ=255)
                                                                                                                           gCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATc                    <  4:214265/149‑1 (MQ=255)
                                                                                                                                aaTCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCgg                 <  6:199613/147‑1 (MQ=255)
                                                                                                                                    gAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCa                      >  1:288372/1‑138 (MQ=255)
                                                                                                                                    gAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCa                      <  2:288372/138‑1 (MQ=255)
                                                                                                                                    gAAGGGGCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATg                                        <  5:252942/120‑1 (MQ=255)
                                                                                                                                    gAAGGGGCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCATCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATg                                        >  6:252942/1‑120 (MQ=255)
                                                                                                                                     aaGGGGCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCa          >  6:263243/1‑149 (MQ=255)
                                                                                                                                       ggggCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAcc        <  4:256816/149‑1 (MQ=255)
                                                                                                                                       ggggCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAcc        >  2:225026/1‑149 (MQ=255)
                                                                                                                                        gggCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCg       <  5:140184/149‑1 (MQ=255)
                                                                                                                                            cGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGc   >  7:283386/1‑149 (MQ=255)
                                                                                                                                             gCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGcc  >  4:141043/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCC  >  NC_002947/863755‑864044

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