Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 480,236 C→G 14.6% E25Q (GAG→CAG)  hha ← modulator of gene expression, with H‑NS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913480,2360CG14.6% 120.5 / 14.0 41E25Q (GAG→CAG) hhamodulator of gene expression, with H‑NS
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base G (3/3);  total (21/20)
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.11e-01

ATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAA  >  NC_000913/480097‑480382
                                                                                                                                           |                                                                                                                                                  
aTAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtgt                                                                                                                                            >  1:983742/1‑148 (MQ=255)
                     gggATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTa                                                                                                                       >  1:341255/1‑148 (MQ=255)
                               cGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTaa                                                                                                             <  2:341255/148‑1 (MQ=255)
                                           tCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCg                                                                                                                                               >  2:943482/1‑102 (MQ=255)
                                           tCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCg                                                                                                                                               <  1:943482/102‑1 (MQ=255)
                                              tGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATATCACGCTCCAGCGTGTCAATTGTc                                                                                                                                    <  2:633695/110‑1 (MQ=255)
                                              tGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAGCACGCTCCAGCGTGTCAATTGTc                                                                                                                                    >  1:633695/1‑110 (MQ=255)
                                                gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataat                                                                                                        >  2:952309/1‑136 (MQ=255)
                                                gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataat                                                                                                        <  1:952309/136‑1 (MQ=255)
                                                    aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacga                                                                                                                           <  2:1025346/113‑1 (MQ=255)
                                                    aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacga                                                                                                                           >  1:1025346/1‑113 (MQ=255)
                                                     ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                                                                            <  1:634993/111‑1 (MQ=255)
                                                     ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                                                                            >  2:634993/1‑111 (MQ=255)
                                                                     tCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTc                                                                       <  2:983742/148‑1 (MQ=255)
                                                                                gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataa                                                                                                         >  4:1098937/1‑103 (MQ=255)
                                                                                gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataa                                                                                                         <  3:1098937/103‑1 (MQ=255)
                                                                                  aaaaTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt                                                          >  4:254449/1‑148 (MQ=255)
                                                                                      tACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgt                                                                                                                                            <  1:340574/62‑1 (MQ=255)
                                                                                      tACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgt                                                                                                                                            >  2:340574/1‑62 (MQ=255)
                                                                                        ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg                                                               >  1:52137/1‑137 (MQ=255)
                                                                                        ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg                                                               <  2:52137/137‑1 (MQ=255)
                                                                                                tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTgg                                                                                                                                 >  3:159215/1‑63 (MQ=255)
                                                                                                tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTgg                                                                                                                                 <  4:159215/63‑1 (MQ=255)
                                                                                                tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCgg                                                                                >  2:581418/1‑112 (MQ=255)
                                                                                                tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCgg                                                                                <  1:581418/112‑1 (MQ=255)
                                                                                                       cTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAAc                                     >  2:841870/1‑148 (MQ=255)
                                                                                                              ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                           <  1:701649/87‑1 (MQ=255)
                                                                                                              ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                           >  2:701649/1‑87 (MQ=255)
                                                                                                               tCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGa                             <  3:254449/148‑1 (MQ=255)
                                                                                                                  tatttattTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATtc                          >  4:135500/1‑148 (MQ=255)
                                                                                                                                 aTAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTa                                                <  3:671387/111‑1 (MQ=255)
                                                                                                                                 aTAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTa                                                >  4:671387/1‑111 (MQ=255)
                                                                                                                                  tgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                           <  2:1033114/65‑1 (MQ=255)
                                                                                                                                  tgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                           >  1:1033114/3‑67 (MQ=255)
                                                                                                                                     cacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaata                                                      <  2:380908/101‑1 (MQ=255)
                                                                                                                                     cacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaata                                                      >  1:380908/1‑101 (MQ=255)
                                                                                                                                      acGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGttttt                                                                                   >  2:751641/1‑71 (MQ=255)
                                                                                                                                      acGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGttttt                                                                                   <  1:751641/71‑1 (MQ=255)
                                                                                                                                      acGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACaaa      <  3:135500/148‑1 (MQ=255)
                                                                                                                                       cGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATa                                                                           >  3:290342/1‑78 (MQ=255)
                                                                                                                                       cGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATa                                                                           <  4:290342/78‑1 (MQ=255)
                                                                                                                                           ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTa                                                                     >  2:347309/1‑80 (MQ=255)
                                                                                                                                           ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTa                                                                     <  1:347309/80‑1 (MQ=255)
                                                                                                                                           ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaacaa  >  2:813295/1‑147 (MQ=255)
                                                                                                                                           ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaacaa  <  1:813295/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
ATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAA  >  NC_000913/480097‑480382

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: