Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,396,138 C→T D21D (GAC→GAT insH‑4 → IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,396,1380CT92.3% 68.6 / ‑1.9 26D21D (GAC→GATinsH‑4IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (18/6);  total (20/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAG  >  NC_000913/1396002‑1396177
                                                                                                                                        |                                       
aGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGt                                      <  2:585717/140‑1 (MQ=255)
  atcatcAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGa                                    >  2:492162/1‑140 (MQ=255)
           aGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTcagcag                           <  1:646917/140‑1 (MQ=255)
            gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAAt                                  <  2:124285/132‑1 (MQ=255)
            gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAAt                                  >  1:124285/1‑132 (MQ=255)
            gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGt                          >  2:527197/1‑140 (MQ=255)
             tCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTa                         >  2:151933/1‑140 (MQ=255)
               aGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAg                       >  1:386723/1‑140 (MQ=255)
               aGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAg                       >  1:1216471/1‑140 (MQ=255)
                 ggaaggaaGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAgcg                     >  1:1437660/1‑140 (MQ=255)
                 ggaaggaaGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAgcg                     >  2:1006944/1‑140 (MQ=255)
                 ggaaggaaGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAgcg                     >  2:1303761/1‑140 (MQ=255)
                  gaaggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAgcgc                    >  1:1230332/1‑140 (MQ=255)
                     ggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGt                 <  1:1303761/140‑1 (MQ=255)
                     ggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGt                 <  1:1006944/140‑1 (MQ=255)
                      gaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTc                >  1:1176477/1‑140 (MQ=255)
                       aagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCCGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCa               >  2:1324756/1‑140 (MQ=255)
                          tAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCAGAc            <  1:747165/140‑1 (MQ=255)
                           aGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCAGAcc           >  1:1150559/1‑140 (MQ=38)
                            gCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCAGACCa          >  1:465822/1‑140 (MQ=37)
                            gCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCa          >  1:1048796/1‑140 (MQ=37)
                             cAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCAGACCAg         >  2:659945/1‑140 (MQ=37)
                                cATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaa      >  1:1544721/1‑140 (MQ=18)
                                cATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaa      >  2:802507/1‑140 (MQ=18)
                                cATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaa      >  2:765504/1‑140 (MQ=33)
                                    aGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGATAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAagag  >  1:818726/1‑140 (MQ=11)
                                                                                                                                        |                                       
AGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAG  >  NC_000913/1396002‑1396177

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