Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,487 G→A A80A (GCC→GCT 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,4870GA97.0% 72.6 / ‑6.2 33A80A (GCC→GCTinsH5IS5 family transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/21);  total (12/21)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATG  >  NC_000913/1428319‑1428727
                                                                                                                                                                        |                                                                                                                                                                                                                                                
catcatGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACATAGCATCTGCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGcc                                                                                                                                                                <  1:278528/251‑1 (MQ=11)
   catGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTc                                                                                                                                                             <  2:93953/251‑1 (MQ=11)
                tCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGccc                                                                                                                                                  <  2:224074/249‑1 (MQ=11)
                  ggCCAGCCAGCGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGTTGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCt                                                                                                                                              >  1:64333/1‑251 (MQ=11)
                   gccagccaGCGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCATCATGGTCTCCAGCGGATAGGGCCGTCGTCCATTGCCCGCCtt                                                                                                                                             <  2:61691/251‑1 (MQ=11)
                   gccagccaGCGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCTTCGGCCATTGCCCGCCtt                                                                                                                                             <  1:303693/251‑1 (MQ=11)
                     cagccagcGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCCGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTgg                                                                                                                                           <  1:143365/251‑1 (MQ=11)
                        ccagcGAATGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAAGAGGCGGAAATTCATGATGGTGGTGCGGTCCGGCCAGGCGCTATCCAGGGATAACCGGGCACACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTACATGGCCCCGTCGCTACGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGAGGATCGGGCCGTCGGCCATTGCCCGCCTTGGGa                                                                                                                                         >  2:428932/1‑250 (MQ=11)
                         cagcGATTGATGGTCTTGAACAATTGGCGGGTCAGTTGATGCTGCACCAGCAGGGGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCACTATCCAGGGCTAACCGGGCAAACAGACACATGGAGGCGATTTCGTCCAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGCTAGGGCCGTCGGCCATTGCCCGCCTTGGGAt                                                                                                                                        >  2:435980/1‑250 (MQ=14)
                         cagcGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGAt                                                                                                                                        >  1:367924/1‑250 (MQ=11)
                           gctattcaTGGTCTTGAACAATTGGCGTGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGTGGTCCGGCAAGGCGTTATCCAGAGATAACCGGGCAATCAGACGCATTGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCGCCAGCGGATAGGGCCGTCGACCATTGCCCGCCTTGTGATa                                                                                                                                       <  1:310286/242‑1 (MQ=11)
                            cGATTGATGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATaaaa                                                                                                                                    <  1:430545/251‑1 (MQ=14)
                                 gATGGTCTTGAACAAGTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATCGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCt                                                                                                                               <  1:190373/251‑1 (MQ=14)
                                         tGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGCTGAc                                                                                                                        >  2:14863/1‑250 (MQ=14)
                                                 gggcGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCat                                                                                                               <  2:235322/248‑1 (MQ=14)
                                                         cAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGcc                                                                                                       >  1:392192/1‑251 (MQ=255)
                                                           gTTGATGCTGCTCCAGCAGTTGGCGGAAATTCATGAGGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCAGGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTTGGATAAAACGGCTCGATGACAGCGGACTTATGCTGCCAt                                                                                                     <  2:64333/251‑1 (MQ=14)
                                                                         aGCAGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTc                                                                                       <  2:48015/251‑1 (MQ=255)
                                                                            aGGTGGCGGAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCAt                                                                                    <  2:77656/251‑1 (MQ=255)
                                                                                   ggAAATTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGgag                                                                             <  1:410746/251‑1 (MQ=255)
                                                                                       aTTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCgg                                                                                >  1:92892/1‑244 (MQ=255)
                                                                                       aTTCATGATGGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCgg                                                                                <  2:92892/244‑1 (MQ=255)
                                                                                             gatgGTGGTGCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAAtc                                                                    <  1:339954/250‑1 (MQ=255)
                                                                                                 gtggtgCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGTATGGAGGCGAGTTCGTACAGAGCAGCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGTCCATTGCCCGCCTTGGGAGAAAACTGCTCGATGACAGCGGGCTTATTTTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCttt                                                               <  1:27630/251‑1 (MQ=17)
                                                                                                     tgCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGTACCGTCGCTCAGGTTGTACCAATGCGGCATGCAGTGAATACGCATCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATTGCAGAATCTGCTCCAt                                                                                    <  2:276664/226‑1 (MQ=18)
                                                                                                     tgCGGTCCGGCAAGGCGCTATCCAGGGATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCAt                                                                                    >  1:276664/1‑226 (MQ=255)
                                                                                                                          ccAGGGAGAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATCGCGCCATCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATGCGTAGCATGGTTTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCGTAGTGCt                                       <  2:395816/250‑1 (MQ=14)
                                                                                                                               gATAACCGGGCAAACAGACGCATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTc                                 >  2:190483/1‑251 (MQ=255)
                                                                                                                                                   cATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGTCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCGTAGGGCTGAATTCACTATCGGCGAAGGTGAGtt             <  1:369301/251‑1 (MQ=255)
                                                                                                                                                   cATGGAGGCGATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGtt             >  1:118397/1‑251 (MQ=255)
                                                                                                                                                            gATTTCGTACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCCCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGAGAAGGTGAGTTGATGGCTct    >  2:37437/1‑250 (MQ=255)
                                                                                                                                                             aTTTCGGACAGAGCATCTTCCATGGCACCGTCGCTCAGGTTGTAACAATGCTGCATGCAGTGAATACGAAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGACTCTGATCCATGCGGGAGAGGGAAATCTCTTTTCGGGTCTGATGGCGCTTAGTGATGAATTCACTATCGACGCAGGTGAGTTGATGGCCCat   >  1:195214/1‑251 (MQ=255)
                                                                                                                                                               ttCGTACAGAGTATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGTCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCCTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCatg  <  2:23268/250‑1 (MQ=255)
                                                                                                                                                                        |                                                                                                                                                                                                                                                
CATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATG  >  NC_000913/1428319‑1428727

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