Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,107 C→T G207D (GGT→GAT)  insH5 ← IS5 family transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,1070CT96.4% 69.2 / ‑3.9 28G207D (GGT→GAT) insH5IS5 family transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/13);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

CACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCA  >  NC_000913/1427860‑1428345
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
cACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCa                                                                                                                                                                                                                                             <  2:62103/251‑1 (MQ=255)
              gCCGCTTGATGATGTGACACGGGTGGTCCGTCCTGGCACGGATGCTGGCTTTCCTGTATTCCATGATGTTGGCAGTTTTGTTATTGCGCGGATTCTGCTTCATGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCCTCGGCTGCGACAAATTGCTCCTCTCCATGACGCAGATTCTCCAGCTGATTGAGGTcc                                                                                                                                                                                                                               >  2:142026/1‑250 (MQ=14)
                cGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATg                                                                                                                                                                                                                             >  2:95145/1‑251 (MQ=255)
                        gatgCGAAAAGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGCTTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAACTTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTgg                                                                                                                                                                                                                     >  2:42321/1‑251 (MQ=255)
                                                     ggATGCTGGCTTTCATGTATTCGATGTTGATGGCGGTTTTGTTCTTGCGCGGATGCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCAGCGGCTGAGACAAATTGCTCCTCTCCATGAATCAGATTATCCAGCTGATTGAGGTCATGCTTGTTGGCGGCGTTGCGGACCAGGCTGTTGGTCAgg                                                                                                                                                                                        <  2:186583/251‑1 (MQ=14)
                                                           tGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCActc                                                                                                                                                                                  >  2:250353/1‑251 (MQ=255)
                                                            ggCTTTCATGTATTCGATGTTGATTGCCGTTTTGTTCTTGCTCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCCCCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCATGCCActct                                                                                                                                                                                 <  1:50797/251‑1 (MQ=255)
                                                              cTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTg                                                                                                                                                                               >  1:202921/1‑251 (MQ=255)
                                                                        ttCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGCCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCATGCCACTCTTGGCATCGacac                                                                                                                                                                     >  1:283525/1‑251 (MQ=255)
                                                                                 gATGGCCGTTTTGTTCTTTCGCGGATTCGGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCGCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCGCCATGAAGCAGATTATCCAGTTGATTGAGGTCATGCTCGTGGGCCGGGGTGTTGACCATGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTgg                                                                                                                                                             <  1:185471/250‑1 (MQ=14)
                                                                                                     gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCAc                                                                                                                                        <  1:241797/251‑1 (MQ=255)
                                                                                                              gCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCt                                                                                                                               <  1:95145/251‑1 (MQ=38)
                                                                                                                       tttttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTc                                                                                                                      >  2:20881/1‑251 (MQ=25)
                                                                                                                         tttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggtgg                                                                                                                                                                                                            >  1:90445/1‑163 (MQ=255)
                                                                                                                         tttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCGCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggtgg                                                                                                                                                                                                            <  2:90445/163‑1 (MQ=255)
                                                                                                                            cccTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGAt                                                                                                                  <  2:259510/249‑1 (MQ=25)
                                                                                                                              cTTGCCGGAACGCTCGGCGATCAGCCAGTCTACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGCCAAATTGCCCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCCCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCCCTCTTGGCATCGACCCCAATGTGGGCCTTCATGCCAAAGGGCCACTGATTGCTTTTCTTGGTCTGATGc                                                                                                                >  2:17404/1‑250 (MQ=11)
                                                                                                                                          cTCGGCGCTCAGCCAGTCCACATCCCCCTCGGCCAGCTCATCGCGCTGTGGCGCTCCTTGGTAGCCGGCACCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGAGGTGACCAGGCTGTGGGTCAGGCCACTCTGGGCATCGACACCAATGTGGGCCTTCATGCCCAAGAGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATcgc                                                                                                   >  2:193012/1‑251 (MQ=11)
                                                                                                                                                     gCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCGCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATTCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGGCTGTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGctct                                                                                        <  1:42321/251‑1 (MQ=255)
                                                                                                                                                        aGTCCACATCCACCTCGGCCAGCTCCTCGAGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTg                                                                                     >  2:274802/1‑251 (MQ=18)
                                                                                                                                                            cacaTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCt                                                                                 <  1:20881/251‑1 (MQ=18)
                                                                                                                                                                          ccAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGt                                                                   >  1:243405/1‑251 (MQ=21)
                                                                                                                                                                                                         gCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGc                                    >  1:246569/1‑251 (MQ=37)
                                                                                                                                                                                                            ggCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCtt                                 >  2:63240/1‑251 (MQ=37)
                                                                                                                                                                                                              cATCGGCTGAGACAAATTGCGCCTCTCCATGAAGCAGAGTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTgg                               <  2:1293/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGt                             <  1:63240/250‑1 (MQ=38)
                                                                                                                                                                                                                                   cctcTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTc           <  1:224102/250‑1 (MQ=39)
                                                                                                                                                                                                                                            gAAGCAGATTATACAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGGATCGACACCAATGTGGGCCGTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGgccagcca  <  2:283525/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
CACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCA  >  NC_000913/1427860‑1428345

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