Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,289 T→C 6.9% Y79Y (TAT→TAC ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2890TC 6.9% 2429.1 / 105.6 707Y79Y (TAT→TACssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (176/482);  new base C (23/26);  total (199/508)
Fisher's exact test for biased strand distribution p-value = 4.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

Only 100 of 743 total aligned reads displayed.

GCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAG  >  CP009273/4264222‑4264356
                                                                   |                                                                   
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                                                                   >  1:19571099/1‑68 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                                                                   <  1:18494125/68‑1 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                                                                   <  1:7187158/68‑1 (MQ=255)
    tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc                                                                  <  1:22129035/67‑1 (MQ=255)
    ggtgCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATc                                                                  <  1:13057531/66‑1 (MQ=255)
    ggtgCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCg                                                                 >  1:6657137/2‑68 (MQ=255)
       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                                              >  1:13499722/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                              >  1:9901657/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                              >  1:7974772/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                              >  1:7602072/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                              >  1:27131102/1‑68 (MQ=255)
        cTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                             <  1:31680788/68‑1 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                             <  1:10002020/67‑1 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                             >  1:28944172/1‑67 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                             <  1:30547736/67‑1 (MQ=255)
          gTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTc                                                           >  1:16820284/1‑68 (MQ=255)
             cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                                                        <  1:4093522/68‑1 (MQ=255)
             cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                                                        <  1:27709730/68‑1 (MQ=255)
              ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                                                       <  1:14642274/68‑1 (MQ=255)
              ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                                                       <  1:340227/68‑1 (MQ=255)
               gCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                                                       >  1:32698587/1‑67 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                                     <  1:4386907/67‑1 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                    >  1:26141305/1‑68 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                    >  1:17808032/1‑68 (MQ=255)
                  aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt                                                   <  1:21057718/68‑1 (MQ=255)
                   aCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt                                                   >  1:12703285/1‑67 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:32214760/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:5296676/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:8162970/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:21888739/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:103376/68‑1 (MQ=255)
                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                            <  1:4653611/67‑1 (MQ=255)
                          gAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTaa                                           >  1:3510301/1‑68 (MQ=255)
                          gAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                            >  1:16600340/1‑67 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:20558385/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:19342040/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:9835481/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:236073/65‑1 (MQ=255)
                            aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagct                                         <  1:8909308/68‑4 (MQ=255)
                            aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                                          <  1:34312140/67‑3 (MQ=255)
                             gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg                                        >  1:14400609/1‑68 (MQ=255)
                             gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg                                        >  1:17367090/1‑68 (MQ=255)
                               ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg                                       <  1:14116344/67‑1 (MQ=255)
                                    gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                                  <  1:1552892/67‑1 (MQ=255)
                                        atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCaa                              <  1:2991899/67‑1 (MQ=255)
                                        atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt                             <  1:29541368/68‑1 (MQ=255)
                                          atCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAtgc                           >  1:15094848/1‑66 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:6287775/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:20884469/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:25554162/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:2414131/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:23877097/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:31131094/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:21408955/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:12340452/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:10895779/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:18214693/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:15874251/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:11315274/68‑1 (MQ=255)
                                             tGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTaaa                                          >  1:30967331/1‑50 (MQ=255)
                                                gTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCa                      <  1:31966698/67‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:13914046/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:11840/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:28362376/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:11495420/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:2324754/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:4983864/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   >  1:5718951/1‑68 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:12109152/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:3960742/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:10546654/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:19868501/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:24764569/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:19988316/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:21221362/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:18991408/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:24543528/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:33199576/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:1625545/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:14808376/68‑1 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 >  1:2676574/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 >  1:25701921/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 >  1:6776928/1‑67 (MQ=255)
                                                       ttCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCg                <  1:9423358/66‑1 (MQ=255)
                                                       ttCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCg                <  1:29245941/66‑1 (MQ=255)
                                                        tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCt              <  1:21698418/67‑1 (MQ=255)
                                                        tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCt              <  1:22552112/67‑1 (MQ=255)
                                                         ctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTac            <  1:27400605/68‑1 (MQ=255)
                                                         ctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTac            <  1:9230237/68‑1 (MQ=255)
                                                           cAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTacac          <  1:25136559/68‑1 (MQ=255)
                                                             ggTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTacaac         >  1:8485589/1‑65 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:28626484/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:30345380/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:5941587/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:31383642/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACcaca      >  1:22810084/1‑68 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACcaca      >  1:11035051/1‑68 (MQ=255)
                                                                ttttatCGAAGGTCAGCTGCGTACCCGTAAATGGACCTATCAATCCGTTCAGGATCGCTACACcaca      <  1:5141568/64‑1 (MQ=255)
                                                                ttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACcaca      <  1:33582352/67‑1 (MQ=255)
                                                                   tatCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAg  <  1:8073272/68‑1 (MQ=255)
                                                                   |                                                                   
GCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAG  >  CP009273/4264222‑4264356

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: