Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,714,244 T→C 7.1% intergenic (+14/+195) srmB → / ← yfiE ATP‑dependent RNA helicase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,714,2440TC 7.1% 161.1 / 3.3 57intergenic (+14/+195)srmB/yfiEATP‑dependent RNA helicase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (20/32);  new base C (2/2);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAA  >  NC_000913/2714100‑2714389
                                                                                                                                                |                                                                                                                                                 
gTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTaa                                                                                                                                               <  2:201717/149‑1 (MQ=255)
gTGCTCCGAAGAAAAAAGCTAACGAGAACGAGAAGCCGCGGGTGAACAAACGCCATCCCGACACCAAAAATATTGGTAAGCGCAGTAAACCACGCGGACCGGGCGTGCCACCGCAAACGACAGAAGCGTAATCTCAATGCCAGGTTaaa                                                                                                                                               <  4:170456/149‑1 (MQ=255)
     gagaAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTg                                                                                                                                          >  4:86921/1‑149 (MQ=255)
      agaagaAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAACCGCCATAGAGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGTTTTACGCATgg                                                                                                                                         >  2:67585/1‑149 (MQ=255)
      agaagaAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCt                                                                                                                                           <  3:175835/147‑1 (MQ=255)
      agaagaAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTgg                                                                                                                                         <  3:55146/149‑1 (MQ=255)
      agaagaAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTgg                                                                                                                                         <  4:109433/149‑1 (MQ=255)
      agaagaAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTg                                                                                                                                          <  7:45516/148‑1 (MQ=255)
          gaaAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCTGGCAtt                                                                                                                                     <  8:243381/149‑1 (MQ=255)
          gaaAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACC‑AAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCTGGCAtt                                                                                                                                     >  6:118348/1‑148 (MQ=255)
           aaaaaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACC‑AAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCTGGCAttga                                                                                                                                   <  5:118348/149‑3 (MQ=255)
                     aagagaaagagaAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTgg                                                                                                                          <  6:156987/149‑1 (MQ=255)
                     aagagaaagagaAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTg                                                                                                                           <  6:105188/148‑1 (MQ=255)
                        agaaagagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCa                                                                                                                        >  7:109520/1‑148 (MQ=255)
                                     gcgGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGcc                                                                                                          >  8:62877/1‑149 (MQ=255)
                                      cgGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCt                                                                                                         <  2:98921/149‑1 (MQ=255)
                                          tGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATg                                                                                                     >  1:85074/1‑149 (MQ=255)
                                            aaaaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATgc                                                                                                    <  4:102638/148‑1 (MQ=255)
                                             aaaaaaGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAAACAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATgcg                                                                                                   <  6:136321/148‑1 (MQ=255)
                                             aaaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATgcgc                                                                                                  >  4:36049/1‑149 (MQ=255)
                                             aaaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATgcgc                                                                                                  >  6:111549/1‑149 (MQ=255)
                                              aaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTaa                                                                                                                                   <  5:130835/115‑1 (MQ=255)
                                              aaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTaa                                                                                                                                   >  6:130835/1‑115 (MQ=255)
                                              aaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGc                                                                                                                         <  1:206812/125‑1 (MQ=255)
                                               aaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTa                                                                                                <  7:62877/149‑1 (MQ=255)
                                                  cGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGc                                                                                             <  3:104695/149‑1 (MQ=255)
                                                         gcgACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGGTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                       <  2:181557/148‑1 (MQ=255)
                                                           gACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                                    <  1:246803/149‑1 (MQ=255)
                                                                 aaaaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTTAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                              >  6:76696/1‑149 (MQ=255)
                                                                 aaaaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                              >  7:16359/1‑149 (MQ=255)
                                                                   aaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGATGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                            <  2:24801/149‑1 (MQ=255)
                                                                   aaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                            >  6:213720/1‑149 (MQ=255)
                                                                    aaTATTGGTAAGCGCCGAAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTTATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGGGCTGGCATACGCAATACTGCCTTATGCGCTACGCTTATCAGGCCTACGTGGc                                                                           <  4:27170/149‑1 (MQ=39)
                                                                        ttGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTTCCACCGCATACGACAGATTAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCACACGCAAAACTGCCTGATGCGCTTCGCTTATCAGGCCTACGTGGCTCa                                                                        <  6:64700/148‑1 (MQ=38)
                                                                           gTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGc                                                                     >  3:229490/1‑148 (MQ=255)
                                                                                cgcCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatata                                                               <  8:16359/149‑1 (MQ=255)
                                                                                       aaCCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAt                                                         <  5:213720/148‑1 (MQ=255)
                                                                                        aCCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAAttt                                                       <  3:86921/149‑1 (MQ=255)
                                                                                            aGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCAc                                                   <  3:36049/149‑1 (MQ=255)
                                                                                                gAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGAAATATATTGaa                                                          <  6:37009/138‑1 (MQ=255)
                                                                                                gAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGAAATATATTGaa                                                          >  5:37009/1‑138 (MQ=255)
                                                                                                     ggCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTa                                          <  5:76696/149‑1 (MQ=255)
                                                                                                           ccACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                                     >  7:94936/1‑148 (MQ=255)
                                                                                                            cACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGa                                   >  3:228711/1‑149 (MQ=255)
                                                                                                            cACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGa                                   <  2:85074/149‑1 (MQ=255)
                                                                                                              ccGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatat                                                                >  4:81447/1‑118 (MQ=255)
                                                                                                              ccGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatat                                                                <  3:81447/118‑1 (MQ=255)
                                                                                                                  aatcGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGc                             <  4:105823/146‑1 (MQ=255)
                                                                                                                    aCGACAGAAGAGTAATCTCAATGACAGGTTTACGCCTGGTATTAAAAGTGCTGGCAACCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGt                           <  1:67585/149‑1 (MQ=255)
                                                                                                                     cGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGtt                          <  5:111549/149‑1 (MQ=255)
                                                                                                                         agaagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                      >  1:31244/1‑149 (MQ=255)
                                                                                                                           aagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgcc                    >  3:195296/1‑149 (MQ=255)
                                                                                                                           aagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgcc                    >  5:178265/1‑149 (MQ=255)
                                                                                                                                cacttttaaTACCAGGCTTAAACCTGGCATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATc               >  8:8706/8‑149 (MQ=18)
                                                                                                                                      cAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTa                                                                                                <  4:40526/62‑1 (MQ=255)
                                                                                                                                      cAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTa                                                                                                >  3:40526/1‑62 (MQ=255)
                                                                                                                                             aGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATaaa  <  4:228711/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAA  >  NC_000913/2714100‑2714389

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