Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,193 C→G 7.9% G47G (GGC→GGG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,1930CG 7.9% 2033.9 / 117.3 582G47G (GGC→GGGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (267/269);  new base G (32/14);  total (299/283)
Fisher's exact test for biased strand distribution p-value = 1.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 615 total aligned reads displayed.

AAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCG  >  CP009273/4264127‑4264260
                                                                  |                                                                   
aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                                    >  1:20329542/1‑68 (MQ=255)
  aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                                  >  1:27596945/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                 >  1:9591661/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                 >  1:3603877/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                 >  1:32305265/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                 >  1:26706860/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                 >  1:18819376/1‑68 (MQ=255)
   tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                                 >  1:16757823/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                                >  1:5904042/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                                >  1:29306215/1‑68 (MQ=255)
    ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                                >  1:6962707/1‑68 (MQ=255)
    ggtgcCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                                >  1:30375594/1‑68 (MQ=255)
     gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                               >  1:39423344/1‑68 (MQ=255)
       ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                             >  1:27219856/1‑68 (MQ=255)
         cgcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                           >  1:23385144/1‑68 (MQ=255)
         cgcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                           >  1:18595065/1‑68 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                           <  1:7360653/67‑1 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                          <  1:2220352/68‑1 (MQ=255)
          gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                          <  1:758466/68‑1 (MQ=255)
            aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGaa                                                         <  1:32613348/67‑1 (MQ=255)
            aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                        >  1:13675886/1‑68 (MQ=255)
            aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                        >  1:38770743/1‑68 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                      <  1:34790226/67‑1 (MQ=255)
               tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                     <  1:25234338/68‑1 (MQ=255)
                  cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                   <  1:39631023/67‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                 <  1:4441417/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                 <  1:26046017/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                 <  1:32866075/68‑1 (MQ=255)
                   aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                                 <  1:10247858/68‑1 (MQ=255)
                    aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                                >  1:29864249/1‑68 (MQ=255)
                      aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                              <  1:36847117/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGTAAGAACAGACTGAATGGc                                            <  1:30147606/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                            <  1:3386525/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                            <  1:11963531/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                            <  1:33551752/68‑1 (MQ=255)
                        tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCCAGATGAAAGAACAGACTGAATGGc                                            <  1:28324822/68‑1 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:5151179/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:8980834/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:15811337/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:38946712/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:16203939/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:31186031/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:2092179/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:21201134/1‑68 (MQ=255)
                           gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:19124188/1‑68 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:37316346/1‑67 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:25594848/1‑67 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:10545012/1‑67 (MQ=255)
                            cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                         >  1:8725799/1‑67 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                       <  1:11001690/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                       <  1:30713897/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                       <  1:28705492/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                       <  1:12999671/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                       <  1:16041056/68‑1 (MQ=255)
                             tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                       <  1:4772044/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                  <  1:13870915/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                  <  1:15569552/68‑1 (MQ=255)
                                  aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                  <  1:1433833/68‑1 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                  >  1:20743325/1‑65 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                  >  1:24104619/1‑65 (MQ=255)
                                     tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                                  >  1:22357720/1‑65 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                             <  1:22625946/68‑1 (MQ=255)
                                       cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                             <  1:21576722/68‑1 (MQ=255)
                                        gAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGtt                            <  1:18225230/68‑1 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTc                           <  1:1213723/67‑1 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                          >  1:37544240/1‑68 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                          >  1:4199158/1‑68 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                          >  1:22015666/1‑68 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                          >  1:35097760/1‑68 (MQ=255)
                                          aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGATGTTCg                          >  1:19492814/1‑68 (MQ=255)
                                           tCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCgg                         <  1:1330965/68‑1 (MQ=255)
                                            gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                        >  1:10078067/2‑68 (MQ=255)
                                            ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                        >  1:2303874/1‑68 (MQ=255)
                                               ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaa                      >  1:6267616/1‑67 (MQ=255)
                                                 cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                    <  1:27994499/67‑3 (MQ=255)
                                                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                  <  1:17120015/68‑1 (MQ=255)
                                                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                  >  1:5001793/1‑68 (MQ=255)
                                                  gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg                  <  1:35601089/68‑1 (MQ=255)
                                                  gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCt                   <  1:35285408/67‑1 (MQ=255)
                                                         acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg           >  1:3197339/4‑68 (MQ=255)
                                                         acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg           <  1:5503092/65‑1 (MQ=255)
                                                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg           >  1:2452551/1‑68 (MQ=255)
                                                         aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg           >  1:17400595/1‑68 (MQ=255)
                                                           cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg         >  1:2654307/1‑68 (MQ=255)
                                                           cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg         >  1:17767142/1‑68 (MQ=255)
                                                           cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGAGCTGTTCGGCAAACTGGCAGAAGTg         >  1:36453245/1‑68 (MQ=255)
                                                           agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg         >  1:19813154/2‑68 (MQ=255)
                                                           agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg         >  1:12584852/2‑68 (MQ=255)
                                                           agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg         <  1:38515328/67‑1 (MQ=255)
                                                            gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt          <  1:2481810/66‑1 (MQ=255)
                                                            gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg        <  1:14757797/68‑1 (MQ=255)
                                                            gACCGGCGAGATGAAAGAACAGACTCAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg        <  1:37907988/68‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcg  >  1:20588530/1‑68 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc   <  1:31586514/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc   <  1:29454679/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc   <  1:7916989/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc   <  1:8333674/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc   <  1:25756272/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc   <  1:9328039/67‑1 (MQ=255)
                                                                  cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc   <  1:11562010/67‑1 (MQ=255)
                                                                  |                                                                   
AAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCG  >  CP009273/4264127‑4264260

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