Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,163,633 C→T 8.9% R569H (CGT→CAT)  PP_5405 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,163,6330CT 8.9% 121.4 / 3.2 45R569H (CGT→CAT) PP_5405transposase
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (2/2);  total (24/21)
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.33e-01

CGTCTCTGAGATCAGCAGCAGCTCCTCGACGTACCTGATCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCTTTCTTATCGCAGCTGGTCG  >  NC_002947/6163484‑6163778
                                                                                                                                                     |                                                                                                                                                 
cGTCTCTGAGATCAGCAGCAGCTCCTCGACGTACCTGATCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGAc                                                                                                                                                   <  2:603843/150‑1 (MQ=255)
cGTCTCTGAGATCAGCAGCAGCTCCTCGACGTACCTGATCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGAc                                                                                                                                                   <  2:123112/150‑1 (MQ=255)
                              gTACCTGATCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAgg                                                                                                                     >  2:1108233/1‑150 (MQ=255)
                                 ccTGATCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTAt                                                                                                                  >  2:1647266/1‑150 (MQ=255)
                                    gATCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATtga                                                                                                               <  2:542784/150‑1 (MQ=255)
                                     aTCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCATTTAACTTCTATCCATATATGACGAACACACATGACCAAATTATATACGGATGTATtgat                                                                                                              <  2:746659/150‑1 (MQ=255)
                                         aTAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGaca                                                                                                          >  1:313897/1‑150 (MQ=255)
                                          tAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGacac                                                                                                         >  1:1123598/1‑150 (MQ=255)
                                              ttGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTg                                                                                                     <  1:969311/150‑1 (MQ=255)
                                                gAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGcc                                                                                                   >  1:591662/1‑150 (MQ=255)
                                                    aTCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCggac                                                                                               >  1:216889/1‑148 (MQ=255)
                                                         gCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGCAGTCTAACTGGTATAAAACCCCTTTAACTTCCATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGccacca                                                                                          >  2:801638/1‑150 (MQ=255)
                                                             ggAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAg                                                                                                                      >  1:423148/1‑118 (MQ=255)
                                                             ggAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAg                                                                                                                      <  2:423148/118‑1 (MQ=255)
                                                                   atCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAatat                                                                                <  2:313897/150‑1 (MQ=255)
                                                                       aCGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCt                                                                            <  1:1108233/150‑1 (MQ=255)
                                                                           aGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGATGAACACACATGACCAAAATATATACGGAGGTATtgatga                                                                                                            >  2:1230739/1‑114 (MQ=255)
                                                                           aGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGATGAACACACATGACCAAAATATATACGGAGGTATtgatga                                                                                                            <  1:1230739/114‑1 (MQ=255)
                                                                              tttGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTttctt                                                                     <  1:147550/150‑1 (MQ=255)
                                                                                       tACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTc                                                                              <  1:50992/132‑1 (MQ=255)
                                                                                       tACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTc                                                                              >  2:50992/1‑132 (MQ=255)
                                                                                       tACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCAc                                                            >  2:697833/1‑150 (MQ=255)
                                                                                                     aTGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACcaaatc                                                                                      >  1:1523036/1‑110 (MQ=255)
                                                                                                     aTGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACcaaatc                                                                                      <  2:1523036/110‑1 (MQ=255)
                                                                                                         aGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGcc                                          <  1:697833/150‑1 (MQ=255)
                                                                                                             tAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGATGAACACACATGACCAAAAtatata                                                                                                                           <  2:1613828/65‑1 (MQ=255)
                                                                                                             tAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGATGAACACACATGACCAAAAtatata                                                                                                                           >  1:1613828/1‑65 (MQ=255)
                                                                                                                 tGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTaaaaa                                  >  1:1234581/1‑150 (MQ=255)
                                                                                                                   gTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATtgatga                                                                                                            >  2:795018/1‑74 (MQ=255)
                                                                                                                   gTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATtgatga                                                                                                            <  1:795018/74‑1 (MQ=255)
                                                                                                                   gTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGac                                                                                                           >  1:465527/1‑75 (MQ=255)
                                                                                                                   gTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGac                                                                                                           <  2:465527/75‑1 (MQ=255)
                                                                                                                    tgtttaaCCCCTTTAACGGCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAg                               >  2:1297419/6‑150 (MQ=255)
                                                                                                                           ccccTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATtgatga                                                                                                            >  2:116989/1‑66 (MQ=255)
                                                                                                                           ccccTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATtgatga                                                                                                            <  1:116989/66‑1 (MQ=255)
                                                                                                                           ccccTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTACTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGt                        >  2:1120523/1‑150 (MQ=255)
                                                                                                                              cTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCt                     >  2:1465228/1‑150 (MQ=255)
                                                                                                                              cTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCt                     <  1:1647266/150‑1 (MQ=255)
                                                                                                                              cTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCt                     <  2:591662/150‑1 (MQ=255)
                                                                                                                                 tAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACAGCATCGGTATTTTAGCCCGTAAAAATAGACACCGTCCTTTc                  <  2:216889/150‑1 (MQ=255)
                                                                                                                                   aCTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCTTTCtt                >  2:1446664/1‑150 (MQ=255)
                                                                                                                                     ttCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATtgat                                                                                                              <  1:815134/54‑1 (MQ=255)
                                                                                                                                     ttCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATtgat                                                                                                              >  2:815134/1‑54 (MQ=255)
                                                                                                                                     ttCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCTTTCTTAt              <  1:1446664/150‑1 (MQ=255)
                                                                                                                                             atatatGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCTTTCTTATCGCAGCTg      >  1:880269/1‑150 (MQ=255)
                                                                                                                                                tatGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCTTTCTTATCGCAGCTGGTc   <  1:1297419/150‑1 (MQ=255)
                                                                                                                                                 ttGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCTTTCTTATCGCAGCTGGTCg  >  1:1196459/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGTCTCTGAGATCAGCAGCAGCTCCTCGACGTACCTGATCAATAGCTTGAAGATCACGCAAGGAAATATCGACGGAGCTTTGCGGACTACGCGATGGCCGCATGAAGTCTAACTGGTATAAAACCCCTTTAACTTCTATCCATATATGACGAACACACATGACCAAAATATATACGGAGGTATTGATGACACCTTGCCGGGCCACCAAATCAAATATTCCTGTTTCTTCCAATTCACGGGATCGGTATTTTCGCCCGTAAAAATAGACACCGTCCTTTCTTATCGCAGCTGGTCG  >  NC_002947/6163484‑6163778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑