Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 104,352 A→G E124G (GAG→GGG)  ftsA → ATP‑binding cell division FtsK recruitment protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913104,3520AG97.7% 151.4 / NA 44E124G (GAG→GGG) ftsAATP‑binding cell division FtsK recruitment protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/21);  minor base C (0/1);  total (22/22)
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

ATCTGGCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATAT  >  NC_000913/104220‑104481
                                                                                                                                    |                                                                                                                                 
atCTGGCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCg                                                                                                                             <  2:434812/139‑1 (MQ=255)
  cTGGCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCgtg                                                                                                                           <  2:370303/139‑1 (MQ=255)
    ggCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  2:193630/1‑139 (MQ=255)
     gcgcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGtgc                                                                                                                        >  1:39605/1‑139 (MQ=255)
      cgcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGtgct                                                                                                                       <  1:282729/139‑1 (MQ=255)
            cTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 <  2:39605/139‑1 (MQ=255)
                  aGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         >  1:163861/1‑125 (MQ=255)
                  aGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                         <  2:163861/125‑1 (MQ=255)
                    ctcaTCAGCTGCCAGAATGAAATTGGTATGTTGCCTATTTCTGAAGAAGAAGTCACGCAATATGATGTGGAAAACGTCGTCCATACCGCGAAATCGCTCCGTGTGCGCTATGCGCATCGTGTGCTGCATGTGATCCCGc                                                                                                         <  2:2521/137‑1 (MQ=255)
                       aTCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAag                                                                                                      >  1:20024/1‑139 (MQ=255)
                         cAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAagag                                                                                                    >  2:97014/1‑139 (MQ=255)
                         cAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAagag                                                                                                    >  2:495979/1‑139 (MQ=255)
                                 aGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGAt                                                                                            <  2:246687/139‑1 (MQ=255)
                                  gAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGAtt                                                                                           <  1:97014/139‑1 (MQ=255)
                                  gAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGAtt                                                                                           <  1:193630/139‑1 (MQ=255)
                                         aTTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGGATCGTGTGCTGCATGCGATCCCGCAAGAGTATGCGATTGACTATc                                                                                    >  2:496710/1‑139 (MQ=255)
                                              tATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                               <  1:231669/139‑1 (MQ=255)
                                                 ggTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAggg                                                                            >  2:718321/1‑139 (MQ=255)
                                                 ggTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGA‑‑TGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGAt                                                                          >  1:557213/1‑139 (MQ=255)
                                                        aTTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATc                                                                                                                              >  2:534468/1‑82 (MQ=255)
                                                        aTTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATc                                                                                                                              <  1:534468/82‑1 (MQ=255)
                                                            cTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 >  1:120217/1‑91 (MQ=255)
                                                            cTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 >  1:192529/1‑91 (MQ=255)
                                                            cTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 <  2:120217/91‑1 (MQ=255)
                                                            cTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                 <  2:192529/91‑1 (MQ=255)
                                                               aagaagaagTGACACTAGAAGATGTGGAAAACGTCAACCATACCGCGAAAACGGTGCGTGTGCGCGACGGGCATCGTGTGCTGCACGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGCATCAgga                                                              >  2:344994/1‑138 (MQ=255)
                                                                    gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                               >  1:463213/1‑85 (MQ=255)
                                                                    gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                               <  2:463213/85‑1 (MQ=255)
                                                                                    aTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTc                                                     >  1:698096/1‑127 (MQ=255)
                                                                                    aTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTc                                                     <  2:698096/127‑1 (MQ=255)
                                                                                      gtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGc                                       <  1:496710/139‑1 (MQ=255)
                                                                                                 cgtcCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTg                            >  2:769560/1‑139 (MQ=255)
                                                                                                  gtcCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGc                           >  1:287460/1‑139 (MQ=255)
                                                                                                    ccATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAc                         <  2:20024/139‑1 (MQ=255)
                                                                                                    ccATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCAc                         <  1:356044/139‑1 (MQ=255)
                                                                                                        aCCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTa                                                             <  1:563057/99‑1 (MQ=255)
                                                                                                        aCCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTa                                                             >  2:563057/1‑99 (MQ=255)
                                                                                                        aCCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGa                     <  2:484188/139‑1 (MQ=255)
                                                                                                         ccGCGAAATCGGTGCGTGTGCGCGTTGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGAt                    >  2:522283/1‑139 (MQ=255)
                                                                                                         ccGCGAAATCGGTGCGTGTGCGCGATGGGCATCgtg                                                                                                                           <  2:388913/36‑1 (MQ=255)
                                                                                                         ccGCGAAATCGGTGCGTGTGCGCGATGGGCATCgtg                                                                                                                           >  1:388913/1‑36 (MQ=255)
                                                                                                                    gTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                               <  1:430221/69‑1 (MQ=255)
                                                                                                                    gTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                               >  2:430221/1‑69 (MQ=255)
                                                                                                                           tgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGatat  >  2:608072/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                 
ATCTGGCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATAT  >  NC_000913/104220‑104481

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