Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,289 T→C 11.1% 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,2890TC11.1% 2210.7 / 186.8 710Y79Y (TAT→TACssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (188/443);  new base C (21/58);  total (209/501)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 745 total aligned reads displayed.

GCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAAG  >  CP009273/4264222‑4264356
                                                                   |                                                                   
 cgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                    >  1:21792194/1‑68 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                    >  1:10414136/1‑67 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                    <  1:878420/67‑1 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                    <  1:29077367/67‑1 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                    <  1:16763100/67‑1 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                                                                   <  1:10846159/68‑1 (MQ=255)
     gtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGa                                                                >  1:9505113/1‑68 (MQ=255)
      tgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                               <  1:28684410/68‑1 (MQ=255)
       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                                              >  1:29309682/1‑68 (MQ=255)
       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                                              >  1:22916679/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                               <  1:4846927/67‑1 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                               <  1:34573709/67‑1 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                               <  1:11952819/67‑1 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                              >  1:14887397/1‑68 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                             <  1:1001959/67‑1 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATACCGAAgg                                                             <  1:34856977/67‑1 (MQ=255)
            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                                          >  1:32486826/1‑67 (MQ=255)
             cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                                                        <  1:29931439/68‑1 (MQ=255)
             cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                                                        <  1:21519290/68‑1 (MQ=255)
              ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                                                       <  1:274068/68‑1 (MQ=255)
                cAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCttc                                                     >  1:9160271/1‑66 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAGGGTCAGCTGCg                                                    >  1:31651821/1‑68 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                                     <  1:25049042/67‑1 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                                     <  1:23230536/67‑1 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                    >  1:18474809/1‑68 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                    >  1:36279258/1‑68 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                    >  1:36637114/1‑68 (MQ=255)
                  aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt                                                   <  1:5497327/68‑1 (MQ=255)
                   aCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt                                                   >  1:25847787/1‑67 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:3409914/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:22091793/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:7211208/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTAAGGTTTATATCGAAGGTCAGCTGCGTAcc                                                <  1:12177132/68‑1 (MQ=255)
                      ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                                               <  1:1154618/68‑1 (MQ=255)
                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt                                             <  1:12697970/66‑1 (MQ=255)
                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                            <  1:35594523/67‑1 (MQ=255)
                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                            <  1:3887579/67‑1 (MQ=255)
                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                            <  1:11297001/67‑1 (MQ=255)
                         ggAAGTGGCAGGTGAATGGCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTTGCACCCGTa                                            <  1:36980362/67‑1 (MQ=255)
                          gAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTaa                                           >  1:26200852/1‑68 (MQ=255)
                          gAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                            >  1:35960873/1‑67 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:15933305/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:6830011/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:30599167/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:6615095/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:6059585/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:18203885/65‑1 (MQ=255)
                            aGTGGCGAGAGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                            <  1:35336346/65‑1 (MQ=255)
                            aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                                          <  1:16601067/67‑3 (MQ=255)
                            aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                                          <  1:20344507/67‑3 (MQ=255)
                             gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg                                        >  1:32253179/1‑68 (MQ=255)
                             gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg                                        >  1:2848892/1‑68 (MQ=255)
                                 cgagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGAcc                                    >  1:19609244/1‑68 (MQ=255)
                                    gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                                  <  1:33438247/67‑1 (MQ=255)
                                    gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                                  <  1:15504690/67‑1 (MQ=255)
                                    gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                                  <  1:19916549/67‑1 (MQ=255)
                                    gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                                  <  1:19304338/67‑1 (MQ=255)
                                        atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt                             <  1:4300518/68‑1 (MQ=255)
                                        atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt                             <  1:13911903/68‑1 (MQ=255)
                                          atCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGCTCAATcc                           >  1:34272713/1‑68 (MQ=255)
                                           tCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCg                          >  1:17272620/1‑68 (MQ=255)
                                           tCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCg                          >  1:16321427/1‑68 (MQ=255)
                                           tCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCg                          >  1:32775397/1‑68 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:5738341/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:26770851/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                         <  1:125030/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:24869086/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:21158197/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:3189965/68‑1 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGGCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                   <  1:20763186/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:4536210/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:26605495/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:2426964/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:8590619/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:23996760/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:17445218/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:2966055/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:23740233/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:14652453/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:3721221/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:31043828/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:33666614/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:22624062/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 <  1:22441995/68‑1 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 >  1:33332594/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 >  1:5169792/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 >  1:13472410/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                 >  1:17790753/1‑67 (MQ=255)
                                                       ttCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCg                <  1:7611317/66‑1 (MQ=255)
                                                        tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCt              <  1:19769504/67‑1 (MQ=255)
                                                        tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCt              <  1:821096/67‑1 (MQ=255)
                                                         ctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTac            <  1:6520123/68‑1 (MQ=255)
                                                            aGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcccc         >  1:35027321/1‑68 (MQ=255)
                                                             ggTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         >  1:31222890/1‑67 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:27803772/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:7458233/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:18831971/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTAACCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc         <  1:13075253/65‑1 (MQ=255)
                                                                ttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCTCTACACCACAg     >  1:25367488/1‑68 (MQ=255)
                                                                ttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAg     >  1:1528824/1‑68 (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: