Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,190 C→G 5.4% T46T (ACC→ACG ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,1900CG 5.4% 1963.0 / 65.4 549T46T (ACC→ACGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (241/277);  new base G (20/10);  total (262/287)
Fisher's exact test for biased strand distribution p-value = 3.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 552 total aligned reads displayed.

TGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGA  >  CP009273/4264123‑4264257
                                                                   |                                                                   
tGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGAcc                                                                     >  1:7134279/1‑68 (MQ=255)
  ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                                    >  1:17817653/1‑67 (MQ=255)
  ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                                    >  1:754427/1‑67 (MQ=255)
  ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                                    >  1:26480609/1‑67 (MQ=255)
  ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                                    >  1:23893042/1‑67 (MQ=255)
  ccAAATGGTGGCGCAGTTGCCAACATCACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                                    >  1:12055379/1‑67 (MQ=255)
   cAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCgg                                                                   <  1:2695594/67‑1 (MQ=255)
   cAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCgg                                                                   <  1:34102231/67‑1 (MQ=255)
   cAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCgg                                                                   <  1:12353437/67‑1 (MQ=255)
      aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                               >  1:17248272/1‑68 (MQ=255)
      aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                               >  1:4742570/1‑68 (MQ=255)
       tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                              >  1:2790288/1‑68 (MQ=255)
       tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                              >  1:25388138/1‑68 (MQ=255)
       tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                              <  1:29195008/68‑1 (MQ=255)
        ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                             >  1:22167401/1‑68 (MQ=255)
          tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            <  1:25981898/67‑1 (MQ=255)
          tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            <  1:24055730/67‑1 (MQ=255)
          tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            <  1:20858619/67‑1 (MQ=255)
              gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                       <  1:6584422/68‑1 (MQ=255)
                aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGaa                                                      <  1:14102237/67‑1 (MQ=255)
                aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                     >  1:3144243/1‑68 (MQ=255)
                   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                   <  1:30946748/67‑1 (MQ=255)
                   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                   <  1:19732524/67‑1 (MQ=255)
                   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                  <  1:28279706/68‑1 (MQ=255)
                   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                  <  1:15188599/68‑1 (MQ=255)
                   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                  <  1:24576218/68‑1 (MQ=255)
                   tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                  <  1:5636833/68‑1 (MQ=255)
                      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                <  1:8796786/67‑1 (MQ=255)
                      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                <  1:4978154/67‑1 (MQ=255)
                      cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                <  1:15046217/67‑1 (MQ=255)
                       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:18153485/68‑1 (MQ=255)
                       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:27626045/68‑1 (MQ=255)
                       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:6408458/68‑1 (MQ=255)
                       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:3215744/68‑1 (MQ=255)
                       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:1305571/68‑1 (MQ=255)
                       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:16768839/68‑1 (MQ=255)
                       aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:21654964/68‑1 (MQ=255)
                        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:21853274/1‑68 (MQ=255)
                        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:12739382/1‑68 (MQ=255)
                        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:24441308/1‑68 (MQ=255)
                        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:14423330/1‑68 (MQ=255)
                        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:14744699/1‑68 (MQ=255)
                        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:15430166/1‑68 (MQ=255)
                        aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:16151457/1‑68 (MQ=255)
                          aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                           >  1:4063573/1‑68 (MQ=255)
                          aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                           >  1:18509282/1‑68 (MQ=255)
                          aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                           >  1:11434968/1‑68 (MQ=255)
                           ttACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                           >  1:22649804/1‑67 (MQ=255)
                            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                         <  1:23366150/68‑1 (MQ=255)
                            tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                         <  1:30302924/68‑1 (MQ=255)
                             aCGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                         >  1:29679329/1‑67 (MQ=255)
                             aCGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                         >  1:24963308/1‑67 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:17575484/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:28585765/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:33352520/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:32770600/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:3164128/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:32593078/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:10178410/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:15976818/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:21042348/1‑68 (MQ=255)
                               gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:16240022/1‑68 (MQ=255)
                                cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACGGAATGGcccc                                      >  1:11034454/1‑67 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTTATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:737522/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:33569568/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:9610235/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:10665092/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:27294169/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:1883155/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:19106736/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:31216866/68‑1 (MQ=255)
                                 tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:20047783/68‑1 (MQ=255)
                                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                  <  1:1830407/67‑1 (MQ=255)
                                    cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                  <  1:34509290/67‑1 (MQ=255)
                                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                               <  1:2024318/68‑1 (MQ=255)
                                      aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                               <  1:1326603/68‑1 (MQ=255)
                                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                               >  1:9074065/1‑65 (MQ=255)
                                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                            <  1:9205831/68‑1 (MQ=255)
                                         tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                            <  1:25770994/68‑1 (MQ=255)
                                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                          <  1:16950254/68‑1 (MQ=255)
                                           cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                          <  1:10002105/68‑1 (MQ=255)
                                              aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                       >  1:10494487/1‑68 (MQ=255)
                                              aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                       <  1:15755489/68‑1 (MQ=255)
                                              aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTCAATGGCACCGCGTTGTGCTGTTCg                       >  1:16445024/1‑68 (MQ=255)
                                                gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                     >  1:26327816/2‑68 (MQ=255)
                                                gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                     >  1:8291167/2‑68 (MQ=255)
                                                gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                     >  1:28900968/2‑68 (MQ=255)
                                                     cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                   >  1:1372324/1‑65 (MQ=255)
                                                     cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                   >  1:33206964/1‑65 (MQ=255)
                                                     cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                 <  1:3915373/67‑3 (MQ=255)
                                                     cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                 <  1:20519404/67‑3 (MQ=255)
                                                      gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCtc               <  1:18641296/68‑2 (MQ=255)
                                                        gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg              <  1:34737214/67‑1 (MQ=255)
                                                         aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa            >  1:23749361/1‑68 (MQ=255)
                                                          tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa            <  1:30468968/67‑1 (MQ=255)
                                                          tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAg           >  1:6015530/1‑68 (MQ=255)
                                                           aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGa          >  1:11710066/1‑68 (MQ=255)
                                                             aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg        >  1:21326380/1‑68 (MQ=255)
                                                              gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt       >  1:5464898/1‑68 (MQ=255)
                                                                gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg     <  1:29426409/68‑1 (MQ=255)
                                                                   |                                                                   
TGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGA  >  CP009273/4264123‑4264257

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: