Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,944,006 A→G 10.8% F138F (TTT→TTC)  splA ← S1B family serine protease SplA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,944,0060AG10.8% 146.0 / 6.6 54F138F (TTT→TTC) splAS1B family serine protease SplA
Reads supporting (aligned to +/- strand):  ref base A (16/31);  new base G (3/3);  total (20/34)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCATCAAATTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGTTTCCTCCGCCTTTACCTTTACT  >  CP000730/1943872‑1944146
                                                                                                                                      |                                                                                                                                            
tACGCATCAAATTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCCAAtttt                                                                                                                                        >  2:490699/1‑141 (MQ=255)
      tCAAATTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTAtaa                                                                                                                                  <  2:380664/141‑1 (MQ=255)
      tCAAATTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTAtaa                                                                                                                                  <  1:152560/141‑1 (MQ=255)
       cAAATTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATACCAGCAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTAtaac                                                                                                                                 <  2:24615/141‑1 (MQ=255)
       cAAATTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTAtaac                                                                                                                                 >  2:10078/1‑141 (MQ=255)
        aaaTTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTAtaact                                                                                                                                <  2:408549/141‑1 (MQ=255)
        aaaTTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTAtaact                                                                                                                                >  2:683687/1‑141 (MQ=255)
          aTTCCATAAACGTTCCACTGAGATGGTTAGTCGTGCCTGTCGATTCAAACATTTTATATTTTGTGTGTGCACCCTTTGGATAACCAAGAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTAtaactaa                                                                                                                              >  2:175732/1‑141 (MQ=255)
              cATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGtt                                                                                                                          >  2:305842/1‑141 (MQ=255)
                tAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCGAATTTTGTATAACTAACGTTCt                                                                                                                        >  1:329085/1‑141 (MQ=255)
                tAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTAAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCt                                                                                                                        <  1:530701/141‑1 (MQ=255)
                 aaaCGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCtt                                                                                                                       >  2:691020/1‑141 (MQ=255)
                 aaaCGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCtt                                                                                                                       <  1:61240/141‑1 (MQ=255)
                  aaCGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCtta                                                                                                                      <  2:421198/141‑1 (MQ=255)
                         ctttcttatGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAAt                                                                                                               <  2:562640/135‑1 (MQ=255)
                           cTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTc                                                                                                             >  1:106985/1‑141 (MQ=255)
                           cTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTc                                                                                                             <  2:302819/141‑1 (MQ=255)
                            tGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCGAATTTTGTATAACTAACGTTCTTATTAAAATTCa                                                                                                            >  1:381507/1‑141 (MQ=255)
                            tGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATGTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCGAATTTTGTATAACTAACGTTCTTATTAAAATTCa                                                                                                            >  1:381438/1‑141 (MQ=255)
                                     tAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTg                                                                                                   <  1:189901/141‑1 (MQ=255)
                                     tAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTg                                                                                                   >  1:140792/1‑141 (MQ=255)
                                     tAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTg                                                                                                   <  2:516804/141‑1 (MQ=255)
                                       aTCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTa                                                                                                 >  2:395024/1‑141 (MQ=255)
                                              cTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGttt                                                                                          >  2:339184/1‑141 (MQ=255)
                                                 tCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTcat                                                                                       >  1:464326/1‑141 (MQ=255)
                                                      tCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTcatgaaca                                                                                  <  2:506901/141‑1 (MQ=255)
                                                       cAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTcatgaacat                                                                                 <  1:115283/141‑1 (MQ=255)
                                                           cATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCatgaacatgaac                                                                             <  1:322347/141‑1 (MQ=255)
                                                              tttATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTAt                                                                          <  2:292422/141‑1 (MQ=255)
                                                                 atatTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGc                                                                       <  1:550358/141‑1 (MQ=255)
                                                                 atatTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGc                                                                       <  2:229826/141‑1 (MQ=255)
                                                                 atatTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGc                                                                       >  1:631076/1‑141 (MQ=255)
                                                                   atTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCaa                                                                     >  2:78124/1‑141 (MQ=255)
                                                                                   cTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCggg                                                     <  1:290920/141‑1 (MQ=255)
                                                                                   cTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCggg                                                     <  2:587340/141‑1 (MQ=255)
                                                                                     ttGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGat                                                   >  1:233030/1‑141 (MQ=255)
                                                                                            aCCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTAc                                            <  1:691020/141‑1 (MQ=255)
                                                                                             ccAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACa                                           >  2:655992/1‑141 (MQ=255)
                                                                                                aTAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATg                                        <  1:347243/141‑1 (MQ=255)
                                                                                                 tAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGt                                       >  1:396725/1‑141 (MQ=255)
                                                                                                  aaCAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTc                                      <  1:395024/141‑1 (MQ=255)
                                                                                                          ttCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAAcgt                              <  1:300276/141‑1 (MQ=255)
                                                                                                          ttCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAAcgt                              <  1:273792/141‑1 (MQ=255)
                                                                                                             tATCTTTCACTTTTGCTCCGTCTGCGAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAAcgtcgt                           <  2:329085/141‑1 (MQ=255)
                                                                                                              aTCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTa                          <  2:58740/141‑1 (MQ=255)
                                                                                                              aTCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTa                          <  1:175732/141‑1 (MQ=255)
                                                                                                                cTTTCACTTTTGCTCCGTCTGCGAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGt                        <  2:381438/141‑1 (MQ=255)
                                                                                                                cTTTCACTTTTGCTCCGTCTGCGAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGt                        <  2:381507/141‑1 (MQ=255)
                                                                                                                 tttCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGtt                       <  1:412365/141‑1 (MQ=255)
                                                                                                                   tCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAGCGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCTTAGTTtc                     >  1:344474/1‑141 (MQ=255)
                                                                                                                     aCTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGTTtcct                   <  1:30495/141‑1 (MQ=255)
                                                                                                                      cTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGTTtcctc                  <  1:490699/141‑1 (MQ=255)
                                                                                                                              ccGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGTTTCCTCCGccttta          <  2:631076/141‑1 (MQ=255)
                                                                                                                                  cTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGTTTCCTCCGcctttacctt      <  2:326213/141‑1 (MQ=255)
                                                                                                                                      aaaTTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGTTTCCTCCGCCTTTACCTTTACt  >  1:576758/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                            
TACGCATCAAATTCCATAAACGTTCCACTGATATGGTTAATCGTTCCTGTCGATTCAAACATTTTATATTTTGTTTGTGCACCCTTTGGATAACCAATAACAGAAATTCTATCTTTCACTTTTGCTCCGTCTGCAAATTTTGTATAACTAACGTTCTTATTAAAATTCAAACCTTCTGTACTTGTTTCATGAACATGAACTATCGCAAGGTCTTCTTTTCCGGGATATTCTACAATGTCTTTAACGTCGTAGTTTCCTCCGCCTTTACCTTTACT  >  CP000730/1943872‑1944146

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