Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,598 A→G V121A (GTC→GCC)  yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,5980AG96.9% 94.0 / ‑5.7 32V121A (GTC→GCC) yaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/16);  total (15/17)
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

AGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGG‑TAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCT  >  NC_000913/7370‑7812
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                       
aGACCGCGAGTGATCGCCAGCAGAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACTGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAgctg                                                                                                                                                                                                   <  1:347145/251‑4 (MQ=255)
           gatcgCCAGCAGAGCAAAGACAGCGAGAATAATACCTGGCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGcgagggaa                                                                                                                                                                                         <  1:239638/246‑8 (MQ=255)
                 cagcagAGCAAAGACAGCGAGAATAATACCTGTCCCCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCgg                                                                                                                                                                                  <  2:331976/251‑1 (MQ=255)
                  agcagAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGCAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGc                                                                                                                                                                                 >  2:92594/1‑251 (MQ=255)
                    cagAGCAAAGACAGCGAGAATAATACCTGTCACCAGCGGGGGCAAATCAAAAGAAACACTCAGGGCGCGCGCAACGGCGTTCGCCTGAATTCCGCTGCAAATTATGCCGTAGGCGATGCGCAAAAAGCCGACGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCGTCTACCATGCCGGTCCGCCACGCAACTGCCCATTGGCATCACGCTCTTTCTaaatctgtgcgaggaaacattctgac                                                                                                                                                                                >  2:211637/1‑227 (MQ=255)
                          aaaGACAGCTAGAATAATACCTGTCACCCGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGTCGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAgg                                                                                                                                                                         <  2:204770/251‑1 (MQ=255)
                                  cGAGAATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAACACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCAATg                                                                                                                                                                 >  2:190682/1‑251 (MQ=255)
                                     gaaaaataccTGCCACCCGCGGGGGAAAAGCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTGCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCg                                                                                                                                                              <  2:3472/243‑1 (MQ=255)
                                     gaATAATACCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCTCGCGCAACGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCCCGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCg                                                                                                                                                              >  1:271435/1‑251 (MQ=255)
                                                                     aaaGAAAAACTCAGGGCTCGCGCAACGGCGTTCGTCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCAAATACCAAGCCGTTACGCCACGAAACTGCCCATTGGCATCACGTTCATTATAAAGCTGAGCGAGGGAAAATTCGGCAACCGCGGTAGCCATGCCGGTAAACGCGGAAAccccacaccagaaacc                                                                                                                                 >  2:220542/1‑235 (MQ=255)
                                                                                 aGGGCGCGCGAAACGGCGTGCGTCTGAATTCTGCTGAAAATTAGGCCGTAGGCGATGAGAAAAAAGCCGGCAAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATATCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCAcc                                                                                                                  <  2:93977/251‑1 (MQ=255)
                                                                                         cgaaaCGGCGTTCGCCTGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGACCGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑Taa                                                                                                         <  2:411608/248‑1 (MQ=255)
                                                                                                          gaattCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATTTACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCCGCGCAACGCCg                                                                                        >  1:381316/5‑251 (MQ=255)
                                                                                                           aattCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGGCCGCGCAAcccccg                                                                                       >  1:456193/4‑246 (MQ=255)
                                                                                                                    gAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGTAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACcgg                                                                                                                <  1:235609/218‑1 (MQ=255)
                                                                                                                    gAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAAACCACATCCAGAAGACGGCTCCAGGTCCACcgg                                                                                                                >  2:235609/1‑218 (MQ=255)
                                                                                                                     aaaaTTAGGCCGTAGGCGATGAGAAAAAAGACGCAGAACAGAACCCCCCTCCAGCGCATCCCCCGCCCGCGCCCCATATACCATGCCGGTCCGCCACGAAACGGCCCATTGGCATCACGTTCGTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGcc                                                                             <  2:393739/251‑1 (MQ=255)
                                                                                                                               cgtAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGGAAATTCCCAGGGCACCGCCGGCCAgttaccgtctcccacc                                                                    >  1:244572/3‑236 (MQ=255)
                                                                                                                                    gCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGAGCAACGCCAGCCAGGGTGCCGCAAcccccccggcc                                                              >  1:347794/1‑243 (MQ=255)
                                                                                                                                     cGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑GAATTGCCAGCGCAACCCCGGCCAGGTGGCCGCTccccccccccgcc                                                             >  2:310622/1‑242 (MQ=255)
                                                                                                                                              aaaaGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAACGCCACCCCGGCCAGGCTGCCGCCACCCCcgcggcccccccaaccg                                                    >  2:179556/1‑238 (MQ=255)
                                                                                                                                                aaGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCCGCGCAACGCCGGCCAGGGTGCCGCTACCCCCGCGCGCCCCCAGACTggg                                                  >  1:443645/1‑250 (MQ=255)
                                                                                                                                                    cGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTacac                                              <  1:92594/251‑1 (MQ=255)
                                                                                                                                                         aacaaaacccccccccAGCGCACCCCCAGCCCGCGCCCCATTTCCCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGa                                         <  2:381316/237‑1 (MQ=255)
                                                                                                                                                               aCGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCCGGTTGCCGCTACCCACGCGCGCCGCACGACTTGTCCACAAAGACTgga                                    >  1:149829/1‑248 (MQ=255)
                                                                                                                                                                          gcgaaTCCCCAGCCCGCGCGCCATATCCCATGTCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAACCGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGG‑TAATAGCCAGCGCAACGCCGGCCAGGTTGACGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAAc                        <  1:15719/247‑1 (MQ=255)
                                                                                                                                                                            ccaaccccCGCCCGCGCCCCAAATCCAATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCg                      <  2:119920/247‑1 (MQ=255)
                                                                                                                                                                               tCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCCGGGCAAAACCGGCCAGGGTGCCGCCACCCCCCCCCGCCGCCAGCCTGGTGCACAATGAAAGAAATGAGGGTAAACCCCCt                   >  2:216041/1‑251 (MQ=255)
                                                                                                                                                                                   cagcccccccccaatatACCATGCCGGCCCCCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAACCTGGGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCCCATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCt               <  1:128903/238‑1 (MQ=255)
                                                                                                                                                                                       ccgcccgcCATATCCCAGGCCGGCCCGCCACGAAACGGCCCATTGGCATCACGTTCTTTAAAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTgg           <  1:310622/251‑1 (MQ=255)
                                                                                                                                                                                                atatCCCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCCCGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCt  <  1:245681/251‑1 (MQ=255)
                                                                                                                                                                                                aaatcccaggCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGCCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGG‑TAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACt                                       <  1:325746/205‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                       
AGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAATGCCTGTCACCAGCGGGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGG‑TAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCT  >  NC_000913/7370‑7812

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