Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,763,451 C→T 10.2% V187M (GTG→ATG)  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,4510CT10.2% 113.4 / 3.7 39V187M (GTG→ATG) sufBcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base C (19/16);  new base T (2/2);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCAC  >  NC_000913/1763308‑1763592
                                                                                                                                               |                                                                                                                                             
cTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGc                                                                                                                                          <  2:208199/149‑1 (MQ=255)
cTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGc                                                                                                                                          <  3:431834/149‑1 (MQ=255)
   cccGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCga                                                                                                                                       <  2:37270/149‑1 (MQ=255)
     cGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCga                                                                                                                                       >  1:140535/1‑147 (MQ=255)
     cGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCga                                                                                                                                       <  2:140535/147‑1 (MQ=255)
                   ttAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCGTGGTGCCggg                                                                                                                                      >  1:426338/1‑132 (MQ=255)
                   ttAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCGTGGTGCCggg                                                                                                                                      <  2:426338/134‑3 (MQ=255)
                      aTGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTc                                                                                                                    <  2:277669/149‑1 (MQ=255)
                           aaaaTAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGt                                                                                                               >  2:76624/1‑149 (MQ=255)
                                    ggAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGa                                                                                                      >  1:338992/1‑149 (MQ=255)
                                         gTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGcc                                                                                                 >  3:193708/1‑149 (MQ=255)
                                              aTCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCAcc                                                                                            >  3:143673/1‑149 (MQ=255)
                                                   gCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACTGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAgg                                                                                       >  4:234898/1‑149 (MQ=255)
                                                     aGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCATGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGaa                                                                                     >  1:5919/1‑149 (MQ=255)
                                                     aGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGaa                                                                                     <  4:143673/149‑1 (MQ=255)
                                                          aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaa                                                                                >  1:298441/1‑149 (MQ=255)
                                                          aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaa                                                                                >  1:288928/1‑149 (MQ=255)
                                                               tttAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataa                                                                           <  3:458914/149‑1 (MQ=255)
                                                                           aataaaCGTACCATCAGAAGCTACCGCCGCATTAAGCGCGGAAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTACGGGTGATCGTGGATCGCCTCACCAAAGGACCAGAAAATAATTCCCTGCTCCg                                                               >  2:316780/1‑149 (MQ=255)
                                                                                   tACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCATGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGtttt                                                       <  2:5919/149‑1 (MQ=255)
                                                                                       aTCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                       >  1:167589/1‑129 (MQ=255)
                                                                                       aTCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                       <  2:167589/129‑1 (MQ=255)
                                                                                       aTCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCATGCTCCGCCAGTTTTTcgc                                                   >  4:324339/1‑149 (MQ=255)
                                                                                                tACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagtag                                          <  3:324339/149‑1 (MQ=255)
                                                                                                  ccgccgCATTAAGCGCGGCAAAGAAGTTGTCATACCCCGGCACCACGGTGCCGAGATAGTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTg                                        >  4:455978/1‑149 (MQ=255)
                                                                                                     ccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCa                                     >  1:322984/1‑149 (MQ=255)
                                                                                                      cgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCaa                                    <  1:76624/149‑1 (MQ=255)
                                                                                                        cATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAAcc                                  >  4:353957/1‑149 (MQ=255)
                                                                                                             aGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAAc                             >  2:399544/1‑149 (MQ=255)
                                                                                                                gcgGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGTACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGa                          <  2:288928/149‑1 (MQ=255)
                                                                                                                     aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagt                                            >  3:414158/1‑126 (MQ=255)
                                                                                                                     aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagt                                            <  4:414158/126‑1 (MQ=255)
                                                                                                                     aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGa                     >  2:56866/1‑149 (MQ=255)
                                                                                                                        gaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGaaaa                  <  2:298441/149‑1 (MQ=255)
                                                                                                                                     cccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAAt                 >  1:384637/1‑137 (MQ=255)
                                                                                                                                     cccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAAt                 <  2:384637/137‑1 (MQ=255)
                                                                                                                                      ccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGcc    <  4:193708/149‑1 (MQ=255)
                                                                                                                                      ccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGGTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGcc    >  2:190152/1‑149 (MQ=255)
                                                                                                                                        ggCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCAc  <  1:56866/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCAC  >  NC_000913/1763308‑1763592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑