Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,109 C→T L206L (CTG→CTA insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,1090CT100.0% 80.4 / NA 31L206L (CTG→CTAinsH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/14);  total (17/14)

TCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATG  >  NC_000913/1427973‑1428234
                                                                                                                                        |                                                                                                                             
tCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAg                                                                                                                             <  2:1011587/139‑1 (MQ=255)
 cAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGc                                                                                                                            >  2:606244/1‑139 (MQ=255)
  aaGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCt                                                                                                                           >  2:109956/1‑139 (MQ=255)
     gTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGAt                                                                                                                        <  1:109956/139‑1 (MQ=255)
                           cGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCcgc                                                                                                  >  1:444997/1‑139 (MQ=21)
                                  cAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGgtggtg                                                                                           <  2:858434/139‑1 (MQ=21)
                                        gTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAgg                                                                                     >  2:855554/1‑139 (MQ=25)
                                            acaTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCtgt                                                                                 <  1:522910/139‑1 (MQ=21)
                                                   ccTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAg                                                                          >  1:1223181/1‑139 (MQ=37)
                                                        gCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCAc                                                                     >  1:996225/1‑139 (MQ=21)
                                                                 tcGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCAt                                                            <  1:860992/139‑1 (MQ=25)
                                                                        gtgGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACAcc                                                     >  2:425861/1‑139 (MQ=25)
                                                                              cTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAAtgtg                                               >  2:389407/1‑139 (MQ=21)
                                                                                   tGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCt                                          >  1:222026/1‑139 (MQ=17)
                                                                                           ggCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGc                                    >  1:127615/1‑137 (MQ=17)
                                                                                           ggCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGc                                    <  2:127615/137‑1 (MQ=17)
                                                                                           ggCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCa                                  >  1:1186150/1‑139 (MQ=17)
                                                                                           ggCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCa                                  >  1:1304234/1‑139 (MQ=18)
                                                                                             cATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCaaa                                >  2:1220344/1‑139 (MQ=17)
                                                                                                cGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTg                             >  1:169389/1‑139 (MQ=17)
                                                                                                   cTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCa                          <  1:855554/139‑1 (MQ=18)
                                                                                                            aTTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGcc                 <  2:444997/139‑1 (MQ=17)
                                                                                                            aTTGCTCCTCTCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGcc                 >  2:796491/1‑139 (MQ=17)
                                                                                                                   ctctCCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTg          <  1:645468/139‑1 (MQ=17)
                                                                                                                     ctcCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCAt                                                            <  2:581869/87‑1 (MQ=25)
                                                                                                                     ctcCATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCAt                                                            >  1:581869/1‑87 (MQ=25)
                                                                                                                       ccATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCcg                                                                                                   <  1:759144/46‑1 (MQ=255)
                                                                                                                       ccATGAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCcg                                                                                                   >  2:759144/1‑46 (MQ=255)
                                                                                                                           gAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATg  <  2:1223181/139‑1 (MQ=14)
                                                                                                                           gAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATg  <  2:169389/139‑1 (MQ=14)
                                                                                                                           gAAGCAGATTACCTAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATg  <  1:389407/139‑1 (MQ=14)
                                                                                                                                        |                                                                                                                             
TCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATG  >  NC_000913/1427973‑1428234

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