Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 977,870 G→A 9.3% intergenic (+36/‑144) cysE → / → iscR serine acetyltransferase/DNA‑binding transcriptional regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947977,8700GA 9.3% 114.8 / 3.0 43intergenic (+36/‑144)cysE/iscRserine acetyltransferase/DNA‑binding transcriptional regulator IscR
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base A (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

ACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCA  >  NC_002947/977730‑978014
                                                                                                                                            |                                                                                                                                                
aCCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTc                                                                                                                                           <  4:71925/148‑1 (MQ=255)
 ccAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCgt                                                                                                                                                 >  6:42259/1‑141 (MQ=255)
 ccAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCgt                                                                                                                                                 <  5:42259/141‑1 (MQ=255)
  cAAGATGGGCAGTGTCTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACg                                                                                                                                        >  7:13906/1‑149 (MQ=255)
            aGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCt                                                                                                                              >  1:46923/1‑149 (MQ=255)
             gTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTg                                                                                                                             <  3:85029/149‑1 (MQ=255)
             gTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTg                                                                                                                             <  4:119816/149‑1 (MQ=255)
                       tgCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCAATGATCGCGCCTCACTCAAGACATGGAGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATcgc                                                                                                                   >  1:61488/1‑149 (MQ=255)
                         ccAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCt                                                                                                                                            <  8:185375/122‑1 (MQ=255)
                         ccAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCt                                                                                                                                            >  7:185375/1‑122 (MQ=255)
                         ccAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCgcgc                                                                                                                 >  8:105506/1‑149 (MQ=255)
                                  tGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTccc                                                                                                        >  2:86544/1‑149 (MQ=255)
                                    ccGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTccc                                                                                                        <  2:46923/147‑1 (MQ=255)
                                    ccGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCt                                                                                                       <  2:151072/148‑1 (MQ=255)
                                     cGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGa                                                                                                     <  7:113217/149‑1 (MQ=255)
                                       gCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATg                                                                                                   >  5:71766/1‑149 (MQ=255)
                                                       acTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCTGCCTCCCTGATGCAATCCCGACTAAAGc                                                                                   <  8:13906/149‑1 (MQ=255)
                                                                         aaGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa                                                                 >  5:56776/1‑149 (MQ=255)
                                                                            tGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                    <  3:208416/111‑1 (MQ=255)
                                                                            tGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                    >  4:208416/1‑111 (MQ=255)
                                                                                ccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAAt                                                          >  3:318402/1‑149 (MQ=255)
                                                                                ccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGCGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAAt                                                          >  5:277676/1‑149 (MQ=255)
                                                                                     gcgACACTCAGCCGGATTGCTCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCg                                                     <  7:92212/149‑1 (MQ=255)
                                                                                       gACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCATGTGCTCTCGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTGTAATAGTTGTCGTAAAAGCTCggg                                                   >  8:206524/1‑149 (MQ=255)
                                                                                       gACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCATGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCgg                                                    >  8:206526/1‑148 (MQ=255)
                                                                                           cTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATc                                               <  4:163746/149‑1 (MQ=255)
                                                                                            tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCCCGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCg                                              >  6:151369/1‑149 (MQ=255)
                                                                                            tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGt                                                                    <  8:242330/127‑1 (MQ=255)
                                                                                            tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGt                                                                    >  7:242330/1‑127 (MQ=255)
                                                                                                ccGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATa                                          >  4:259768/1‑149 (MQ=255)
                                                                                                ccGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATa                                          >  5:100281/1‑149 (MQ=255)
                                                                                                    aTTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCg                                      <  3:220225/149‑1 (MQ=255)
                                                                                                            cgccgcACTGAAGTCATGGTGGCATGCGGGGCATGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATc                              <  7:206526/149‑1 (MQ=255)
                                                                                                            cgccgcACTGAAGTCATGGTGGCATGCCGGGCATGTGCTCACGCCCCTCGTCTGCTACAATCCCGCCCCCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGcacacacc                              <  7:206524/149‑3 (MQ=255)
                                                                                                                     gAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTc                     >  7:99952/1‑149 (MQ=255)
                                                                                                                      aaGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGca                                    <  7:73238/133‑1 (MQ=255)
                                                                                                                      aaGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTcc                    >  8:87372/1‑149 (MQ=255)
                                                                                                                      aaGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTcc                    >  3:254854/1‑149 (MQ=255)
                                                                                                                      aaGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGAATTAAATGCTCGGGATTCGCATACTCgaa                                    >  8:73238/1‑131 (MQ=255)
                                                                                                                         tCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGt                 <  7:87372/149‑1 (MQ=255)
                                                                                                                           aTGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGa                          <  3:99760/138‑1 (MQ=255)
                                                                                                                           aTGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGa                          >  4:99760/1‑138 (MQ=255)
                                                                                                                              gtggCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTa             <  7:105506/148‑1 (MQ=255)
                                                                                                                                  cATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCTGTAGGTCTAATAGTTGACTCAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCaa        <  1:86544/149‑1 (MQ=255)
                                                                                                                                        gggCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCa  <  6:277676/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
ACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCA  >  NC_002947/977730‑978014

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑