Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,205 A→G 8.6% E51E (GAA→GAG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2050AG 8.6% 1791.9 / 105.5 544E51E (GAA→GAGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (211/286);  new base G (32/15);  total (243/301)
Fisher's exact test for biased strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 576 total aligned reads displayed.

CAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTA  >  CP009273/4264138‑4264272
                                                                   |                                                                   
 aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                                    >  1:2532188/1‑68 (MQ=255)
    tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                                  <  1:12524865/67‑1 (MQ=255)
    tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                                 <  1:14783848/68‑1 (MQ=255)
       cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:4497522/67‑1 (MQ=255)
       cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                               <  1:17349223/67‑1 (MQ=255)
        aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGTCCGGCGAGATGAAAGAACAGACTGa                                                             <  1:26592216/68‑1 (MQ=255)
        aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:20678180/68‑1 (MQ=255)
        aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:3013834/68‑1 (MQ=255)
        aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                             <  1:10247858/68‑1 (MQ=255)
         aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:16227126/1‑68 (MQ=255)
         aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:26171671/1‑68 (MQ=255)
         aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:24483334/1‑68 (MQ=255)
         aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                            >  1:29236712/1‑68 (MQ=255)
           aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                          >  1:30421097/1‑68 (MQ=255)
           aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                                          >  1:30748957/1‑68 (MQ=255)
             tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGCGATGAAAGAACAGACTGAATGGc                                                        <  1:5022056/68‑1 (MQ=255)
             tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                                        <  1:10687820/68‑1 (MQ=255)
             tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCCAGATGAAAGAACAGACTGAATGGc                                                        <  1:28324822/68‑1 (MQ=255)
               cGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAc                                                      <  1:12107596/68‑1 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:37843650/1‑68 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:25644986/1‑68 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:13275816/1‑68 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:38204684/1‑68 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:38891970/1‑68 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:34960628/1‑68 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:5151179/1‑68 (MQ=255)
                gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:7056172/1‑68 (MQ=255)
                 cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:37316346/1‑67 (MQ=255)
                 cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:7750946/1‑67 (MQ=255)
                 cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                                     >  1:36832943/1‑67 (MQ=255)
                  tGGCTAGTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:15787038/68‑1 (MQ=255)
                  tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:28705492/68‑1 (MQ=255)
                  tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:34409498/68‑1 (MQ=255)
                  tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:5395130/68‑1 (MQ=255)
                  tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:16041056/68‑1 (MQ=255)
                  tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:39492629/68‑1 (MQ=255)
                  tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:20346903/68‑1 (MQ=255)
                  tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                                   <  1:22008869/68‑1 (MQ=255)
                     cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                                 <  1:18129649/67‑1 (MQ=255)
                       aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:13870915/68‑1 (MQ=255)
                       aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:1433833/68‑1 (MQ=255)
                       aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              <  1:15569552/68‑1 (MQ=255)
                          tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              >  1:24104619/1‑65 (MQ=255)
                          tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              >  1:31318927/1‑65 (MQ=255)
                          tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                              >  1:24875774/1‑65 (MQ=255)
                          tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                                           <  1:19373263/68‑1 (MQ=255)
                            cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:33252262/68‑1 (MQ=255)
                            cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:20776007/68‑1 (MQ=255)
                            cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:21360490/68‑1 (MQ=255)
                            cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:22625946/68‑1 (MQ=255)
                            cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                                         <  1:21576765/68‑1 (MQ=255)
                             gAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGtt                                        <  1:32863879/68‑1 (MQ=255)
                              aaTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTc                                       >  1:8921264/1‑68 (MQ=255)
                               aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                      >  1:33673883/1‑68 (MQ=255)
                               aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                                      >  1:4801264/1‑68 (MQ=255)
                                tCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTGCg                                      >  1:9651076/1‑67 (MQ=255)
                                 gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                    >  1:10078067/2‑68 (MQ=255)
                                 gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                    >  1:29848451/2‑68 (MQ=255)
                                 gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                    >  1:12830041/2‑68 (MQ=255)
                                 gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                                    >  1:11869239/2‑68 (MQ=255)
                                 ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                                    >  1:37439161/1‑68 (MQ=255)
                                      cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                                  >  1:6233722/1‑65 (MQ=255)
                                      cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                <  1:27994499/67‑3 (MQ=255)
                                      cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                                >  1:7443300/1‑65 (MQ=255)
                                       gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                              <  1:18832458/68‑1 (MQ=255)
                                       gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                              >  1:23062370/1‑68 (MQ=255)
                                       gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCt                               <  1:35285408/67‑1 (MQ=255)
                                           tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                           <  1:18589018/67‑1 (MQ=255)
                                              acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg                       >  1:3241962/4‑68 (MQ=255)
                                              aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                       >  1:5861016/1‑68 (MQ=255)
                                               gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                      >  1:35951238/1‑68 (MQ=255)
                                                cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                     >  1:23358042/1‑68 (MQ=255)
                                                cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                     >  1:26929541/1‑68 (MQ=255)
                                                cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                     >  1:17767142/1‑68 (MQ=255)
                                                cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGAGCTGTTCGGCAAACTGGCAGAAGTg                     >  1:36453245/1‑68 (MQ=255)
                                                  aCCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgc                   >  1:25921321/1‑68 (MQ=255)
                                                   cgggCGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                     <  1:1964432/63‑1 (MQ=255)
                                                    cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga                 <  1:34786586/68‑1 (MQ=255)
                                                    cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga                 <  1:20996704/68‑1 (MQ=255)
                                                     ggCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgag                <  1:872502/68‑1 (MQ=255)
                                                      gCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               >  1:22353499/1‑68 (MQ=255)
                                                       cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:38744759/67‑1 (MQ=255)
                                                       cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:8333674/67‑1 (MQ=255)
                                                       cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:11085379/67‑1 (MQ=255)
                                                       cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:13650858/67‑1 (MQ=255)
                                                       cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:9371400/67‑1 (MQ=255)
                                                       cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc               <  1:22157582/67‑1 (MQ=255)
                                                       cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGGGgcgagc               <  1:39041707/67‑1 (MQ=255)
                                                        gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga             <  1:29459729/68‑1 (MQ=255)
                                                             gAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATc        <  1:27454402/68‑1 (MQ=255)
                                                               taGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAAGTGGCAGAAGTGGCGAGCGAATATCt       <  1:24657667/66‑1 (MQ=255)
                                                                aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:16625589/68‑1 (MQ=255)
                                                                aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:4223186/68‑1 (MQ=255)
                                                                aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:6540848/68‑1 (MQ=255)
                                                                aGAACAGAATGAATGGCACCGCGTTGTGCTGTTCGGCTAACTGGCAGAAGTGGCGAGCGAATATCTGc     <  1:4074721/68‑1 (MQ=255)
                                                                 gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg    >  1:30911658/1‑68 (MQ=255)
                                                                  aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt   >  1:17190829/1‑68 (MQ=255)
                                                                  aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt   >  1:19910571/1‑68 (MQ=255)
                                                                   aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa  >  1:35668264/1‑68 (MQ=255)
                                                                   aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa  >  1:31798951/1‑68 (MQ=255)
                                                                   |                                                                   
CAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTA  >  CP009273/4264138‑4264272

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: