Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,877 A→G 27.2% E176G (GAG→GGG)  yabI → DedA family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,8770AG27.2% 47.0 / 3.5 31E176G (GAG→GGG) yabIDedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/19);  new base G (7/2);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 5.42e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  NC_000913/71660‑72099
                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGTATGCTGGATCTGCCGGTGGCCCAATTTATTCCGCCGACTATTTACGGCTGCCTGCTGTGGCCGCCGTTTTTCTTCTTGCCAGGGATTTTGGCGGGCGCGGCCGTCGGATTTTCTTCCCGGAATCAGAGCGGTGGGGTTTaaag                                                                                                                                                                                                                        >  2:209057/1‑226 (MQ=255)
 cGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTCGCTATTCCTGCCGGAATGCCGAGCGGGGGGTTTAAATggc                                                                                                                                                                                                                     >  1:155575/1‑227 (MQ=255)
     gCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAAGATTATCGGCTGCCTGCTGTTGCCGCCGTTTTACTTCCTGCCCGGGATTCTGGCGGGGGCGGCGGTCCATATTCCTGCCGGAATGCAGAGGGGTGGGGGTAAATggggggg                                                                                                                                                                                                                 >  2:293514/1‑223 (MQ=255)
                gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGGTTTTGGCGGGGGCGGCGGCCGATATTTCTTCCCGCAAGAGGAGCGGGGGGGTTAAAGGGTGGTGGGTGGCAAAAg                                                                                                                                                                                                      >  2:304889/1‑228 (MQ=255)
                  ttGTTCCCCGTTTTGGGGGGTTTTTTTGTTGGCCGGGGGGGTCGGCTGGTGCCATGGGTGGCGGGAATGCTGGACCTCCGGGGGGTTGAGTTTATTCCGCCGAATATTACCGGGTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGAGTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCg                                                                                                                                                                                                    <  2:325177/227‑1 (MQ=255)
                              cTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGCTTCTGGCGGGCGCGGGGATCGATATTCCTGCCGGAATGCAGGGCGGGTGGTTTGAATGGTTGCTGCTGGCAACAGCGGGGGTTTTgggg                                                                                                                                                                                        >  1:104262/1‑228 (MQ=255)
                               tGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTTTGGGCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGGGATCGCAATTCCTTCCGGCATGCAGAGCCGGGCGGTTACATGGTGGCTGCTGGTAACAGCGGGGtgttttgggg                                                                                                                                                                                       >  2:189722/1‑223 (MQ=255)
                                 gtcggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTTGCGGGCGCGGGGGTCCGTTATCCTCCCCGACAGCAGCGGGGTGGGGTTAAAGGGGTGTCGCGTGCAACACCGGGGttttttgggggg                                                                                                                                                                                     >  2:379966/1‑221 (MQ=255)
                                    ggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGTTTTTGGGGGTtggt                                                                                                                                                                                  >  1:448242/1‑228 (MQ=255)
                                      tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTTGATCTGCCGGTGGCTAAATTTATTACGCCGGATATTATCGGCTGCCCGCTGTGGCCCCCGTTTTTATTCCTGCCCGGGATTCTGGCGGGCCCGGGGATTGCTATTTCTGCCCGAATCGCGAACGGGGGGTTTTAATGGTTGTTTCATGCACCCACGGGGTTTTTTTGGGTtggggg                                                                                                                                                                                >  2:175273/1‑228 (MQ=255)
                                                  cccGACGGGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCt                                                                                                                                                                     <  1:432211/227‑1 (MQ=255)
                                                              gtgggTGCAAATGGGGGCGGGAAGGCTGGATCGGCCGGTGGCTAAATTTATTACGCCGAATATTGTCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                        <  2:488268/225‑1 (MQ=255)
                                                              ggtGGTGCCATGGTGGGGGGTATGTCTGGTCCTGCGGGTGCGAAAATTTTTTAGCCCAATATTTATGGGGTGCCGTTGTGGCCCGCGGTTTTTCTTCCGGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                        <  2:345358/226‑1 (MQ=255)
                                                                                        ggggtacTGCCGGGGGCTAAATTTTTTGCGCCGAATATTATCGGCTGCCTGCTGGGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                <  1:209057/221‑1 (MQ=255)
                                                                                        cTGGATCTGCCGGTGGCTAAATTTATTACGGCGAATATTAGCGGCTGCCTGGTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                <  1:452875/227‑1 (MQ=255)
                                                                                                         tatttttATTACGCAGAATTGTAGCGGCTGCCTGCTGGGGCCGCAGTTTTACTTTCTGCCAGTGGTTCTGGCGGGCGCGGCGATCGATACTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGATGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                              <  2:330623/224‑1 (MQ=255)
                                                                                                            aTTTATTACGCCGAATTTTACCGGCTGCCGGCGGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcgcg                                                                                                            <  1:341547/227‑1 (MQ=255)
                                                                                                                 ttaCGCCGAATATTAGCGGCTGCCTGCTGTGGCCGGCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGt                                                                                                       <  1:586567/227‑1 (MQ=255)
                                                                                                                                        cctgcTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTTTGGGGCAGCGGGTAAGCGGCGGACCGGTTTAGTCATTATTTTTGCCCCGGGGGGTTGTTGTGGGTGGCGCCCtttgt                                                                                >  1:520821/1‑225 (MQ=255)
                                                                                                                                                           ttACCTCCTGCCGGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  1:201832/228‑1 (MQ=255)
                                                                                                                                                            ttctgccggCCAGGGATTTTGGGGGGCGGGGGTATCGATTGTCCTGCCGGAATGTAGAGGGGGGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTGGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  2:380458/219‑1 (MQ=255)
                                                                                                                                                                  cTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtg                                                       <  1:806/228‑1 (MQ=255)
                                                                                                                                                                       aGGGATTTTTGGTGGCGCGGGGGGTGATTTTCCGGCCGGAATGCAGGGGGGTGAGGTTAAATGGTTGCTGCTGGCAAGAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt                                                  <  2:104262/228‑1 (MQ=255)
                                                                                                                                                                                        cggcgATCGATATTGCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt                                 <  1:379966/228‑1 (MQ=255)
                                                                                                                                                                                              tcgatATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTTGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTAATTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                            <  1:428020/227‑1 (MQ=255)
                                                                                                                                                                                                       ccTGCGGGAATGGAGGGTGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTTGTTTGGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGta                    <  1:608665/227‑1 (MQ=255)
                                                                                                                                                                                                         tGCCGGAATGCAGGGCGGTGATTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGtata                  <  1:304889/227‑1 (MQ=255)
                                                                                                                                                                                                           ccGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGGGGTGTTTTTGTGGGTTGGGTGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCg               <  1:189722/228‑1 (MQ=255)
                                                                                                                                                                                                                aTGGAGAGTGGTGGGTTTAAATGGTTGCGTCGGGCAACAGGGGGTTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATAtt          <  2:493559/228‑1 (MQ=255)
                                                                                                                                                                                                                     ggggggTGGGTTTAAATGGTTGGTGCTGGCAACAGGGGGGTTTTTGGGGGTTGGTGGCTGGTTGTGCGGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCg     <  2:590194/224‑1 (MQ=255)
                                                                                                                                                                                                                        cGGTGAGTTTAAATGGGTGCTGCTGGCAACAGCGGTGGTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaa  <  1:293514/228‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  NC_000913/71660‑72099

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