Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,291,030 A→T 12.2% L263L (CTA→CTT rssB → PcnB‑degradosome interaction factor; response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,291,0300AT12.2% 78.8 / 5.1 33L263L (CTA→CTTrssBPcnB‑degradosome interaction factor; response regulator
Reads supporting (aligned to +/- strand):  ref base A (16/13);  new base T (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CCTTGTTATTACGCGCATTGTTTAACGGATTATTACAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATGGTGTTCCGTTAGGCACTTTAGGTAACGCTTATTTGAATCAAT  >  NC_000913/1290888‑1291175
                                                                                                                                              |                                                                                                                                                 
ccTTGTTATTACGCGCATTGTTTAACGGATTATTACAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGt                                                                                                                                              <  1:192015/148‑1 (MQ=255)
    gTTATTACGCGCATTGTTTAACGGATTATTACAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTg                                                                                                                                            <  4:243666/146‑1 (MQ=255)
           cgcgCATTGTTTAACGGATTATTACAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCtattat                                                                                                                                    >  2:248350/1‑147 (MQ=255)
            gcgcATTGTTTAACGGATTGTTACAGGAACAGCATGCACACCAAAATCAACGGTTGCCAGAGGTGGGCGCGTTATTGAAGCAGGTAAACCAATTACGTCGTCAGGCCAAACTGCCGGGGCACGTTCCGCTATTAGTTGGCTAATatca                                                                                                                                  >  2:79639/1‑148 (MQ=255)
                  tGTTTAACGGATTATTACAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATcgcg                                                                                                                             >  1:183176/1‑147 (MQ=255)
                              tattaCAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTatcat                                                                                                                                 >  7:119491/1‑131 (MQ=255)
                              tattaCAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTatcat                                                                                                                                 <  8:119491/131‑1 (MQ=255)
                               attaCAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAAtct                                                                                                                <  7:324363/147‑1 (MQ=255)
                                    aGGAACAGCTCGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCt                                                                                                          <  7:250975/148‑1 (MQ=255)
                                       aaCAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGt                                                                                                       >  8:57830/1‑148 (MQ=255)
                                         cAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCtatt                                                                                                                                      >  7:148381/1‑115 (MQ=255)
                                          aGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTc                                                                                                    >  8:165708/1‑148 (MQ=255)
                                              tGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCg                                                                                                >  5:234375/1‑148 (MQ=255)
                                                  cacCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTatc                                                                                                                                   <  2:191716/109‑1 (MQ=255)
                                                  cacCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTatc                                                                                                                                   >  1:191716/1‑109 (MQ=255)
                                                    ccAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCt                                                                                           >  4:169901/1‑147 (MQ=255)
                                                                 ttGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTTTTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTa                                                                              >  5:242554/1‑147 (MQ=255)
                                                                 ttGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTTTTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTa                                                                              >  5:238080/1‑147 (MQ=255)
                                                                 ttGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTaa                                                                             <  2:183176/148‑1 (MQ=255)
                                                                    ccAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAAt                                                                           <  2:78781/147‑1 (MQ=255)
                                                                     cAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATAc                                                                         <  6:234375/148‑1 (MQ=255)
                                                                        agTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTTTTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCgg                                                                      <  6:238080/148‑1 (MQ=255)
                                                                        agTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTTTTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCgg                                                                      <  6:242554/148‑1 (MQ=255)
                                                                                   ttATTGAAGCAGGTAAACCATTTACTTCGTCAGGCTAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGa                                                                                         >  2:70466/1‑118 (MQ=255)
                                                                                   ttATTGAAGCAGGTAAACCATTTACTTCGTCAGGCTAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGa                                                                                         <  1:70466/118‑1 (MQ=255)
                                                                                                tAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTaa                                                                             >  7:150758/1‑117 (MQ=255)
                                                                                                tAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTaa                                                                             <  8:150758/117‑1 (MQ=255)
                                                                                                   aCCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATgg                                           >  8:84133/1‑148 (MQ=255)
                                                                                                        ttACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATGGTGTTc                                      >  4:11035/1‑148 (MQ=255)
                                                                                                                   ggCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATGGTGTTCCGTTAGGCACt                           >  6:132228/1‑148 (MQ=255)
                                                                                                                    gCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATGGTGTTCCGTTAGGCACtt                          <  3:169901/148‑1 (MQ=255)
                                                                                                                     ccAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATGGTGTTCCGTTAGGCACtt                          >  7:72202/1‑147 (MQ=255)
                                                                                                                                            cTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATGGTGTTCCGTTAGGCACTTTAGGTAACGCTTATTTGAATCAAt  <  7:165708/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CCTTGTTATTACGCGCATTGTTTAACGGATTATTACAGGAACAGCTTGCACACCAAAATCAACGGTTGCCAGAGTTGGGCGCGTTATTGAAGCAGGTAAACCATTTACTTCGTCAGGCCAATCTGCCGGGGCAGTTTCCGCTATTAGTTGGCTATTATCATCGCGAACTGAAAAATCTCATTCTGGTTTCTGCGGGTCTGAATGCGACGTTAAATACCGGCGAACACCAGGTGCAAATCAGTAATGGTGTTCCGTTAGGCACTTTAGGTAACGCTTATTTGAATCAAT  >  NC_000913/1290888‑1291175

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: