Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 10,865 A→G G164G (GGT→GGC
T12T (ACA→ACG
yaaW ←
mbiA →
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,8650AG95.8% 69.2 / NA 24G164G (GGT→GGC
T12T (ACA→ACG
yaaW
mbiA
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/12);  minor base T (0/1);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTC  >  NC_000913/10622‑11106
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
aatTTTCCCCGGGTTCCGTTGTCAGTCCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGGCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAATTGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTg                                                                                                                                                                                                                                            <  2:66500/251‑1 (MQ=255)
 atTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGc                                                                                                                                                                                                                                           >  1:219371/1‑251 (MQ=255)
   tttCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGcc                                                                                                                                                                                                                                          >  2:141061/1‑250 (MQ=255)
      ccccGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCaat                                                                                                                                                                                                                                       >  2:204973/1‑250 (MQ=255)
           gTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCCCAGGGCCTCCCAGACCCGCGCCGCGCAGCAACCCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTg                                                                                                                                                                                                                                 >  2:213044/1‑251 (MQ=255)
                         aCCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGGGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACTCCTTTTGCCAATAATTTGTCTTTCATcagcag                                                                                                                                                                                                                   <  1:140677/251‑1 (MQ=255)
                                                       caggcGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCt                                                                                                                                                                                     >  2:191306/1‑251 (MQ=255)
                                                             aTTTGCAGTACGTCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCTCCCTGGCATc                                                                                                                                                                               <  2:144769/251‑1 (MQ=255)
                                                                tGCAGAACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGAc                                                                                                                                                                            >  2:242669/1‑251 (MQ=255)
                                                                       cGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGc                                                                                                                                                                     >  1:194150/1‑251 (MQ=255)
                                                                                                         gACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATAttc                                                                                                                                    <  2:102853/250‑1 (MQ=255)
                                                                                                            cGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTcc                                                                                                                                >  1:84553/1‑251 (MQ=255)
                                                                                                                                           gggCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTg                                                                                                 <  2:55433/251‑1 (MQ=255)
                                                                                                                                            ggCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGTAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGt                                                                                                <  2:219371/251‑1 (MQ=255)
                                                                                                                                              cctccaagagccgcgcagctCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATCTTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTg                                                                                              >  2:196208/11‑251 (MQ=255)
                                                                                                                                                   cagacccGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCaa                                                                                         <  1:191283/247‑1 (MQ=255)
                                                                                                                                                           cgccgcgcAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTCGCTCAATCTCAAAc                                                                                 >  2:36803/1‑251 (MQ=255)
                                                                                                                                                                             gCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTCGCTCAATCTCAAACGTAGACATCTCTTTGTc                                                                <  1:141061/250‑1 (MQ=255)
                                                                                                                                                                                              gctgcGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCg                                              <  1:59995/251‑1 (MQ=255)
                                                                                                                                                                                                  cGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTCCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGtttt                                          <  1:46231/251‑1 (MQ=255)
                                                                                                                                                                                                          cGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATc                                  <  1:191306/251‑1 (MQ=255)
                                                                                                                                                                                                              aactGCGGGTCAATTGGCTGGAAAGAAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGAAGCGGCAGTAGCTCTTCCAGCTCATACACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCCTATTATGCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATATATTTGTCGGATTTCAACTTCACTCGTTTTGAAACATctcgc                              >  2:229840/4‑247 (MQ=255)
                                                                                                                                                                                                                        caattGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGcc                     <  1:187792/247‑1 (MQ=255)
                                                                                                                                                                                                                                          aTGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTc  <  1:178301/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
AATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTC  >  NC_000913/10622‑11106

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