Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,413,092 C→T A48A (GCG→GCA intR ← putative integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,0920CT96.5% 102.7 / ‑6.0 29A48A (GCG→GCAintRputative integrase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (12/16);  total (12/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

TATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAGT  >  NC_000913/1412845‑1413334
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
tatCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGct                                                                                                                                                                                                                                                 >  2:331473/1‑251 (MQ=255)
tatCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGct                                                                                                                                                                                                                                                 >  2:382505/1‑251 (MQ=255)
             aaTATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGCGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCATCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTCTTTCATGAACTATCAGTGCTCTTAAATTTCCgg                                                                                                                                                                                                                                    >  2:229837/1‑251 (MQ=255)
               tatGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGcc                                                                                                                                                                                                                                  >  2:71280/1‑251 (MQ=255)
                     atatCGGGGTACTTTCACAAACTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTATCTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCGCTTAAATTTCCGGCCTTTTTa                                                                                                                                                                                                                            <  1:8292/251‑1 (MQ=255)
                       atCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAat                                                                                                                                                                                                                          <  2:236967/251‑1 (MQ=255)
                       atCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAat                                                                                                                                                                                                                          <  1:445779/251‑1 (MQ=255)
                                  ttCACAAATTATTAATCTGAGTGTTTTTAATTGTTATTTCGTTTTCTTCCTTGTGTTTGTTGTTTACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGttt                                                                                                                                                                                                                <  1:26620/250‑1 (MQ=255)
                                    cacaAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGt                                                                                                                                                                                                              <  2:285927/250‑1 (MQ=255)
                                                        gtTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTACTTTTTGTCCAGATATCACAAAGCTCTCCGAAAGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGt                                                                                                                                                                                         >  2:356990/1‑251 (MQ=255)
                                                                  gtgATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCTCAGCCTTTTAATGTTTCCCGACAt                                                                                                                                                                               >  1:9963/1‑251 (MQ=255)
                                                                                  cATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAAc                                                                                                                                                                <  1:277074/250‑1 (MQ=255)
                                                                                              tgttAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGATTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATg                                                                                                                                                    <  2:253491/250‑1 (MQ=255)
                                                                                               gttAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCg                                                                                                                                                  <  1:131158/251‑1 (MQ=255)
                                                                                                                         ccAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGacac                                                                                                                        >  1:294409/1‑251 (MQ=255)
                                                                                                                                      aGCTCTACGAACGTTTTTATGACTCTCGTTGTCACCCTGTTTGCCCCAGTGCTGGACTTGGGAAAACGCCTTAAATACTCACATTCACCGGCGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCATCGAAACATGAACCAGAGGCGAATGCATCTACCTCTAAACTCGACACCTGTTGGTAAttt                                                                                                           >  2:156982/1‑251 (MQ=255)
                                                                                                                                                             tctcGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGCCATATCATGAGTCTTTGAt                                                                                     >  2:29734/1‑250 (MQ=255)
                                                                                                                                                              ctcGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGAt                                                                                     <  1:185061/249‑1 (MQ=255)
                                                                                                                                                                         cATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATc                                                                        <  1:331473/251‑1 (MQ=255)
                                                                                                                                                                                           ggACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCa                                                       <  1:316510/250‑1 (MQ=255)
                                                                                                                                                                                           ggACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCa                                                       <  2:299881/250‑1 (MQ=255)
                                                                                                                                                                                                      aaCGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTc                                                                                                                                                                        <  2:75404/126‑1 (MQ=255)
                                                                                                                                                                                                      aaCGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTc                                                                                                                                                                        >  1:75404/1‑126 (MQ=255)
                                                                                                                                                                                                             tAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTg                                       <  2:273577/248‑1 (MQ=255)
                                                                                                                                                                                                                  aCTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGAAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTg                                <  1:382505/250‑1 (MQ=255)
                                                                                                                                                                                                                           cACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTaa                      <  2:40323/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTa                       >  1:145910/1‑247 (MQ=255)
                                                                                                                                                                                                                              cGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGGTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTa                       <  2:145910/247‑1 (MQ=255)
                                                                                                                                                                                                                                               aCTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAgg  >  2:254696/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
TATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAGT  >  NC_000913/1412845‑1413334

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