Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,763,987 Δ1 bp coding (23/1488 nt) sufB ← component of SufBCD Fe‑S cluster assembly scaffold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,763,9870G.98.2% 222.5 / ‑2.6 56coding (23/1488 nt)sufBcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (24/31);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCATTACGCTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGT  >  NC_000913/1763854‑1764119
                                                                                                                                     |                                                                                                                                    
cTCATTACGCTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCtt                                                                                                                                <  1:1232563/139‑1 (MQ=255)
   aTTACGCTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAg                                                                                                                             >  1:1051186/1‑139 (MQ=255)
         cTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTAc                                                                                                                       >  1:1549495/1‑139 (MQ=255)
          ttCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACg                                                                                                                      <  1:924658/139‑1 (MQ=255)
          ttCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACg                                                                                                                      >  2:1292723/1‑139 (MQ=255)
          ttCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACg                                                                                                                      <  1:1613796/139‑1 (MQ=255)
                                           aTCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGCTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCg                                                                                     >  1:1472103/1‑139 (MQ=255)
                                            tCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgc                                                                                    <  2:1180855/139‑1 (MQ=255)
                                              cccTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                  <  2:25567/139‑1 (MQ=255)
                                                  ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                              <  2:1051186/139‑1 (MQ=255)
                                                  ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                              <  2:1472103/139‑1 (MQ=255)
                                                  ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                              <  2:1549495/139‑1 (MQ=255)
                                                       aGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcat                                                                         <  1:1009103/139‑1 (MQ=255)
                                                               ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCATGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                 <  1:839436/139‑1 (MQ=255)
                                                               ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATGCAGCGGGCCGCCGGTCCAGGTTTTGTCATCGTCAGTTCTTCAGTATTACGAGACATAGTACCGCTTATACCCTAAAGCTTTCGCCACAGCCACTTTCATTCTg                                                                 <  1:502516/139‑1 (MQ=255)
                                                                   gCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACg                                                                                                                      >  1:1928272/1‑82 (MQ=255)
                                                                   gCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACg                                                                                                                      <  2:1928272/82‑1 (MQ=255)
                                                                     tAACTGGGTGAATAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACAGCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCttt                                                           <  2:1946130/139‑1 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCGCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                 >  1:1336599/1‑127 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                              <  1:1691105/98‑1 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                              >  2:1691105/1‑98 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACAGTCATTCTg                                                                 <  2:1336599/127‑1 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTa                                                                                                            >  1:1674856/1‑82 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTa                                                                                                            <  2:1674856/82‑1 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCtt                                                            >  1:26386/1‑130 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCtt                                                            <  2:26386/130‑1 (MQ=255)
                                                                                      cttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCtt                                                                                                                                >  1:841059/1‑53 (MQ=255)
                                                                                      cttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCtt                                                                                                                                <  2:841059/53‑1 (MQ=255)
                                                                                            tataATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACAt                                                                                                                <  1:403532/63‑1 (MQ=255)
                                                                                            tataATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACAt                                                                                                                >  2:403532/1‑63 (MQ=255)
                                                                                               aaTTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGAt                                 >  1:787597/1‑139 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTAc                                                                                                           >  1:374849/1‑63 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTAc                                                                                                           <  2:374849/63‑1 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                        >  1:1162071/1‑130 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                        <  2:1162071/130‑1 (MQ=255)
                                                                                                   cAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAg                             <  1:1292723/139‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACg                                                                                                                      <  2:1424286/47‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACg                                                                                                                      >  1:1424286/1‑47 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTAc                                                                                                           <  1:608457/58‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTAc                                                                                                           >  2:608457/1‑58 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                          <  2:1554280/107‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                          >  1:1554280/1‑107 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                >  1:170679/1‑117 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                <  2:170679/117‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                          >  1:502773/1‑139 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATATGATTAAGTcc                          >  2:1332390/1‑139 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCGCAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCt                         >  1:1215654/1‑139 (MQ=255)
                                                                                                              cgGTCCAGGTTTTGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacgaa                  <  1:1332390/139‑1 (MQ=255)
                                                                                                                          tGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                  >  1:1519075/1‑79 (MQ=255)
                                                                                                                          tGACATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                  <  2:1519075/79‑1 (MQ=255)
                                                                                                                           gACATCGTCATTGCTTCAGTATTACGAGACATAGCACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTc                       <  2:1787063/121‑1 (MQ=255)
                                                                                                                           gACATCGTCATTGCTTCAGTATTACGAGACATAGCACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTc                       >  1:1787063/1‑121 (MQ=255)
                                                                                                                            aCATCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTcc    <  2:1215654/139‑1 (MQ=255)
                                                                                                                              aTCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGt  <  2:787597/139‑1 (MQ=255)
                                                                                                                               tCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacga                   >  2:1749806/1‑121 (MQ=255)
                                                                                                                               tCGTCATTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacga                   <  1:1749806/121‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CTCATTACGCTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCAGTTGCTTCAGTATTACGAGACATAGTACCGCCTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGT  >  NC_000913/1763854‑1764119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: