Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,220 C→T 14.2% H56H (CAC→CAT ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2200CT14.2% 1663.7 / 232.9 623H56H (CAC→CATssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (224/308);  new base T (46/42);  total (270/353)
Fisher's exact test for biased strand distribution p-value = 8.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 687 total aligned reads displayed.

CGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTT  >  CP009273/4264153‑4264286
                                                                   |                                                                  
cGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAc                                                                    <  1:12107596/68‑1 (MQ=255)
 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACCGACTGAATGGCAcc                                                                   >  1:15957136/1‑68 (MQ=255)
 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:38891970/1‑68 (MQ=255)
 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:25836654/1‑68 (MQ=255)
 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:8980834/1‑68 (MQ=255)
 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:7857073/1‑68 (MQ=255)
 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:12984497/1‑68 (MQ=255)
 gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:13783529/1‑68 (MQ=255)
  cTGGCTACTTCCGAATCCTGGCGTGGTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:26600610/1‑67 (MQ=255)
  cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:25594848/1‑67 (MQ=255)
  cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:29385484/1‑67 (MQ=255)
  cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:31842745/1‑67 (MQ=255)
  cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                                   >  1:10545012/1‑67 (MQ=255)
   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                                 <  1:39111412/68‑1 (MQ=255)
   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                                 <  1:22465396/68‑1 (MQ=255)
   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                                 <  1:28853430/68‑1 (MQ=255)
   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                                 <  1:30649473/68‑1 (MQ=255)
   tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                                 <  1:11001690/68‑1 (MQ=255)
      cTACTTCCGTATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACATACTTAATGGCACCGCGt                                                               <  1:28058027/67‑1 (MQ=255)
      cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                              >  1:39518520/1‑68 (MQ=255)
       tACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                                              <  1:620189/67‑1 (MQ=255)
        aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                            <  1:15569552/68‑1 (MQ=255)
        aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                            <  1:34922212/68‑1 (MQ=255)
        aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                            <  1:18640835/68‑1 (MQ=255)
           tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                                            >  1:31318927/1‑65 (MQ=255)
             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                       <  1:35188106/68‑1 (MQ=255)
             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                       <  1:21576765/68‑1 (MQ=255)
             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                       <  1:21360490/68‑1 (MQ=255)
             cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                                       <  1:17849076/68‑1 (MQ=255)
              gAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGtt                                                      <  1:32863879/68‑1 (MQ=255)
                aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                                    >  1:37544240/1‑68 (MQ=255)
                 tCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                                    >  1:14228075/1‑67 (MQ=255)
                 agcTGGCGTGTCAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCgg                                                   <  1:14519398/66‑1 (MQ=255)
                  gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                                  >  1:10078067/2‑68 (MQ=255)
                  ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                                  >  1:36830940/1‑68 (MQ=255)
                     ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaa                                                >  1:38684716/1‑67 (MQ=255)
                       cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                                >  1:30159347/1‑65 (MQ=255)
                       cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAAc                                              <  1:24461592/67‑1 (MQ=255)
                       cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                              <  1:33780846/67‑3 (MQ=255)
                        gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                            <  1:3789926/68‑1 (MQ=255)
                        gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                                            <  1:35601089/68‑1 (MQ=255)
                          gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                                           <  1:24584936/67‑1 (MQ=255)
                            tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                         <  1:37306320/67‑1 (MQ=255)
                            tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                         <  1:9675666/67‑1 (MQ=255)
                               acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                                     <  1:35945244/65‑1 (MQ=255)
                               acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                                     >  1:26038428/4‑68 (MQ=255)
                               aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                     >  1:17400595/1‑68 (MQ=255)
                                 cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                   >  1:2481152/1‑68 (MQ=255)
                                 cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                   >  1:38323452/1‑68 (MQ=255)
                                 cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                   >  1:26929541/1‑68 (MQ=255)
                                  gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                                  <  1:14757797/68‑1 (MQ=255)
                                    cgggCGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                                   <  1:1964432/63‑1 (MQ=255)
                                     cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga                               <  1:20996704/68‑1 (MQ=255)
                                     cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga                               <  1:20591527/68‑1 (MQ=255)
                                      ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag                              <  1:450888/68‑1 (MQ=255)
                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcg                            >  1:16198420/1‑68 (MQ=255)
                                        cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                             <  1:23033549/67‑1 (MQ=255)
                                         gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                           <  1:13610547/68‑1 (MQ=255)
                                                 aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                   <  1:6540848/68‑1 (MQ=255)
                                                  gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                   <  1:8494136/64‑1 (MQ=255)
                                                  gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                  >  1:30911658/1‑68 (MQ=255)
                                                   agcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgc                 >  1:18306041/3‑67 (MQ=255)
                                                   aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                 >  1:17190829/1‑68 (MQ=255)
                                                    gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                >  1:23207229/2‑66 (MQ=255)
                                                    aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                >  1:5117639/1‑68 (MQ=255)
                                                    aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                >  1:19195049/1‑68 (MQ=255)
                                                    aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                >  1:27794576/1‑68 (MQ=255)
                                                     cAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaa               >  1:24108487/1‑68 (MQ=255)
                                                        aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg            >  1:28451600/1‑68 (MQ=255)
                                                         cagAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCaag              <  1:23597922/63‑2 (MQ=255)
                                                         cTGTATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           <  1:22104829/68‑1 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCACACTGGCAGAAGTGGCGCGCGAATATCTGCGTAAAGGt           >  1:1268527/1‑68 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           >  1:10243439/1‑68 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           <  1:494907/68‑1 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           <  1:413044/68‑1 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           >  1:34545899/1‑68 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           <  1:16953024/68‑1 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           <  1:20099492/68‑1 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           <  1:32251743/68‑1 (MQ=255)
                                                         cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt           >  1:21978524/1‑68 (MQ=255)
                                                              tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa      >  1:5680956/1‑68 (MQ=255)
                                                              tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa      >  1:3414565/1‑68 (MQ=255)
                                                               ggCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAg     >  1:679535/1‑68 (MQ=255)
                                                                gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    >  1:25061887/1‑68 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:29918040/67‑1 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:20393797/67‑1 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:536199/67‑1 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:7452903/67‑1 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:15763208/67‑1 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:8775388/67‑1 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:12349369/67‑1 (MQ=255)
                                                                 cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:11571143/67‑1 (MQ=255)
                                                                 cACCGCGTGGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAATGTTCTCAgg    <  1:7022258/67‑1 (MQ=255)
                                                                 cACCGCGGTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg    <  1:36341320/67‑1 (MQ=255)
                                                                  aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt  >  1:3350318/1‑68 (MQ=255)
                                                                  aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt  >  1:18881347/1‑68 (MQ=255)
                                                                  aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt  >  1:22295806/1‑68 (MQ=255)
                                                                  aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt  >  1:2736932/1‑68 (MQ=255)
                                                                  aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt  >  1:811109/1‑68 (MQ=255)
                                                                  aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt  >  1:9437658/1‑68 (MQ=255)
                                                                   |                                                                  
CGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTT  >  CP009273/4264153‑4264286

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: