Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,314,801 G→T 9.4% Q438H (CAG→CAT) ‡ rcsD → RcsD phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,314,8010GT 9.4% 92.6 / 2.2 42Q438H (CAG→CAT) ‡rcsDRcsD phosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base T (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

ATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTA  >  NC_000913/2314659‑2314943
                                                                                                                                              |                                                                                                                                              
aTTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAg                                                                                                                                                >  1:11716/1‑143 (MQ=255)
    gactCTGGAAATCATCACCTCCATGGCGGAACAGCATGAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATc                                                                                                                                            <  2:28906/140‑1 (MQ=255)
      aTCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcgc                                                                                                                                          >  1:329874/1‑143 (MQ=255)
       tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGCCGCTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCGGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGCTAGGGATcgcg                                                                                                                                         >  1:480028/1‑143 (MQ=255)
         tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGATATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaa                                                                                                                                       <  2:290295/143‑1 (MQ=255)
         tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGCTCAGGAGCGCGaa                                                                                                                                       >  2:679081/1‑143 (MQ=255)
         tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                                       >  1:620861/1‑143 (MQ=255)
         ggCAAAACATCACCACCAGGGGGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                                       <  1:731163/142‑1 (MQ=255)
          gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAAGTTTCACTATTCGCGATCAGCACCGCGAAg                                                                                                                                      >  2:812850/1‑143 (MQ=255)
             aaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCCGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaataat                                                                                                                                   <  1:420456/143‑5 (MQ=255)
               aCATCACCACCATGGCGGAACAGCATCAGGGGAATATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGGACACAAATTTTCATTATTCGCGATCAGGATCGGGGAGTGCTg                                                                                                                                 >  1:235213/1‑143 (MQ=255)
                 aTCCCCACCATGCCGGAACCGCATCAGGGGGCTATTCAGGCGACGATCAATAACGAGCTGTCTGAAATCCGCATGTTCCGCAGCCAGGCCGCGCCGAGCACACAAACTTGCATTATTCGCGATCAGGATCGCGCAGTGCTggg                                                                                                                               >  2:647214/1‑142 (MQ=255)
                   caccacCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaa                                                                                                                             >  1:356006/1‑143 (MQ=255)
                     ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAAc                                                                                                                           >  2:442286/1‑143 (MQ=255)
                     ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCAGCAGCCAGGTCGCGCCGCGCACACAAATTGTCATTATTCGCGCACGGGCTCGCGAAgtggtgacaaag                                                                                                                           >  2:237928/1‑137 (MQ=255)
                       accaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTCTTCGCGATCAGGATCGCGAAGTGCTGGTAAACaa                                                                                                                         >  2:245443/1‑143 (MQ=255)
                           tGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                                     >  1:529205/1‑143 (MQ=255)
                           tGGAGGAACAGCATCAGGGGATTATTCATGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                                     <  2:63668/143‑1 (MQ=255)
                               ggAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCa                                                                                                                 <  1:395313/143‑1 (MQ=255)
                                       aTCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGACACTCAAGCAGgcg                                                                                                         >  1:598250/1‑143 (MQ=255)
                                        tCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCGGATCGCGAAGTGCTGGTAAACAAGAAACTCTAGCAGgcgc                                                                                                        >  1:381761/1‑143 (MQ=255)
                                                     cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTTTga                                                                                           <  1:696977/143‑1 (MQ=255)
                                                     cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATga                                                                                           <  1:442286/143‑1 (MQ=255)
                                                           acgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGACGCGCCGCGCACACAAATTTTCATTCTTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGGATGACaaaaa                                                                                     >  2:570636/1‑143 (MQ=255)
                                                             gaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAcc                                                                                   <  1:748751/143‑1 (MQ=255)
                                                                           tGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatg                                                                     <  2:784168/143‑1 (MQ=255)
                                                                              aTGATATCCGCATGTTCCGCAGCCAGGCCGCGCCGCGCACACAAATTTGCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAATCTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                  <  2:284086/143‑1 (MQ=255)
                                                                              aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                  <  2:268132/143‑1 (MQ=255)
                                                                              aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                  <  1:621478/143‑1 (MQ=255)
                                                                               tGAGATCCGCATGTTCCGCAGCCAGGTAGCGCGGCGCACACAAATTTTCATTATTCGCGCTCAGGATCGCGAAGTGCAGGTAAACAAGAAACTCAAGCAGGCGCAGCGCCTGCATGAGAAAAACAAGCCGCGGCGGATGATcc                                                                 >  2:149121/1‑142 (MQ=255)
                                                                                  gctcggCATGTTCCGCAGCCAGGTCGCGGCGCGCACACAAATTTGCATTATTCGCGATCAGGAGCGCGAAGTGCTGGTAAACAAGAAGCTCAGGCAGGGGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                              <  2:185075/138‑1 (MQ=255)
                                                                                  gctCCGCATTGTGCGCAGCCGGGTCGCGCCGCGCACACAAACTATCCTGCTTGGCGTGCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                              <  2:757745/141‑1 (MQ=255)
                                                                                                       ggTCGCGCCGCGCACACAAATTGTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCCGATGATCTTTAGGAAAAACATTGGCGATgcgc                                         <  2:430021/143‑1 (MQ=255)
                                                                                                        gTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTGAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGgcgctgcgct                                        >  1:573993/1‑143 (MQ=255)
                                                                                                                   cacaaaAATGTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCg                             <  1:625448/143‑1 (MQ=255)
                                                                                                                        aaaTTTTCATTATTCGCGATCAGGATGGCGAAGTGCTGGTAAACAAGAGACGCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGCTGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAg                        <  2:527137/143‑1 (MQ=255)
                                                                                                                              tCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAGAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTg                  >  1:302951/1‑143 (MQ=255)
                                                                                                                                attattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGc                >  2:124494/1‑143 (MQ=255)
                                                                                                                                 cacttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCg               >  1:479431/4‑143 (MQ=255)
                                                                                                                                  tattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCACAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTAGCGATGCGCTGAAAGAAACCGCACAGTCCCTGCCgg              >  2:295563/1‑143 (MQ=255)
                                                                                                                                   atgcgcgATCAGGATCGCGAAGGGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGga             <  1:516972/140‑1 (MQ=255)
                                                                                                                                    ttCGCGATCCTGATCGCGAAGTGCTGGTAAACGAGTAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAGACATTGGCGATGCGCTGAAAGAACCCGCACAGTCACAGGCGgag            >  2:626435/1‑143 (MQ=255)
                                                                                                                                     tCGCGATTAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCCGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGACCCCGCACAGTCCCTGGCGgaga           >  1:326650/1‑143 (MQ=255)
                                                                                                                                              ggATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTa  >  2:120968/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTA  >  NC_000913/2314659‑2314943

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