Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,701 G→A 11.1% intergenic (‑315/+12) 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,7010GA11.1% 178.3 / 10.3 62intergenic (‑315/+12)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base A (4/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATC  >  NC_000913/1764554‑1764842
                                                                                                                                                    |                                                                                                                                             
aTCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                                                                              <  3:97270/149‑1 (MQ=255)
     aTTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                                                                         <  1:138599/149‑1 (MQ=255)
                aTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                                                                         >  1:1152812/1‑138 (MQ=255)
                aTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                                                                         <  2:1152812/138‑1 (MQ=255)
                                            tgtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                                                                              <  3:461364/105‑1 (MQ=255)
                                            tgtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                                                                              >  4:461364/1‑105 (MQ=255)
                                               gatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTaa                                                                                                                                          >  4:1283991/1‑106 (MQ=38)
                                                    aaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTa                                                                                          >  4:247770/1‑149 (MQ=255)
                                                     aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCtt                                                                                                                                            <  3:312261/98‑1 (MQ=255)
                                                     aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCtt                                                                                                                                            >  4:312261/1‑98 (MQ=255)
                                                               ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                                                                         <  1:1309478/91‑1 (MQ=255)
                                                               ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                                                                         >  2:1309478/1‑91 (MQ=255)
                                                                tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                              <  3:247770/149‑1 (MQ=255)
                                                                     tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                                                                                     <  3:102881/105‑1 (MQ=255)
                                                                     tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                                                                                     >  4:102881/1‑105 (MQ=255)
                                                                                      gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAc                                                                                                                                <  3:94102/77‑1 (MQ=255)
                                                                                      gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAc                                                                                                                                >  4:94102/1‑77 (MQ=255)
                                                                                         aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                                                                          >  3:151594/1‑128 (MQ=255)
                                                                                         aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGCCGa                                                                          <  4:151594/128‑1 (MQ=255)
                                                                                           aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGa                                                                                  >  2:234239/1‑118 (MQ=255)
                                                                                           aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGa                                                                                  <  1:234239/118‑1 (MQ=255)
                                                                                            ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                                  >  3:966345/1‑149 (MQ=255)
                                                                                              tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc                                                >  2:653756/1‑149 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTaa                                                                                                                                          <  3:1337212/49‑1 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTaa                                                                                                                                          >  4:1337212/1‑49 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                      <  3:590049/149‑1 (MQ=255)
                                                                                                        cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                      <  1:653756/149‑1 (MQ=255)
                                                                                                             ttCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                             <  3:1127983/121‑1 (MQ=255)
                                                                                                             ttCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                             >  4:1127983/1‑121 (MQ=255)
                                                                                                                  aaaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                               <  2:295839/98‑1 (MQ=255)
                                                                                                                  aaaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                               >  1:295839/1‑98 (MQ=255)
                                                                                                                    aTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                              >  1:573097/1‑97 (MQ=255)
                                                                                                                    aTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAggg                                                                              <  2:573097/97‑1 (MQ=255)
                                                                                                                        gAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGa                                                                     >  3:198664/1‑102 (MQ=255)
                                                                                                                        gAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGa                                                                     <  4:198664/102‑1 (MQ=255)
                                                                                                                              ccTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat                >  3:724390/1‑149 (MQ=255)
                                                                                                                               cTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                                                            >  2:941718/1‑88 (MQ=255)
                                                                                                                               cTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                                                                            <  1:941718/88‑1 (MQ=255)
                                                                                                                               cTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc                                                               <  4:953304/101‑1 (MQ=37)
                                                                                                                               cTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATACGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc                                                               >  3:953304/1‑101 (MQ=37)
                                                                                                                                ttAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                            <  1:1012377/71‑1 (MQ=255)
                                                                                                                                ttAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                            >  2:1012377/1‑71 (MQ=255)
                                                                                                                                ttAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGttt                                         >  3:348133/1‑122 (MQ=255)
                                                                                                                                ttAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGttt                                         <  4:348133/122‑1 (MQ=255)
                                                                                                                                ttAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                                                                                      <  4:1181833/77‑1 (MQ=25)
                                                                                                                                ttAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                                                                                      >  3:1181833/1‑77 (MQ=25)
                                                                                                                                ttAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                                                                <  4:606378/83‑1 (MQ=37)
                                                                                                                                ttAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                                                                >  3:606378/1‑83 (MQ=25)
                                                                                                                                  aaGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACTCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt            <  4:966345/149‑1 (MQ=255)
                                                                                                                                  aaGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                                                                               <  2:772552/66‑1 (MQ=255)
                                                                                                                                  aaGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                                                                               >  1:772552/1‑66 (MQ=255)
                                                                                                                                    gTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                                  <  4:598450/45‑1 (MQ=255)
                                                                                                                                    gTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                                  >  3:598450/1‑45 (MQ=255)
                                                                                                                                     tCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                                  >  4:896625/1‑44 (MQ=255)
                                                                                                                                     tCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                                                                                  <  3:896625/44‑1 (MQ=255)
                                                                                                                                             gACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc  <  4:724390/149‑1 (MQ=255)
                                                                                                                                              acacAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTa                       >  4:930677/1‑127 (MQ=255)
                                                                                                                                              acacAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTa                       <  3:930677/127‑1 (MQ=255)
                                                                                                                                                acaGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                             <  1:638667/87‑1 (MQ=255)
                                                                                                                                                acaGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                             >  2:638667/1‑87 (MQ=255)
                                                                                                                                                acaGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatt                  <  3:977394/130‑1 (MQ=255)
                                                                                                                                                acaGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatt                  >  4:977394/1‑130 (MQ=255)
                                                                                                                                                    gCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGta                    >  1:728563/1‑124 (MQ=255)
                                                                                                                                                    gCTTAAGGGTTTCTACCCCATCCGGCGCTAATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGta                    <  2:728563/124‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATC  >  NC_000913/1764554‑1764842

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