Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,763,451 C→T 7.7% 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,4510CT 7.7% 163.0 / 3.9 52V187M (GTG→ATG) ‡sufBcomponent of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base T (2/2);  total (27/25)
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.23e-01

TCAAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCC  >  NC_000913/1763303‑1763590
                                                                                                                                                    |                                                                                                                                           
tCAAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGcaccac                                                                                                                                             <  1:207875/149‑1 (MQ=255)
     cTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGc                                                                                                                                        <  3:336273/149‑1 (MQ=255)
       gCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCg                                                                                                                                      <  1:605115/149‑1 (MQ=255)
            gTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGata                                                                                                                                 >  4:265721/1‑149 (MQ=255)
            gTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGata                                                                                                                                 >  3:51529/1‑149 (MQ=255)
              tttttCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATAtt                                                                                                                               >  2:111997/1‑149 (MQ=255)
                    tGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCa                                                                                                                         >  4:128794/1‑149 (MQ=255)
                              cGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCgg                                                                                                               <  3:265721/149‑1 (MQ=255)
                                        tGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTgg                                                                                                     >  2:578688/1‑149 (MQ=255)
                                         ggAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTg                                                                                                      >  1:117584/1‑147 (MQ=255)
                                         ggAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTg                                                                                                      <  2:117584/147‑1 (MQ=255)
                                         ggAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGa                                                                                                    >  2:578796/1‑149 (MQ=255)
                                               ttCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCt                                                                                              >  1:215309/1‑149 (MQ=255)
                                                  cATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCAc                                                                                           <  4:81063/149‑1 (MQ=255)
                                                          aGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGaa                                                                                   >  1:520787/1‑149 (MQ=255)
                                                          aGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGaa                                                                                   <  3:456830/149‑1 (MQ=255)
                                                           gcgcACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                                                                  <  2:520787/149‑1 (MQ=255)
                                                                 gCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaaaa                                                                            >  2:511510/1‑149 (MQ=255)
                                                                   cTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaata                                                                          >  4:607193/1‑149 (MQ=255)
                                                                    tttAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataa                                                                         <  4:96592/149‑1 (MQ=255)
                                                                      tAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGGACAGAAAATAAtt                                                                       <  4:68922/149‑1 (MQ=255)
                                                                       aGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCa                                                                                                                         >  1:151311/1‑98 (MQ=255)
                                                                       aGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCa                                                                                                                         <  2:151311/98‑1 (MQ=255)
                                                                        ggCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                     <  1:261060/149‑1 (MQ=255)
                                                                           acaTAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTg                                                                  >  1:291947/1‑149 (MQ=255)
                                                                           acaTAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTg                                                                  >  1:464559/1‑149 (MQ=255)
                                                                                 ataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGc                                                            >  1:565190/1‑149 (MQ=255)
                                                                                      cGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGtt                                                       <  1:511510/149‑1 (MQ=255)
                                                                                           cATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcg                                                  >  3:67487/1‑149 (MQ=255)
                                                                                            aTCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTcgc                                                 <  2:565190/149‑1 (MQ=255)
                                                                                              cAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGa                                               >  3:67481/1‑149 (MQ=255)
                                                                                               agagGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGAt                                              <  1:578688/149‑1 (MQ=255)
                                                                                                 agGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                            <  2:291947/149‑1 (MQ=255)
                                                                                                      aCCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagtagt                                       <  4:67487/149‑1 (MQ=255)
                                                                                                       ccgccgCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTg                                      >  2:563063/1‑149 (MQ=255)
                                                                                                         gccgcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGc                                    >  2:26135/1‑149 (MQ=255)
                                                                                                                      cgGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAg                       <  3:235433/149‑1 (MQ=255)
                                                                                                                       ggCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGt                      <  3:128794/149‑1 (MQ=255)
                                                                                                                             gaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGaaaa                <  2:215309/149‑1 (MQ=255)
                                                                                                                              aagTTGTCATTCCCCGGCACCACGGTGCCGAGATAGTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAAt               >  4:383751/1‑149 (MQ=255)
                                                                                                                               agTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATg              <  3:129910/149‑1 (MQ=255)
                                                                                                                               agTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATg              <  4:51529/149‑1 (MQ=255)
                                                                                                                                  tGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGACTCGAAAATGGCa           >  2:212390/1‑149 (MQ=255)
                                                                                                                                   gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGaaaa                <  3:227765/143‑1 (MQ=255)
                                                                                                                                   gTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGaaaa                >  4:227765/1‑143 (MQ=255)
                                                                                                                                       ttCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCAc      >  1:242233/1‑149 (MQ=255)
                                                                                                                                        tCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAaataa                                                                         >  1:433117/1‑81 (MQ=255)
                                                                                                                                        tCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTCGATCGCCTCACCCAAGGAACAGAAaagaa                                                                         <  2:433117/81‑1 (MQ=255)
                                                                                                                                        atctcGGCACCGTGGTGCCGGGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGaaaa                >  4:340688/5‑138 (MQ=255)
                                                                                                                                        atctcGGCACCGTGGTGCCGGGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGaaaa                <  3:340688/134‑1 (MQ=255)
                                                                                                                                           tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTccc                                                                    <  1:419128/82‑1 (MQ=255)
                                                                                                                                           tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTccc                                                                    >  2:419128/2‑83 (MQ=255)
                                                                                                                                           ccGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGcc  >  1:521802/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
TCAAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCC  >  NC_000913/1763303‑1763590

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: ‑