Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,687 T→C 11.5% intergenic (‑301/+26) sufA ← / ← rydB Fe‑S cluster assembly protein/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,6870TC11.5% 178.2 / 10.9 61intergenic (‑301/+26)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (4/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTT  >  NC_000913/1764554‑1764832
                                                                                                                                     |                                                                                                                                                  
aTCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                    <  3:97270/149‑1 (MQ=255)
     aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                               <  1:138599/149‑1 (MQ=255)
                aTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                               >  1:1152812/1‑138 (MQ=255)
                aTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                               <  2:1152812/138‑1 (MQ=255)
                          tCTAAGGAGCATTAGGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCt                                                                                                                                                  <  4:897861/110‑1 (MQ=255)
                          tCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTACTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCt                                                                                                                                                  >  3:897861/1‑110 (MQ=255)
                                          gCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCg                                                                                                                                            >  4:592006/1‑99 (MQ=255)
                                          gCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCg                                                                                                                                            <  3:592006/99‑1 (MQ=255)
                                            tgtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCtgt                                                                                                                                                <  1:564005/94‑1 (MQ=255)
                                            tgtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCtgt                                                                                                                                                >  2:564005/1‑94 (MQ=255)
                                            tgtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                    <  3:461364/105‑1 (MQ=255)
                                            tgtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                    >  4:461364/1‑105 (MQ=255)
                                               gatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTaa                                                                                                                                >  4:1283991/1‑106 (MQ=38)
                                                atgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAg                                                                                                                                      <  4:30814/99‑1 (MQ=255)
                                                atgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAg                                                                                                                                      >  3:30814/1‑99 (MQ=255)
                                                    aaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTa                                                                                >  4:247770/1‑149 (MQ=255)
                                                     aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCtt                                                                                                                                  >  4:312261/1‑98 (MQ=255)
                                                     aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCtt                                                                                                                                  <  3:312261/98‑1 (MQ=255)
                                                               ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                               >  2:1309478/1‑91 (MQ=255)
                                                               ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAg                                                                                                                               <  1:1309478/91‑1 (MQ=255)
                                                                tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                    <  3:247770/149‑1 (MQ=255)
                                                                     tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                                                                           >  4:102881/1‑105 (MQ=255)
                                                                     tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                                                                           <  3:102881/105‑1 (MQ=255)
                                                                                      gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAc                                                                                                                      >  4:94102/1‑77 (MQ=255)
                                                                                      gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTAc                                                                                                                      <  3:94102/77‑1 (MQ=255)
                                                                                         aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                                                                >  3:151594/1‑128 (MQ=255)
                                                                                         aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGCCGa                                                                <  4:151594/128‑1 (MQ=255)
                                                                                           aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGa                                                                        <  1:234239/118‑1 (MQ=255)
                                                                                           aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGa                                                                        >  2:234239/1‑118 (MQ=255)
                                                                                            ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                        >  3:966345/1‑149 (MQ=255)
                                                                                              tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc                                      >  2:653756/1‑149 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTaa                                                                                                                                <  3:1337212/49‑1 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTaa                                                                                                                                >  4:1337212/1‑49 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                            <  3:590049/149‑1 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                            <  1:653756/149‑1 (MQ=255)
                                                                                                             ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                   <  3:1127983/121‑1 (MQ=255)
                                                                                                             ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                   >  4:1127983/1‑121 (MQ=255)
                                                                                                                  aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                     >  1:295839/1‑98 (MQ=255)
                                                                                                                  aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                     <  2:295839/98‑1 (MQ=255)
                                                                                                                    aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                    >  1:573097/1‑97 (MQ=255)
                                                                                                                    aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                    <  2:573097/97‑1 (MQ=255)
                                                                                                                        gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGa                                                           <  4:198664/102‑1 (MQ=255)
                                                                                                                        gAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGa                                                           >  3:198664/1‑102 (MQ=255)
                                                                                                                              ccTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat      >  3:724390/1‑149 (MQ=255)
                                                                                                                               cTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                                                  >  2:941718/1‑88 (MQ=255)
                                                                                                                               cTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                                                  <  1:941718/88‑1 (MQ=255)
                                                                                                                               cTTAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc                                                     <  4:953304/101‑1 (MQ=37)
                                                                                                                               cTTAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATACGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc                                                     >  3:953304/1‑101 (MQ=37)
                                                                                                                                ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                  >  2:1012377/1‑71 (MQ=255)
                                                                                                                                ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                  <  1:1012377/71‑1 (MQ=255)
                                                                                                                                ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGttt                               <  4:348133/122‑1 (MQ=255)
                                                                                                                                ttAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGttt                               >  3:348133/1‑122 (MQ=255)
                                                                                                                                ttAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                                                                            <  4:1181833/77‑1 (MQ=25)
                                                                                                                                ttAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                                                                            >  3:1181833/1‑77 (MQ=25)
                                                                                                                                ttAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                                                      >  3:606378/1‑83 (MQ=25)
                                                                                                                                ttAAGCCTGTGTCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                                                      <  4:606378/83‑1 (MQ=37)
                                                                                                                                  aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACTCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt  <  4:966345/149‑1 (MQ=255)
                                                                                                                                  aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                                                                     >  1:772552/1‑66 (MQ=255)
                                                                                                                                  aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                                                                     <  2:772552/66‑1 (MQ=255)
                                                                                                                                    gTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                        <  4:598450/45‑1 (MQ=255)
                                                                                                                                    gTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                        >  3:598450/1‑45 (MQ=255)
                                                                                                                                     tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                        >  4:896625/1‑44 (MQ=255)
                                                                                                                                     tCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                        <  3:896625/44‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                  
ATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTT  >  NC_000913/1764554‑1764832

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: