Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,289 T→C 8.6% 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 8.6% 2307.7 / 137.0 699Y79Y (TAT→TACssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (160/477);  new base C (19/41);  total (180/519)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 729 total aligned reads displayed.

GCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAA  >  CP009273/4264222‑4264355
                                                                   |                                                                  
gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAGTATCTGCGTAAAGGTTCTCAGGTTtat                                                                    <  1:33850674/68‑1 (MQ=255)
  gTTGTGCTGTTCGGCACACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                                                                  <  1:14307417/68‑1 (MQ=255)
  gTTGTGCTGTTCGGCAAGCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                   <  1:4426663/67‑1 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                   <  1:13873916/67‑1 (MQ=255)
  gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                                                   <  1:11962571/67‑1 (MQ=255)
    tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc                                                                 <  1:17183168/67‑1 (MQ=255)
    ggtgCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATc                                                                 <  1:2284043/66‑1 (MQ=255)
       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                                             >  1:3829422/1‑68 (MQ=255)
       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                                             >  1:30296128/1‑68 (MQ=255)
       gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                                             >  1:601136/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                                              <  1:18654898/67‑1 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                             >  1:1650261/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                             >  1:32683342/1‑68 (MQ=255)
       gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                                                             >  1:32536051/1‑68 (MQ=255)
        cTGTTCGGCAAACTGGCAGAAGTCGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                            >  1:19068767/1‑68 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                            <  1:24409250/67‑1 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                            <  1:27351290/67‑1 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                            <  1:13297127/67‑1 (MQ=255)
         tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                                            <  1:12469510/67‑1 (MQ=255)
            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                                         >  1:21873797/1‑67 (MQ=255)
            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                                         >  1:22389183/1‑67 (MQ=255)
            tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGAGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                                         >  1:27138207/1‑67 (MQ=255)
             cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                                                       <  1:6288502/68‑1 (MQ=255)
             cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                                                       <  1:10434287/68‑1 (MQ=255)
             cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                                                       <  1:11536555/68‑1 (MQ=255)
              ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                                                      <  1:11655264/68‑1 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                                                    <  1:8024950/67‑1 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                   >  1:28360385/1‑68 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                   >  1:13146200/1‑68 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                   >  1:31061395/1‑68 (MQ=255)
                 aaaCTGGCAGAAGTGGCGAGCGAAGATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                                                   >  1:24937153/1‑68 (MQ=255)
                     tGGCAGAAGTGGCGAGCTAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:3476814/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:23585551/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:21336786/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:27984203/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:4189101/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:914772/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:8850588/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:12672491/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:649307/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:1535451/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:15977608/68‑1 (MQ=255)
                     tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                                               <  1:2030731/68‑1 (MQ=255)
                      ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                                              <  1:7433871/68‑1 (MQ=255)
                      ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                                              <  1:29588697/68‑1 (MQ=255)
                         ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                           <  1:9663687/67‑1 (MQ=255)
                           aaGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                                           >  1:25694155/1‑66 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                           <  1:33523632/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                           <  1:15510202/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                           <  1:17411863/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                           <  1:7757200/65‑1 (MQ=255)
                            aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                                           <  1:34296051/65‑1 (MQ=255)
                            aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                                         <  1:31550759/67‑3 (MQ=255)
                             gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg                                       >  1:4696304/1‑68 (MQ=255)
                               ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg                                      <  1:5149997/67‑1 (MQ=255)
                                 cgagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGAcc                                   >  1:16841359/1‑68 (MQ=255)
                                 cgagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGAcc                                   >  1:26453730/1‑68 (MQ=255)
                                  gagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCg                                  <  1:23203549/68‑1 (MQ=255)
                                    gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                                 <  1:25505974/67‑1 (MQ=255)
                                    gcgaATATCTGCGTAAAGGCTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa                                 <  1:10006388/67‑1 (MQ=255)
                                        atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGGACCCGTAAATGGACCGATCAAt                            <  1:29183339/68‑1 (MQ=255)
                                            cTGCGTAACGGTTCTCAGGTTTATATCGAAGGGCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                        <  1:2112275/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                        <  1:5910449/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                        >  1:8492957/1‑68 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                        <  1:12989373/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGCACCGATCAATCCgg                        <  1:6979548/68‑1 (MQ=255)
                                            cTGCGTAAAGGTTCTCAGGTTTATATCGAAGGGCAGCTGCGTACCCGTAAATGGACCGATCAATCCgg                        <  1:32859803/68‑1 (MQ=255)
                                                gTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAg                    >  1:9422297/1‑68 (MQ=255)
                                                  aaaGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGa                  <  1:30863620/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:19312043/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:1828081/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:15001789/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:14768681/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:24601761/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:25138390/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:13653654/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:30599179/68‑1 (MQ=255)
                                                    aGGTTCTCAGGTTTATATCGAAGGGCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                <  1:29910766/68‑1 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                >  1:11205356/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                >  1:10700440/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                >  1:26938969/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                >  1:20740925/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                >  1:21674677/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATc                >  1:33081637/1‑67 (MQ=255)
                                                     ggTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAACCGGTCAGGATc                >  1:5610935/1‑67 (MQ=255)
                                                       ttCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCTATCAATCCGGTCAGGATCg               <  1:32036345/66‑1 (MQ=255)
                                                       ttCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCg               <  1:25258875/66‑1 (MQ=255)
                                                       ttCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCg               <  1:14432968/66‑1 (MQ=255)
                                                        tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCt             <  1:19681947/67‑1 (MQ=255)
                                                        tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCt             <  1:17954209/67‑1 (MQ=255)
                                                        tctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCt             <  1:22109937/67‑1 (MQ=255)
                                                         ctcAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTac           <  1:26694814/68‑1 (MQ=255)
                                                            aGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc        >  1:26232889/1‑68 (MQ=255)
                                                             ggTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc        >  1:8676125/1‑67 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc        <  1:31758536/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc        <  1:20017433/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc        <  1:27593012/65‑1 (MQ=255)
                                                               tttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTAcacc        <  1:10112236/65‑1 (MQ=255)
                                                                ttATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACcaca     <  1:24203349/67‑1 (MQ=255)
                                                                  atatCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGaa  <  1:28842223/68‑1 (MQ=255)
                                                                   |                                                                  
GCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCCGGTCAGGATCGCTACACCACAGAA  >  CP009273/4264222‑4264355

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: