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

GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  NC_000913/71660‑72099
                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGTATG‑CTGGATCTGCCGGTGGCCCAATTTATTCCGCCGACTA‑TTTACGGCTGCCTGCTGTGG‑CCGCCGTTTTTCTTCTTGCCAGGGATTTTGGCGGGCGCGGCCGTCGGATTTTCTTCCCGGAAT‑CAGAGCGGTGGGGTTTaaag                                                                                                                                                                                                                        >  2:209057/1‑226 (MQ=255)
 cGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGTC‑GCTATTCCTGCCGGAATGCCGAGCGGG‑GGGTTTAAATggc                                                                                                                                                                                                                     >  1:155575/1‑227 (MQ=255)
     gCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAAGA‑TTATCGGCTGCCTGCTGTTG‑CCGCCGTTTTACTTCCTGCCCGGGATTCTGGCGGGGGCGGCGGTC‑CATATTCCTGCCGGAATGCAGAGGGGT‑GGGGGTAAATggggggg                                                                                                                                                                                                                 >  2:293514/1‑223 (MQ=255)
                gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGGTTTTGGCGGGGGCGGCGGCC‑GATATTTCTTCCCGCAAGAGGAGCGGG‑GGGGTTAAAGGGTGGTGGGTGGCAAAAg                                                                                                                                                                                                      >  2:304889/1‑228 (MQ=255)
                  ttGTTCCCCGTTTTGGGGGGTTTTTTTGTTGGCCGGGGGGGTCGGCTGGTGCCATGGGTGGCGGGAATG‑CTGGACCTCCGGGGGGTTGAGTTTATTCCGCCGAATA‑TTACCGGGTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGAGTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCg                                                                                                                                                                                                    <  2:325177/227‑1 (MQ=255)
                              cTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGCTTCTGGCGGGCGCGGGGATC‑GATATTCCTGCCGGAATGCAGGGCGGG‑TGGTTTGAATGGTTGCTGCTGGCAACAGCGGGGGTTTTgggg                                                                                                                                                                                        >  1:104262/1‑228 (MQ=255)
                               tGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTTTGG‑GCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGGGATC‑GCAATTCCTTCCGGCATGCAGAGCCGG‑GCGGTTACATGGTGGCTGCTGGTAACAGCGGGGtgttttgggg                                                                                                                                                                                       >  2:189722/1‑223 (MQ=255)
                                 gtcggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTTGCGGGCGCGGGGGTC‑CGTTATCCTCCCCGACAGCAGCGGGGT‑GGGGTTAAAGGGGTGTCGCGTGCAACACCGGGGttttttgggggg                                                                                                                                                                                     >  2:379966/1‑221 (MQ=255)
                                    ggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGG‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGTTTTTGGGGGTtggt                                                                                                                                                                                  >  1:448242/1‑228 (MQ=255)
                                      tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTTGATCTGCCGGTGGCTAAATTTATTACGCCGGATA‑TTATCGGCTGCCCGCTGTGG‑CCCCCGTTTTTATTCCTGCCCGGGATTCTGGCGGGCCCGGGGATT‑GCTATTTCTGCCCGAATCGCGAACGGG‑GGGTTTTAATGGTTGTTTCATGCACCCACGGGGTTTTTTTGGGTtggggg                                                                                                                                                                                >  2:175273/1‑228 (MQ=255)
                                                  cccGACGGGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCt                                                                                                                                                                     <  1:432211/227‑1 (MQ=255)
                                                              gtgggTGCAAATGGGGGCGGGAAGG‑CTGGATCGGCCGGTGGCTAAATTTATTACGCCGAATA‑TTGTCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                        <  2:488268/225‑1 (MQ=255)
                                                              ggtGGTGCC‑ATGGTGGGGGGTATGTCTGGTCCTGCGGGTGCGAAAATTTTTTA‑GCCCAATATTTATGGGGTGCC‑GTTGTGGCCCGCGGTTTTTCTTCCGGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                        <  2:345358/226‑1 (MQ=255)
                                                                                        ggggtacTGCCGGGGGCTAAATTTTTTGCGCCGAATA‑TTATCGGCTGCCTGCTGGGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                <  1:209057/221‑1 (MQ=255)
                                                                                        cTGGATCTGCCGGTGGCTAAATTTATTACGGCGAATA‑TTAGCGGCTGCCTGGTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGttt                                                                                                                                <  1:452875/227‑1 (MQ=255)
                                                                                                         tatttttATTACGCAGAATT‑GTAGCGGCTGCCTGCTGGGG‑CCGCAGTTTTACTTTCTGCCAGTGGTTCTGGCGGGCGCGGCGATC‑GATACTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGATGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                              <  2:330623/224‑1 (MQ=255)
                                                                                                            aTTTATTACGCCGAATT‑TTACCGGCTGCCGGCGGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCcgcg                                                                                                            <  1:341547/227‑1 (MQ=255)
                                                                                                                 ttaCGCCGAATA‑TTAGCGGCTGCCTGCTGTGG‑CCGGCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGt                                                                                                       <  1:586567/227‑1 (MQ=255)
                                                                                                                                        cctgcTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTTTGGGGCAGCGGGTAAGCGGCGGACCGGTTTAGTCATTATTTTTGCCCCGGGGGGTTGTTGTGGGTGGCGCCCtttgt                                                                                >  1:520821/1‑225 (MQ=255)
                                                                                                                                                           ttACCTCCTGCCGGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  1:201832/228‑1 (MQ=255)
                                                                                                                                                            ttctgccggCCAGGGATTTTGGGGGGCGGGGGTATC‑GATTGTCCTGCCGGAATGTAGAGGGGG‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTGGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  2:380458/219‑1 (MQ=255)
                                                                                                                                                                  cTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtg                                                       <  1:806/228‑1 (MQ=255)
                                                                                                                                                                       aGGGATTTTTGGTGGCGCGGGGGGT‑GATTTTCCGGCCGGAATGCAGGGGGGT‑GAGGTTAAATGGTTGCTGCTGGCAAGAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt                                                  <  2:104262/228‑1 (MQ=255)
                                                                                                                                                                                        cggcgATC‑GATATTGCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt                                 <  1:379966/228‑1 (MQ=255)
                                                                                                                                                                                              tc‑gatATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTTGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTAATTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                            <  1:428020/227‑1 (MQ=255)
                                                                                                                                                                                                       ccTGCGGGAATGGAGGGTGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTTGTTTGGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGta                    <  1:608665/227‑1 (MQ=255)
                                                                                                                                                                                                         tGCCGGAATGCAGGGCGGT‑GATTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGtata                  <  1:304889/227‑1 (MQ=255)
                                                                                                                                                                                                           ccGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGGGGTGTTTTTGTGGGTTGGGTGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCg               <  1:189722/228‑1 (MQ=255)
                                                                                                                                                                                                                aTGGAGAGTGGT‑GGGTTTAAATGGTTGCGTCGGGCAACAGGGGGTTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATAtt          <  2:493559/228‑1 (MQ=255)
                                                                                                                                                                                                                     ggggggT‑GGGTTTAAATGGTTGGTGCTGGCAACAGGGGGGTTTTTGGGGGTTGGTGGCTGGTTGTGCGGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCg     <  2:590194/224‑1 (MQ=255)
                                                                                                                                                                                                                        cGGT‑GAGTTTAAATGGGTGCTGCTGGCAACAGCGGTGGTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaa  <  1:293514/228‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATC‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  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: ‑