Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,719,306:1 +T 17.4% coding (390/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,3061.T17.4% 126.1 / 29.4 46coding (390/585 nt)nfuAFe/S biogenesis protein NfuA
Reads supporting (aligned to +/- strand):  ref base . (21/17);  new base T (5/3);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGC  >  NC_002947/2719160‑2719453
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                               >  1:128953/1‑150 (MQ=255)
                                   cATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCAt                                                                                                               <  2:96920/150‑1 (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                                                                                                                                             >  1:244473/1‑108 (MQ=255)
                                              cAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCGCCATGGCTGGCCAggc                                                                                                                                             <  2:244473/109‑2 (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                                                                                                                        >  2:211294/1‑96 (MQ=255)
                                                                                cccGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGa                                                                                                                        <  1:211294/96‑1 (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                                     <  1:548752/150‑1 (MQ=255)
                                                                                                             tacAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAgc                                     <  2:409407/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                                           <  2:135446/122‑1 (MQ=255)
                                                                                                                                   ccTGGCCGCCATGGCTGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCa                                           >  1:135446/1‑122 (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)
                                                                                                                                              aTGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAgg    >  1:260426/1‑150 (MQ=255)
                                                                                                                                               tGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAggg   >  2:48907/1‑150 (MQ=255)
                                                                                                                                                ggcggcCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGc  <  2:195703/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GACTACGCCACCGACCGCATGGGTGGCCAGCTGACCATCAAGGCGCCAAACGCCAAGGTGCCGATGGTCAACGAAGACAGCCCGATCAACGAACGTATCAACTACTACCTACAGACCGAGATCAACCCCGGCCTGGCCAGCCATGGCGGCCAGGTGAGCCTGGTGGACGTGGTCGATGACGGCATTGCCGTACTGCAGTTCGGTGGTGGCTGCCAGGGTTGCGGCCAGGCCGATGTCACCCTGAAGGAAGGCATCGAGCGCACCCTGCTCGAGCGTATTCCGGAGCTGAAGGGC  >  NC_002947/2719160‑2719453

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: