Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 979,025 A→G E900E (GAA→GAG mukB → chromosome partitioning protein MukB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913979,0250AG97.6% 137.9 / ‑4.5 41E900E (GAA→GAGmukBchromosome partitioning protein MukB
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (18/22);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAACGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGAACAGTTAAAAGAAGATTACGCGTACT  >  NC_000913/978892‑979161
                                                                                                                                     |                                                                                                                                        
aaaaTGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTc                                                                                                                                      <  1:582992/138‑1 (MQ=255)
aaaaTGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTc                                                                                                                                      >  2:582992/1‑138 (MQ=255)
  aaTGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGa                                                                                                                                  >  1:522661/1‑140 (MQ=255)
  aaTGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGa                                                                                                                                  >  1:461574/1‑140 (MQ=255)
    tGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATg                                                                                                                                <  1:243148/140‑1 (MQ=255)
    tGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATg                                                                                                                                <  2:461574/140‑1 (MQ=255)
      aTAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGaa                                                                                                                              <  2:820502/140‑1 (MQ=255)
      aTAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGaa                                                                                                                              >  1:406408/1‑140 (MQ=255)
      aTAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGaa                                                                                                                              <  1:745694/140‑1 (MQ=255)
                 cagcGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGc                                                                                                                   >  1:428413/1‑140 (MQ=255)
                    cGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCcgc                                                                                                                <  1:8518/140‑1 (MQ=255)
                    cGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCcgc                                                                                                                >  2:241069/1‑140 (MQ=255)
                    cGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCcgc                                                                                                                >  2:286724/1‑140 (MQ=255)
                    cGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCcgc                                                                                                                <  2:406408/140‑1 (MQ=255)
                               tGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTc                                                                                                     <  1:504473/140‑1 (MQ=255)
                                gAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCcgc                                                                                                                <  2:682483/128‑1 (MQ=255)
                                gAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCcgc                                                                                                                >  1:682483/1‑128 (MQ=255)
                                gAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGtt                                                                                                      <  2:513415/138‑1 (MQ=255)
                                gAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTACCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGtt                                                                                                      >  1:513415/1‑138 (MQ=255)
                                      gCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGtt                                                                                              >  2:155048/1‑140 (MQ=255)
                                                gCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCaa                                                                                    <  2:522661/140‑1 (MQ=255)
                                                  gTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACt                                                                                  >  1:212973/1‑140 (MQ=255)
                                                       ggCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCa                                                                                     <  1:373291/132‑1 (MQ=255)
                                                       ggCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCa                                                                                     >  2:373291/1‑132 (MQ=255)
                                                         cgcTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGaaa                                                                           >  1:255496/1‑140 (MQ=255)
                                                              aaCCACATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGa                                                                      >  2:293552/1‑140 (MQ=255)
                                                                            gcgTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCgg                                                        <  1:293552/140‑1 (MQ=255)
                                                                             cgTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGt                                                       >  1:634163/1‑140 (MQ=255)
                                                                                      cTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCagag                                              <  1:241069/140‑1 (MQ=255)
                                                                                      cTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCagag                                              <  2:212973/140‑1 (MQ=255)
                                                                                      cTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCagag                                              <  1:286724/140‑1 (MQ=255)
                                                                                      cTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCagag                                              <  2:428413/140‑1 (MQ=255)
                                                                                               gatgaCAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGa                                     <  1:725475/140‑1 (MQ=255)
                                                                                                 tgaCAGCCTGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAAc                                   <  2:634163/140‑1 (MQ=255)
                                                                                                         tGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGaa                           <  1:155048/140‑1 (MQ=255)
                                                                                                         tGGCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGaa                           <  1:852513/140‑1 (MQ=255)
                                                                                                         tGGCGGATCGCGTCGATGAAATCCGCGAACGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGaa                           <  1:329670/140‑1 (MQ=255)
                                                                                                          ggCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGaa                           <  2:600814/139‑1 (MQ=255)
                                                                                                          ggCGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGaa                           >  1:600814/1‑139 (MQ=255)
                                                                                                            cGGATCGCGTCGATGAAATCCGCGAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGAACAg                        >  1:106056/1‑140 (MQ=255)
                                                                                                                                  cgAGCGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGAAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGAACAGTTAAAAGAAGATTACGCGTAct  >  1:270782/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                        
AAAATGATAACCAGCAGCAGCGTATTCAGTTTGAGCAGGCGAAAGAGGGCGTTACGGCGCTGAACCGCATTCTGCCGCGTCTCAACCTGTTGGCTGATGACAGCCTGGCGGATCGCGTCGATGAAATCCGCGAACGTCTGGATGAAGCCCAGGAAGCCGCGCGTTTTGTTCAGCAGTTTGGCAATCAACTGGCGAAACTGGAACCGATTGTTTCGGTATTGCAGAGCGACCCGGAACAGTTCGAACAGTTAAAAGAAGATTACGCGTACT  >  NC_000913/978892‑979161

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