Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 977,876 T→A 8.0% intergenic (+39/‑141) 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,8760TA 8.0% 145.5 / 3.6 50intergenic (+39/‑141)cysE/iscRserine acetyltransferase/DNA‑binding transcriptional regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base A (2/2);  total (25/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.58e-01

CGCTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCGACTG  >  NC_002947/977728‑978025
                                                                                                                                                    |                                                                                                                                                      
cccggACCAAGATGGACAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGtgtg                                                                                                                                                       <  1:857760/146‑1 (MQ=255)
 gcTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGc                                                                                                                                                      <  2:331918/150‑1 (MQ=255)
   tGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGcac                                                                                                                                                    >  1:696016/1‑148 (MQ=255)
     aCCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCAc                                                                                                                                                  >  1:508176/1‑150 (MQ=255)
       cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCaagg                                                                                                                                                >  2:824266/1‑147 (MQ=255)
       cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGc                                                                                                                                                >  2:533966/1‑150 (MQ=255)
       cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGc                                                                                                                                                <  2:148787/150‑1 (MQ=255)
          gATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGcccc                                                                                                                                             >  1:1065013/1‑150 (MQ=255)
              ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACGccccgcat                                                                                                                                         <  2:696016/150‑5 (MQ=255)
              ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGAGGGGCGTGAGCACACGccccgcat                                                                                                                                         <  2:936881/150‑5 (MQ=255)
                  gTGACTATTGTGCCAAGGAACGGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGc                                                                                                                                     <  1:718587/150‑1 (MQ=255)
                        aTTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAAt                                                                                                                               >  2:84563/1‑150 (MQ=255)
                            tgCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGc                                                                                                                                                      >  2:59817/1‑123 (MQ=255)
                            tgCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGc                                                                                                                                                      <  1:59817/123‑1 (MQ=255)
                              ccAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcc                                                                                                                         >  2:739642/1‑150 (MQ=255)
                               cAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgccc                                                                                                                        <  2:149722/150‑1 (MQ=255)
                               cAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgccc                                                                                                                        <  2:809768/150‑1 (MQ=255)
                                aaGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGtgtg                                                                                                                                                       >  1:470345/1‑118 (MQ=255)
                                aaGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGtgtg                                                                                                                                                       <  2:470345/118‑1 (MQ=255)
                                         ccGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGa                                                                                                              <  1:739642/150‑1 (MQ=255)
                                         ccGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACCACACTTCTGAAGGCATGTGGGGCGTTTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGa                                                                                                              >  1:152068/1‑150 (MQ=255)
                                          cGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                             <  1:767748/150‑1 (MQ=255)
                                                    ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCt                                                                                                                                            <  2:445732/109‑1 (MQ=255)
                                                    ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCt                                                                                                                                            >  1:445732/1‑109 (MQ=255)
                                                    ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGa                                                                                                   <  2:508176/150‑1 (MQ=255)
                                                     gaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGAc                                                                                                  <  1:533966/150‑1 (MQ=255)
                                                           gacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGGTCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATc                                                                                                                              <  2:823348/116‑1 (MQ=255)
                                                           gacTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATc                                                                                                                              >  1:823348/1‑116 (MQ=255)
                                                                    gAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGc                                                                                                                                                >  2:154613/1‑89 (MQ=255)
                                                                    gAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGc                                                                                                                                                <  1:154613/89‑1 (MQ=255)
                                                                           cccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGtt                                                                           >  1:658156/1‑150 (MQ=255)
                                                                            ccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                             >  1:957353/1‑116 (MQ=255)
                                                                            ccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                             <  2:957353/116‑1 (MQ=255)
                                                                             cAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa                                                                         <  1:84563/150‑1 (MQ=255)
                                                                                           ttactcTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCggg                                                           <  2:1235310/145‑1 (MQ=255)
                                                                                                     ccGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTaa                                                                                <  2:1174928/119‑1 (MQ=255)
                                                                                                     ccGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTaa                                                                                >  1:1174928/1‑119 (MQ=255)
                                                                                                        cATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGc                                                     <  2:688782/143‑1 (MQ=255)
                                                                                                        cATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGc                                                     >  1:688782/1‑143 (MQ=255)
                                                                                                         aTTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAAt                                                                                                                               <  1:418097/69‑1 (MQ=255)
                                                                                                         aTTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAAt                                                                                                                               >  2:418097/1‑69 (MQ=255)
                                                                                                          ttGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGca                                            <  2:1065013/150‑1 (MQ=255)
                                                                                                              tCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCa                                                    <  2:409324/138‑1 (MQ=255)
                                                                                                              tCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCa                                                    >  1:409324/1‑138 (MQ=255)
                                                                                                              tCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGcacaca                                        >  2:975532/1‑150 (MQ=255)
                                                                                                                  gccgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTTCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCt                                    >  1:669814/1‑150 (MQ=255)
                                                                                                                    cgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGa                                  <  2:152068/150‑1 (MQ=255)
                                                                                                                     gcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGaa                                 >  2:617814/1‑150 (MQ=255)
                                                                                                                         tGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTc                             >  2:453300/1‑150 (MQ=255)
                                                                                                                              tCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTg                        >  1:533274/1‑150 (MQ=255)
                                                                                                                                aTGGTGGCATGCGGGGCGTGTGCTCACGCCCATCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCACGGCTACACCAGATTGATTCCAAGGATGAACAAACGCGCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGt                      >  1:151421/1‑150 (MQ=255)
                                                                                                                                       gacgaggggCGTGAGCACACGCCCCGCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGATCGGGACTCGAATAAAGGCAAACAGCCGCAAACGCCCGTGGTACCCAAt               >  2:781364/6‑150 (MQ=255)
                                                                                                                                         tGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAg             <  1:453300/150‑1 (MQ=255)
                                                                                                                                                  tgtgCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCgac    >  1:51643/1‑150 (MQ=255)
                                                                                                                                                  tgtgCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCga     >  2:1099656/1‑149 (MQ=255)
                                                                                                                                                  tgtgCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCga     <  1:1099656/149‑1 (MQ=255)
                                                                                                                                                    tgCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCgactg  <  2:51643/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGCTGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCCGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGCCATGCGACTG  >  NC_002947/977728‑978025

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