Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,264 A→G 18.2% Y71C (TAT→TGT)  ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2640AG18.2% 1853.4 / inf 719Y71C (TAT→TGT) ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (228/360);  new base G (33/98);  total (261/458)
Fisher's exact test for biased strand distribution p-value = 3.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 740 total aligned reads displayed.

GAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCC  >  CP009273/4264199‑4264331
                                                                 |                                                                   
   aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                                <  1:7339422/68‑1 (MQ=255)
   aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                                >  1:1724368/1‑68 (MQ=255)
   aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                                >  1:5835674/1‑68 (MQ=255)
   aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                                <  1:12184096/68‑1 (MQ=255)
    gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                                <  1:23204716/64‑1 (MQ=255)
    gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                                <  1:20175436/64‑1 (MQ=255)
    gagcagacagaATGGCATCGCGGGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                                <  1:14970024/64‑1 (MQ=255)
     aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                              >  1:16177542/1‑68 (MQ=255)
      gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                                                             >  1:16380081/2‑66 (MQ=255)
      aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                             >  1:20576663/1‑68 (MQ=255)
      aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                             >  1:13574585/1‑68 (MQ=255)
      aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                             >  1:15804879/1‑68 (MQ=255)
       cAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaa                                                            >  1:31668747/1‑68 (MQ=255)
        aGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                                                           >  1:19484313/1‑68 (MQ=255)
          aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                                         >  1:4122961/1‑68 (MQ=255)
           cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                                                           <  1:33994133/65‑1 (MQ=255)
           cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                        >  1:6825927/1‑68 (MQ=255)
               aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCaagg                                                          <  1:345335/62‑3 (MQ=255)
               aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                                      <  1:3022927/66‑4 (MQ=255)
               aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                                      <  1:6308930/66‑4 (MQ=255)
                tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                   >  1:13919600/1‑68 (MQ=255)
                tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                                   <  1:7028536/68‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:26718204/67‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:25368861/67‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:2238902/67‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:11740124/67‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:32294664/67‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:18070681/67‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:37134970/67‑1 (MQ=255)
                   cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                                 <  1:531246/67‑1 (MQ=255)
                    aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                               >  1:31997187/1‑68 (MQ=255)
                    aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                               >  1:32989097/1‑68 (MQ=255)
                    aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                               >  1:35528701/1‑68 (MQ=255)
                    aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                               >  1:35974285/1‑68 (MQ=255)
                    aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                               <  1:27557482/68‑1 (MQ=255)
                    aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                               >  1:22658941/1‑68 (MQ=255)
                    aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                               >  1:24403239/1‑68 (MQ=255)
                     tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                              <  1:509590/67‑3 (MQ=255)
                     tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                              <  1:9558709/67‑3 (MQ=255)
                         gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                           >  1:3661688/1‑67 (MQ=255)
                             tgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                      >  1:2505853/1‑68 (MQ=255)
                              gCTGTTCGGCAATCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATAGAAg                                     >  1:8749594/1‑68 (MQ=255)
                              gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                                     >  1:29309682/1‑68 (MQ=255)
                              gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                                      <  1:11952819/67‑1 (MQ=255)
                                tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                    <  1:1001959/67‑1 (MQ=255)
                                tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                    >  1:21563972/1‑67 (MQ=255)
                                tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                                    <  1:25848727/67‑1 (MQ=255)
                                   tCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCa                                 <  1:29798768/67‑1 (MQ=255)
                                   tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                 >  1:34850090/1‑67 (MQ=255)
                                   tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                 >  1:12512249/1‑67 (MQ=255)
                                   tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                 >  1:22862438/1‑67 (MQ=255)
                                   tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                                 >  1:36733363/1‑67 (MQ=255)
                                    cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAATGATCTCAGGTTTATATCGAAGGTCAGc                               <  1:34265136/68‑1 (MQ=255)
                                    cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                               <  1:21930646/68‑1 (MQ=255)
                                    cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc                               <  1:16881835/68‑1 (MQ=255)
                                     ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                              <  1:19724361/68‑1 (MQ=255)
                                     ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt                              <  1:19071472/68‑1 (MQ=255)
                                        aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                            <  1:35648817/67‑1 (MQ=255)
                                        aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                           >  1:5631883/1‑68 (MQ=255)
                                        aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                           >  1:24808974/1‑68 (MQ=255)
                                        aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGGTCTCAGGTTTATATCGAAGGTCAGCTGCg                           >  1:13139234/1‑68 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTTCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:27908424/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATTTCGAAGGTCAGCTGCGTAcc                       <  1:29117369/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:36924096/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:11320540/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:10993703/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:33547110/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:22091793/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:19875760/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:20777150/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:27164656/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:21245384/68‑1 (MQ=255)
                                            tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGCTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:26105946/68‑1 (MQ=255)
                                            tGGCAGAAGGGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:23916691/68‑1 (MQ=255)
                                             tggaGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTATATATCGAAGGTCAGCTGCGTAcc                       <  1:26417928/64‑1 (MQ=255)
                                             ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc                      <  1:29453100/68‑1 (MQ=255)
                                             ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc                       <  1:236066/67‑1 (MQ=255)
                                                tgAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                   <  1:10701023/67‑1 (MQ=255)
                                                ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt                    <  1:15367549/66‑1 (MQ=255)
                                                ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt                    <  1:28319323/66‑1 (MQ=255)
                                                ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                   <  1:4578718/67‑1 (MQ=255)
                                                ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                   <  1:28640924/67‑1 (MQ=255)
                                                  aaGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa                   >  1:166845/1‑66 (MQ=255)
                                                   aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                   <  1:6615095/65‑1 (MQ=255)
                                                   aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                   <  1:33380753/65‑1 (MQ=255)
                                                   aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                   <  1:31073796/65‑1 (MQ=255)
                                                   aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                   <  1:3610987/65‑1 (MQ=255)
                                                   aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                   <  1:32555251/65‑1 (MQ=255)
                                                   aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                 <  1:12776691/67‑3 (MQ=255)
                                                   aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                 <  1:17767226/67‑3 (MQ=255)
                                                   aGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTagc                 <  1:9150782/67‑3 (MQ=255)
                                                   aGTGCCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa                   <  1:4367469/65‑1 (MQ=255)
                                                    gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg               >  1:28794523/1‑68 (MQ=255)
                                                    gTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATg               >  1:7823850/1‑68 (MQ=255)
                                                        cgagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGAcc           >  1:35061486/1‑68 (MQ=255)
                                                         gagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCg          <  1:33184756/68‑1 (MQ=255)
                                                           gcgaATATCTGCGTAATGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGa         <  1:31408557/67‑1 (MQ=255)
                                                           gcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAATGTCAGCTGCGTACCCGTAAATGGACCGa         <  1:18374125/67‑1 (MQ=255)
                                                               atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt    <  1:14768358/68‑1 (MQ=255)
                                                               atatCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAAt    <  1:8285258/68‑1 (MQ=255)
                                                                 |                                                                   
GAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCGATCAATCC  >  CP009273/4264199‑4264331

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: