Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,669 G→T R20R (CGA→AGA)  insH‑5 ← Rac prophage; IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,6690GT98.2% 170.4 / NA 56R20R (CGA→AGA) insH‑5Rac prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (28/27);  total (28/28)
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

AGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGG  >  NC_000913/1428534‑1428800
                                                                                                                                       |                                                                                                                                   
aGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTATGGTc                                                                                                                                 >  1:526470/1‑140 (MQ=255)
            gcaTGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTg                                                                                                                               >  1:656455/1‑130 (MQ=255)
            gcaTGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTg                                                                                                                               <  2:656455/130‑1 (MQ=255)
             caTGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCt                                                                                                                       >  1:555602/1‑137 (MQ=255)
             caTGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCt                                                                                                                       <  2:555602/137‑1 (MQ=255)
              aTGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTact                                                                                                                   >  2:234568/1‑138 (MQ=255)
                  tctcCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTActgctg                                                                                                               >  1:16390/1‑140 (MQ=255)
                   ctcCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCt                                                                                                                <  2:478523/138‑1 (MQ=255)
                   ctcCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCt                                                                                                                >  1:478523/1‑138 (MQ=255)
                     ccAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTg                                                                                                                  <  1:905461/134‑1 (MQ=255)
                     ccAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTg                                                                                                                  >  2:905461/1‑134 (MQ=255)
                            aTAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTAt                                                                                                     >  2:648304/1‑140 (MQ=255)
                            aTAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTAt                                                                                                     >  2:780414/1‑140 (MQ=255)
                            aTAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTAt                                                                                                     >  1:766804/1‑140 (MQ=255)
                                ggCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGc                                                                                                 <  2:232423/140‑1 (MQ=17)
                                     tCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAgg                                                                                            >  1:802367/1‑140 (MQ=255)
                                        gCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGtaa                                                                                         >  2:785676/1‑138 (MQ=14)
                                                 cGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATCCACTATCGGCGAAGGTGAGTTGATGGc                                                                                <  1:125799/140‑1 (MQ=14)
                                                  gCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACt                                                                               >  1:116760/1‑140 (MQ=11)
                                                  gCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCt                                                                               >  1:11556/1‑140 (MQ=255)
                                                   ccTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCt                                                                               >  1:762423/1‑139 (MQ=255)
                                                   ccTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCGGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCt                                                                               <  2:762423/139‑1 (MQ=255)
                                                         gATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                        <  1:234568/140‑1 (MQ=255)
                                                         gATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                        <  2:526470/140‑1 (MQ=255)
                                                         gATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCatgatg                                                                        <  2:324432/140‑1 (MQ=255)
                                                           tAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTc                                                                      <  2:841419/140‑1 (MQ=255)
                                                              aaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGAt                                                                                   >  1:259499/1‑124 (MQ=255)
                                                              aaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGAt                                                                                   <  2:259499/124‑1 (MQ=14)
                                                              aaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCt                                                                               >  1:227384/1‑128 (MQ=255)
                                                              aaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCt                                                                               <  2:227384/128‑1 (MQ=255)
                                                                cGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCctct                                                                 <  2:767669/140‑1 (MQ=255)
                                                                        tGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATgcg                                                         >  2:395274/1‑140 (MQ=255)
                                                                           cAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGAGATGCGCTc                                                      <  1:656011/140‑1 (MQ=17)
                                                                              cGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAGTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCgg                                                   <  1:212287/140‑1 (MQ=17)
                                                                              cGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCgg                                                   <  2:16390/140‑1 (MQ=255)
                                                                              cGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCgg                                                   <  1:780414/140‑1 (MQ=255)
                                                                              cGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCgg                                                   >  1:812385/1‑140 (MQ=255)
                                                                                gTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGAt                                                 <  2:16182/140‑1 (MQ=255)
                                                                                gTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGAt                                                 <  2:836076/140‑1 (MQ=255)
                                                                                gTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGAt                                                 <  2:802367/140‑1 (MQ=255)
                                                                                agcagatTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGAt                                                 >  2:415505/6‑140 (MQ=17)
                                                                                    gatTCTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAat                                             <  1:415505/139‑1 (MQ=255)
                                                                                          gCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATatgatg                                       <  1:88333/140‑1 (MQ=255)
                                                                                           ccATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATAtgatga                                      >  1:845313/1‑140 (MQ=17)
                                                                                           ccATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATAtgatga                                      >  1:231750/1‑140 (MQ=17)
                                                                                              tGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCgg                                                   <  2:87334/124‑1 (MQ=255)
                                                                                              tGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCgg                                                   >  1:87334/1‑124 (MQ=255)
                                                                                                          cTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCa                                 <  1:77593/130‑1 (MQ=255)
                                                                                                          cTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCa                                 >  2:77593/1‑130 (MQ=18)
                                                                                                          cTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGaa                        >  2:493355/1‑139 (MQ=18)
                                                                                                          cTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGaa                        <  1:493355/139‑1 (MQ=255)
                                                                                                           tCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                       <  2:440432/139‑1 (MQ=255)
                                                                                                           tCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAAc                       >  1:440432/1‑139 (MQ=18)
                                                                                                              aTGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGt                   >  1:809549/1‑140 (MQ=18)
                                                                                                               tGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGtt                  >  1:539100/1‑140 (MQ=255)
                                                                                                                gCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTc                 >  2:80520/1‑140 (MQ=18)
                                                                                                                  gggACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGc               <  2:231750/140‑1 (MQ=255)
                                                                                                                     acaagAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCAcc            <  1:395274/136‑1 (MQ=255)
                                                                                                                               ctctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAggg  >  1:860135/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGG  >  NC_000913/1428534‑1428800

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