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
                                                                                                                                     |                                                                                                                                    
cTCATTACGCTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCtt                                                                                                                                <  1:1232563/139‑1 (MQ=255)
   aTTACGCTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAg                                                                                                                             >  1:1051186/1‑139 (MQ=255)
         cTTCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTAc                                                                                                                       >  1:1549495/1‑139 (MQ=255)
          ttCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACg                                                                                                                      <  1:924658/139‑1 (MQ=255)
          ttCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACg                                                                                                                      >  2:1292723/1‑139 (MQ=255)
          ttCGCCGAAATTGCGCGCACCACCTCTTCGTTTATCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACg                                                                                                                      <  1:1613796/139‑1 (MQ=255)
                                           aTCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGCTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCg                                                                                     >  1:1472103/1‑139 (MQ=255)
                                            tCCCCTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgc                                                                                    <  2:1180855/139‑1 (MQ=255)
                                              cccTTTGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgcca                                                                                  <  2:25567/139‑1 (MQ=255)
                                                  ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                              <  2:1051186/139‑1 (MQ=255)
                                                  ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                              <  2:1472103/139‑1 (MQ=255)
                                                  ttGCCAGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCgccacagc                                                                              <  2:1549495/139‑1 (MQ=255)
                                                       aGCTCATCGGTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCAcat                                                                         <  1:1009103/139‑1 (MQ=255)
                                                               ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCATGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                 <  1:839436/139‑1 (MQ=255)
                                                               ggTGGCTAACTGGGTGAAGAATCCTTCTTTATAATGCAGCGGGCCGCCGGTCCAGGTTTTGTCATCGTCAGTT‑CTTCAGTATTACGAGACATAGTACCGCTTATACCCTAAAGCTTTCGCCACAGCCACTTTCATTCTg                                                                 <  1:502516/139‑1 (MQ=255)
                                                                   gCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACg                                                                                                                      >  1:1928272/1‑82 (MQ=255)
                                                                   gCTAACTGGGTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACg                                                                                                                      <  2:1928272/82‑1 (MQ=255)
                                                                     tAACTGGGTGAATAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACAGCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCttt                                                           <  2:1946130/139‑1 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCGCAAAGCTTTCGCCACAGCCACATTCATTCTg                                                                 >  1:1336599/1‑127 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                              <  1:1691105/98‑1 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCa                                                                                              >  2:1691105/1‑98 (MQ=255)
                                                                           ggTGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACAGTCATTCTg                                                                 <  2:1336599/127‑1 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTa                                                                                                            >  1:1674856/1‑82 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTa                                                                                                            <  2:1674856/82‑1 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCtt                                                            >  1:26386/1‑130 (MQ=255)
                                                                             tGAAGAATCCTTCTTTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCtt                                                            <  2:26386/130‑1 (MQ=255)
                                                                                      cttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCtt                                                                                                                                >  1:841059/1‑53 (MQ=255)
                                                                                      cttcttTATAATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCtt                                                                                                                                <  2:841059/53‑1 (MQ=255)
                                                                                            tataATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACAt                                                                                                                <  1:403532/63‑1 (MQ=255)
                                                                                            tataATTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACAt                                                                                                                >  2:403532/1‑63 (MQ=255)
                                                                                               aaTTCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGAt                                 >  1:787597/1‑139 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTAc                                                                                                           >  1:374849/1‑63 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTAc                                                                                                           <  2:374849/63‑1 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                        >  1:1162071/1‑130 (MQ=255)
                                                                                                 ttCAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAat                                        <  2:1162071/130‑1 (MQ=255)
                                                                                                   cAGCGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAg                             <  1:1292723/139‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACg                                                                                                                      <  2:1424286/47‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACg                                                                                                                      >  1:1424286/1‑47 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTAc                                                                                                           <  1:608457/58‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTAc                                                                                                           >  2:608457/1‑58 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                          <  2:1554280/107‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTa                                                          >  1:1554280/1‑107 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                >  1:170679/1‑117 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGaa                                                <  2:170679/117‑1 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTcc                          >  1:502773/1‑139 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATATGATTAAGTcc                          >  2:1332390/1‑139 (MQ=255)
                                                                                                      cGGGCCGCCGGTCCAGG‑TTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCGCAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCt                         >  1:1215654/1‑139 (MQ=255)
                                                                                                              cgGTCCAGGTTTTGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacgaa                  <  1:1332390/139‑1 (MQ=255)
                                                                                                                          tGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                  >  1:1519075/1‑79 (MQ=255)
                                                                                                                          tGACATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCt                                                                  <  2:1519075/79‑1 (MQ=255)
                                                                                                                           gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGCACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTc                       <  2:1787063/121‑1 (MQ=255)
                                                                                                                           gACATCGTCA‑TTGCTTCAGTATTACGAGACATAGCACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTc                       >  1:1787063/1‑121 (MQ=255)
                                                                                                                            aCATCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTcc    <  2:1215654/139‑1 (MQ=255)
                                                                                                                              aTCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCACGAACGAAATCGACTTCCGt  <  2:787597/139‑1 (MQ=255)
                                                                                                                               tCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacga                   >  2:1749806/1‑121 (MQ=255)
                                                                                                                               tCGTCA‑TTGCTTCAGTATTACGAGACATAGTACCGCTTATACCCCAAAGCTTTCGCCACAGCCACATTCATTCTGGGCTTTAGGGTTGTGAAATTTGAATATCTGATTAAGTCCTTCacga                   <  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: ‑