Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,675 C→T Q103Q (CAG→CAA yagA ← integrase core domain‑containing protein YagA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,6750CT95.7% 69.2 / NA 23Q103Q (CAG→CAAyagAintegrase core domain‑containing protein YagA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/17);  minor base G (0/1);  total (5/18)
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

CGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTG  >  NC_000913/281431‑281911
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
cgcgCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTc                                                                                                                                                                                                                                         <  1:27773/250‑1 (MQ=255)
 gcgccggtaatcggaacatggagcCAGGAACAGGAAAAAACGGGAGTGGTCGACCAGCAGGGTGAGCGAATGGCAGAGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCt                                                                                                                                                                                                                                       <  2:144899/229‑1 (MQ=255)
        tCATCGGTACAGGGCGCCAGGCACAGGGAAAAACGGGAGTGGGCGTCCAGCAGGGTGAGCGGATGGCAGCTTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGcc                                                                                                                                                                                                                                 <  1:92685/250‑1 (MQ=255)
         cagcGGTACAGTGCGCAAGGCACAGGGAAAAACGGGAGTGGTAGTCCAGCAGGGTGAGCGGATGGCAGCGACCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAgaca                                                                                                                                                                                                                                <  2:152526/247‑4 (MQ=255)
                 cAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTGAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCAGGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAAt                                                                                                                                                                                                                       <  2:102934/251‑1 (MQ=255)
                                                              gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACgcag                                                                                                                                                                          <  1:45786/251‑1 (MQ=255)
                                                                  gCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAgg                                                                                                                                                                      <  1:117833/251‑1 (MQ=255)
                                                                            gagTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCAGCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATgt                                                                                                                                                            <  2:242904/249‑1 (MQ=255)
                                                                                            ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCATCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGtt                                                                                                                                            <  1:228345/251‑1 (MQ=255)
                                                                                                   ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGaa                                                                                                                                     <  1:82792/251‑1 (MQ=255)
                                                                                                    gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAt                                                                                                                                    <  2:102288/251‑1 (MQ=255)
                                                                                                                               gtgcGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCGTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGaa                                                                                                         <  2:33015/248‑1 (MQ=255)
                                                                                                                                               tCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCtt                                                                                          <  1:247174/250‑1 (MQ=255)
                                                                                                                                                          ggTGGCGGGAATGCCCGGGGAAGGGCCCGGCAGCAGGGCATGGCGGGCCATCAGGTGATGGACGGTGCGGAAGGCGGGCAGGGTGTGCCCGTGGTCCTCGGGCCGGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGCCGTGGGCCAGGCGGAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTGGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAgc                                                                              <  2:113539/250‑1 (MQ=255)
                                                                                                                                                                ggtaaTGACCGGTGAAGCGACAGGCAGCAGGCCATGGCGGGCCAGCGGGTTATGGACGGTGCTTCGGGCGGTCATGGTGTGCCCTTGATACTCGGGCCAGCGCTTAATCTTGCGGGCTCCCCAGGGTTCATGACGGTCATGTGCTATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGTAATGGTGCGGAATGCGCGGTCGGTCCTGAAGTCCGTCGGCACCTTCCTTAGCCCAGCGCTGGa                                                                        <  2:198166/248‑1 (MQ=255)
                                                                                                                                                                   aaTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAt                                                                                                                                    >  1:121835/1‑188 (MQ=255)
                                                                                                                                                                   aaTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAt                                                                                                                                    <  2:121835/188‑1 (MQ=255)
                                                                                                                                                                        ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTg                                                                <  1:53794/251‑1 (MQ=255)
                                                                                                                                                                           gTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAg                                                             <  2:162036/251‑1 (MQ=255)
                                                                                                                                                                             gAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGcc                                                           >  2:97054/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTGGACGGTGCTGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAAc                     >  1:144038/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGAAGGCGGGCATGGTGTGCCCTTGGTCCTCGAGCCCGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCCTGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGCGCCCAGCGCTGGAGCCCCTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCCGAGGGCACGGATGGTCGcc          >  2:232419/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggCATGGTGTGCCCTTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTg  >  1:99406/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
CGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTG  >  NC_000913/281431‑281911

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