Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,719,298 Δ1 bp 18.2% coding (382/585 nt) nfuA → Fe/S biogenesis protein NfuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,719,2980A.18.2% 110.7 / 23.0 44coding (382/585 nt)nfuAFe/S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base . (5/3);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

GCGCTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGC  >  NC_002947/2719151‑2719445
                                                                                                                                                   |                                                                                                                                                    
gcgcTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGgccagc                                                                                                                                                    >  1:409407/1‑150 (MQ=255)
         gACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATggcggc                                                                                                                                          >  1:106429/1‑150 (MQ=255)
            tACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAg                                                                                                                                       <  1:583361/150‑1 (MQ=255)
                       ccgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtggt                                                                                                                            >  1:96920/1‑150 (MQ=255)
                       ccgCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCtggt                                                                                                                            <  1:464959/150‑1 (MQ=255)
                               gTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTgg                                                                                                                    <  1:543308/150‑1 (MQ=255)
                                     aGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATg                                                                                                              >  1:215929/1‑150 (MQ=255)
                                      cgtGACCATGAAGGCGCCAAACGCAAATGTGGCTATGGTCAAGAACGACAGCCCTATCACCGACCGTATCAACTACTAACTAAAGACCGCGATCCGCCCCGGCCTAGCCAGCGATGTCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGa                                                                                                             <  2:95399/148‑1 (MQ=255)
                                            cATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCAt                                                                                                       <  2:96920/150‑1 (MQ=255)
                                            cATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCAt                                                                                                       >  1:128953/1‑150 (MQ=255)
                                             aTCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCAtt                                                                                                      >  2:169223/1‑150 (MQ=255)
                                                   gcgcCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTa                                                                                                >  2:551352/1‑150 (MQ=255)
                                                      ccAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGAc                                                                                                                        <  1:348033/123‑1 (MQ=255)
                                                      ccAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGAc                                                                                                                        >  2:348033/1‑123 (MQ=255)
                                                       cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGc                                                                                            <  2:106429/150‑1 (MQ=255)
                                                       cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAggc                                                                                                                                     <  2:244473/109‑2 (MQ=255)
                                                       cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAggc                                                                                                                                     >  1:244473/1‑108 (MQ=255)
                                                        aaaCGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTAGCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCa                                                                                           >  2:548752/1‑150 (MQ=255)
                                                               aaGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggt                                                                                    <  1:184747/150‑1 (MQ=255)
                                                                aGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCggtg                                                                                   <  2:215929/150‑1 (MQ=255)
                                                                     ccGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGc                                                                              <  1:527435/150‑1 (MQ=255)
                                                                     ccGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGc                                                                              >  1:70575/1‑150 (MQ=255)
                                                                                cGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCTCCGACTAGCCAGCCAAGAGACTCACGGAATCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTg                                                                   >  1:75505/1‑150 (MQ=255)
                                                                                 gAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACAGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGACTGGTACAGAGGGGAGAGTTCTTCAGTGCCGTACTGCAGTACGGTGGTGGCTGCAAGGGTTGc                                                                  >  2:173287/1‑150 (MQ=255)
                                                                                  aaGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCg                                                                 <  2:128953/150‑1 (MQ=255)
                                                                                  aaGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCAATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCg                                                                 <  2:70575/150‑1 (MQ=255)
                                                                                         cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGa                                                                                                                <  1:211294/96‑1 (MQ=255)
                                                                                         cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGa                                                                                                                >  2:211294/1‑96 (MQ=255)
                                                                                             aTCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGCGAAGCTTGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTCCCAGGGTTGCGGCCAGGCCGAt                                                      >  1:195703/1‑150 (MQ=255)
                                                                                                       gTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTga                                            <  1:404787/150‑1 (MQ=255)
                                                                                                                  tacctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGAGTTCATTGACATACTAAAATCACACGTCTGAACTCAAGTTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATc                                 >  1:599282/1‑150 (MQ=255)
                                                                                                                    cctacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGa                               >  1:305861/1‑150 (MQ=255)
                                                                                                                      tacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGa                               >  1:400730/1‑148 (MQ=255)
                                                                                                                      tacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGa                               <  2:400730/148‑1 (MQ=255)
                                                                                                                      tacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                             <  2:409407/150‑1 (MQ=255)
                                                                                                                      tacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                             <  1:548752/150‑1 (MQ=255)
                                                                                                                                  tCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCg                 >  1:511925/1‑150 (MQ=255)
                                                                                                                                      ccccGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCg             >  1:489993/1‑150 (MQ=255)
                                                                                                                                         tcaccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCAGGTTTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTAt          >  1:42175/4‑150 (MQ=255)
                                                                                                                                           accTGGCCGCCATGGCTGGCCAGGCGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTc        >  1:253325/2‑150 (MQ=255)
                                                                                                                                            ccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCa                                   >  1:135446/1‑122 (MQ=255)
                                                                                                                                            ccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCa                                   <  2:135446/122‑1 (MQ=255)
                                                                                                                                            ccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTcc       >  1:352752/1‑150 (MQ=255)
                                                                                                                                             cTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCg      >  1:410188/1‑150 (MQ=255)
                                                                                                                                               ggCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGa    <  2:305861/150‑1 (MQ=255)
                                                                                                                                                 ccagccaTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGc  >  2:543243/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGCTGGTCGACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGC  >  NC_002947/2719151‑2719445

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: