Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,225 C→T M44I (ATG→ATA
R132R (CGC→CGT
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_00091311,2250CT80.7% 51.6 / 9.4 31M44I (ATG→ATA
R132R (CGC→CGT
yaaW
mbiA
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (12/13);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCG  >  NC_000913/10979‑11465
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
ccAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGGCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATc                                                                                                                                                                                                                                                 <  1:248831/250‑1 (MQ=255)
    gTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATAAGTAc                                                                                                                                                                                                                                            >  2:75828/1‑251 (MQ=255)
           gCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCTATACTATCGCCACCAAAATGCTGTAATTCTTCAGCAATCAGCTGCCAGTTGCGGCGATGGTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGgag                                                                                                                                                                                                                                     <  1:163722/251‑1 (MQ=255)
            cAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGGTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGATTGTTGGCGATACTCTCGCCACCAAAATGCTGTAATTCTCAAGCAATCAGATGACAGATGCGGCGAATCTGCACGGGATGCCCTACCATCGATTTAAACAGTTCGTTGCGTATCAGTAAGCTGgagg                                                                                                                                                                                                                                    >  2:96702/1‑250 (MQ=255)
              gAAAGTGTTCCACTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAAGTTCCCGTGTCCGCGCAGTTTGTTGGCGATACTTTCGCCACCAAAATGCTGTAAGTCTCCAGCAATCATCTGCCAGTTGCGGCGTTGTTGCTCGGGATGCCCTTCCATCGATTTAAACATTTCGTTGCGTATCAGTTCGCTGGAGAgg                                                                                                                                                                                                                                 <  2:73749/251‑1 (MQ=255)
                aaTTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGGTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCg                                                                                                                                                                                                                               <  1:359698/251‑1 (MQ=255)
                 aTTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCg                                                                                                                                                                                                                               >  1:26147/1‑250 (MQ=255)
                   tGTTCCAGTAACTGTTGCTCACTTTCAAAAGTAGACATCTCTTTGTCGGCTTTCAACTTCAATCGCTTTGAAACATCGAGCAACATGGCCCGATAAAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGACACCAAAATGCTGTAATTCTCCGGCAATCATCTGCCAGTTGCGGCGATGTTGCTAGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCCGTACGATGGAGCGGCGAgg                                                                                                                                                                                                                            >  2:44055/1‑250 (MQ=255)
                    gTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGGCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCTAGtt                                                                                                                                                                                                                           <  1:254103/251‑1 (MQ=255)
                           tAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGCTACAATTTACCGTGTCCGCGCAGTTTGTTGGCGCTACTATCGCCACCAAAATGCTGTACTTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGAGGCCCTTCCATCGATTTAAACAGTTCGTTGTGCATCATTACGCTGGAGAGGCGAGTTTTACCtt                                                                                                                                                                                                                    >  2:110411/1‑251 (MQ=255)
                           tAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCt                                                                                                                                                                                                                                         >  1:107447/1‑231 (MQ=255)
                                             aaaCGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAACCAGTTCGTTGCGTATCAGTACGCTGGGACGGCGAGTTTTGCCTTTTTCATTTTGGGGGAGc                                                                                                                                                                                                  >  2:270080/1‑250 (MQ=255)
                                               aCGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAAt                                                                                                                                                                                               <  1:270080/251‑1 (MQ=255)
                                                         ctctTTGTCGGCTTTCAGCCTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACCGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTTAGCAATCggggggaa                                                                                                                                                                                      >  2:181540/1‑245 (MQ=255)
                                                                 cGGCTTTCAGCTTCAAGCGGTTTGAAACATCGAGCAGAATGGCCCGATACAAGTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCaa                                                                                                                                                                              <  1:184675/250‑1 (MQ=255)
                                                                               aaTCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACa                                                                                                                                                               <  1:64789/251‑1 (MQ=255)
                                                                                 tcgttttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGATTGTTGGCGATAATATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGAGAGCAATCGGGCGAAATTTGCCAACTGTGCCTCataacaa                                                                                                                                                              >  2:99405/5‑244 (MQ=255)
                                                                                 ccGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAAt                                                                                                                                                             <  1:181540/250‑1 (MQ=255)
                                                                                  cgttttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATg                                                                                                                                                            >  1:96716/4‑251 (MQ=255)
                                                                                             aTCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAAt                                                                                                                                                 >  1:63835/1‑251 (MQ=255)
                                                                                                      aaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTg                                                                                                                                        <  1:270683/251‑1 (MQ=255)
                                                                                                         ggCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGaa                                                                                                                                      <  1:316429/250‑1 (MQ=255)
                                                                                                                       gaCCGTGTTCGCGCAGTATGTTGTCGATACTATCGCTACCAAAATGCTGTAATGCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGACTTATGCATTATGGGGGAGCAATCGGGCGAAATTTGACAACTGGTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCATGTAAtt                                                                                                                        <  2:70130/249‑1 (MQ=255)
                                                                                                                        aCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTa                                                                                                                       <  1:44055/250‑1 (MQ=255)
                                                                                                                          cGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAAc                                                                                                                     >  1:203656/1‑250 (MQ=255)
                                                                                                                                         tgttgGCGATACTATCGCCACGAACATGCTGTAATTCGCAAGCAATCAGCTGCCAGATGCGGCGATGATGCTCGGAATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCACTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGATGTGGCAATCGGGCGAAGATTGCAACCTGTTCCTCACTACAATGCTGAAGCAATGACAGCTCTGAATCTTTCAGGAAATTAACATTAATTTTTTTTggt                                                                                                     >  2:64219/1‑250 (MQ=255)
                                                                                                                                                         gCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCACTCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGACGAAAATCCAGATCTGAATCATTCAGGACATTAACATTACTTTTTTGGGGTTTCTATCTTCTGGCGt                                                                                     >  1:236353/1‑251 (MQ=255)
                                                                                                                                                                                           tGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGAAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCatat                                                   >  1:280281/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCTAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCg  <  1:100202/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAAc                                                                                                                     <  2:133374/131‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAAc                                                                                                                     >  1:133374/1‑131 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
CCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTGGCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCG  >  NC_000913/10979‑11465

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