Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,211 T→A 8.3% T53T (ACT→ACA ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2110TA 8.3% 1950.9 / 121.1 599T53T (ACT→ACAssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (258/291);  new base A (35/15);  total (293/306)
Fisher's exact test for biased strand distribution p-value = 1.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 631 total aligned reads displayed.

CAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGT  >  CP009273/4264145‑4264277
                                                                  |                                                                  
cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAATGAACAGACt                                                                    <  1:22298790/67‑1 (MQ=255)
 aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGTCCGGCGAGATGAAAGAACAGACTGa                                                                  <  1:26592216/68‑1 (MQ=255)
 aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                                  <  1:10247858/68‑1 (MQ=255)
 aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                                  <  1:35554196/68‑1 (MQ=255)
 aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                                  <  1:14852642/68‑1 (MQ=255)
 aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                                  <  1:16243843/68‑1 (MQ=255)
  aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                                 >  1:7564668/1‑68 (MQ=255)
  aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGAATGaa                                                                 >  1:17823267/1‑68 (MQ=255)
    aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                               >  1:29517767/1‑68 (MQ=255)
    aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                               >  1:25019740/1‑68 (MQ=255)
    aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                               >  1:32566536/1‑68 (MQ=255)
    aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                               >  1:34840057/1‑68 (MQ=255)
     ttACGCTGGCTACTTCCGAATACTGGCGTGATAAAGCGACCGCCGAGATGAAAGCACAGACAGAATg                                                               >  1:20733941/1‑67 (MQ=255)
      tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                             <  1:19426916/68‑1 (MQ=255)
      tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                             <  1:33551752/68‑1 (MQ=255)
      tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                             <  1:39549710/68‑1 (MQ=255)
      tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                             <  1:72072/68‑1 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:31186031/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:15811337/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:24434683/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:29321980/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:7977062/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:12984497/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:2092179/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:21201134/1‑68 (MQ=255)
         gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:8980834/1‑68 (MQ=255)
          cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:36832943/1‑67 (MQ=255)
          cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                          >  1:4107690/1‑67 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:3022899/68‑1 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:23852205/68‑1 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:12108898/68‑1 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:22465396/68‑1 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:30649473/68‑1 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:11001690/68‑1 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:39111412/68‑1 (MQ=255)
           tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                        <  1:27865474/68‑1 (MQ=255)
              cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                     >  1:35707368/1‑68 (MQ=255)
              cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                     >  1:15595623/1‑68 (MQ=255)
              cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                                      <  1:31677443/67‑1 (MQ=255)
                aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                   <  1:17026821/68‑1 (MQ=255)
                aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                   <  1:9377131/68‑1 (MQ=255)
                aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                   <  1:18640835/68‑1 (MQ=255)
                   tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                   >  1:13755039/1‑65 (MQ=255)
                   tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                   >  1:16109189/1‑65 (MQ=255)
                     cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                              <  1:21721735/68‑1 (MQ=255)
                     cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                              <  1:33845979/68‑1 (MQ=255)
                     cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                              <  1:25845838/68‑1 (MQ=255)
                     cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                              <  1:35188106/68‑1 (MQ=255)
                         tCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTGCg                                           >  1:9651076/1‑67 (MQ=255)
                         agcTGGCGTGTCAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCgg                                          <  1:14519398/66‑1 (MQ=255)
                          gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                         >  1:10078067/2‑68 (MQ=255)
                          gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                         >  1:11869239/2‑68 (MQ=255)
                            tGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCa                                        <  1:38418732/67‑1 (MQ=255)
                             ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaa                                       >  1:6267616/1‑67 (MQ=255)
                               cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                       >  1:4820425/1‑65 (MQ=255)
                               cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAAc                                     <  1:13358568/67‑1 (MQ=255)
                                gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                   <  1:3789926/68‑1 (MQ=255)
                                  gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                                  <  1:24739385/67‑1 (MQ=255)
                                    tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                <  1:37306320/67‑1 (MQ=255)
                                       acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                            <  1:26167898/65‑1 (MQ=255)
                                       acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                            <  1:5503092/65‑1 (MQ=255)
                                       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                            >  1:2452551/1‑68 (MQ=255)
                                       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                            >  1:17400595/1‑68 (MQ=255)
                                       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                            >  1:3005545/1‑68 (MQ=255)
                                       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                            >  1:5861016/1‑68 (MQ=255)
                                        gCGACCGGCGAGATGAAAGAACAGACTGATTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                           >  1:17316317/1‑68 (MQ=255)
                                        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                           >  1:18990744/1‑68 (MQ=255)
                                         cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                          >  1:26929541/1‑68 (MQ=255)
                                         cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                          >  1:28157843/1‑68 (MQ=255)
                                         cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGAGCTGTTCGGCAAACTGGCAGAAGTg                          >  1:36453245/1‑68 (MQ=255)
                                         agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                          >  1:12584852/2‑68 (MQ=255)
                                         agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                          >  1:19813154/2‑68 (MQ=255)
                                              ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag                     <  1:28805496/68‑1 (MQ=255)
                                                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                    <  1:23033549/67‑1 (MQ=255)
                                                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                    <  1:4510776/67‑1 (MQ=255)
                                                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                    <  1:38752761/67‑1 (MQ=255)
                                                     tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata               <  1:23222216/67‑1 (MQ=255)
                                                      gAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCTAACTGGCAGAAGTGgcgagc                    <  1:20212813/61‑1 (MQ=255)
                                                      gAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATc             <  1:27454402/68‑1 (MQ=255)
                                                        aaGAACAGACTGAATGGCACCGCGTTGTGATGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTg           <  1:33159092/68‑1 (MQ=255)
                                                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc          <  1:38897949/68‑1 (MQ=255)
                                                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc          <  1:37496219/68‑1 (MQ=255)
                                                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc          <  1:4223186/68‑1 (MQ=255)
                                                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc          <  1:25611785/68‑1 (MQ=255)
                                                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc          <  1:21561086/68‑1 (MQ=255)
                                                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc          <  1:20525561/64‑1 (MQ=255)
                                                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc          <  1:32079242/64‑1 (MQ=255)
                                                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc          <  1:8794430/64‑1 (MQ=255)
                                                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc          <  1:8494136/64‑1 (MQ=255)
                                                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc          <  1:14258309/64‑1 (MQ=255)
                                                          gagcagacagaATGGCATAGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc          <  1:28437266/64‑1 (MQ=255)
                                                          gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg         >  1:30861467/1‑68 (MQ=255)
                                                           agcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgc        >  1:18306041/3‑67 (MQ=255)
                                                            gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca       >  1:512860/2‑66 (MQ=255)
                                                            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa       >  1:6657546/1‑68 (MQ=255)
                                                             cAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaa      >  1:5070996/1‑68 (MQ=255)
                                                               gACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCaagg    >  1:10680949/1‑66 (MQ=255)
                                                                 cTGTATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt  <  1:22104829/68‑1 (MQ=255)
                                                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt  >  1:34545899/1‑68 (MQ=255)
                                                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt  <  1:36097512/68‑1 (MQ=255)
                                                                  |                                                                  
CAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGT  >  CP009273/4264145‑4264277

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