Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,256 G→A 18.9% A68A (GCG→GCA ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2560GA18.9% 1617.0 / inf 694A68A (GCG→GCAssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (222/340);  new base A (33/98);  total (255/439)
Fisher's exact test for biased strand distribution p-value = 2.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 761 total aligned reads displayed.

CCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCG  >  CP009273/4264189‑4264323
                                                                   |                                                                   
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:10000779/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:25306465/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:23581191/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:33517494/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:22814850/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:20734894/67‑1 (MQ=255)
    cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGAGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                                  <  1:26887547/67‑1 (MQ=255)
     gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                <  1:10523618/68‑1 (MQ=255)
     gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                                <  1:27984693/68‑1 (MQ=255)
             aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        >  1:5835674/1‑68 (MQ=255)
             aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                                        <  1:26028532/68‑1 (MQ=255)
              gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                        <  1:2893224/64‑1 (MQ=255)
              gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                                        <  1:31111888/64‑1 (MQ=255)
               aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                                      >  1:15033156/1‑68 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:13574585/1‑68 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:1619048/1‑68 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:1698068/1‑68 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:6807260/1‑68 (MQ=255)
                aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                                     >  1:20576663/1‑68 (MQ=255)
                    aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                                 >  1:8026473/1‑68 (MQ=255)
                    aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                                 >  1:22068023/1‑68 (MQ=255)
                     cTGAATTGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:15383295/68‑1 (MQ=255)
                     cTGAATGGCACGGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:35120426/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                                                   <  1:33994133/65‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:20125744/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:18219301/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:31468493/68‑1 (MQ=255)
                     cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                >  1:6907615/1‑68 (MQ=255)
                     cTGAATGGCACCGCGGTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                                <  1:34608977/68‑1 (MQ=255)
                          tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                           >  1:36033114/1‑68 (MQ=255)
                          tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                           >  1:12965381/1‑68 (MQ=255)
                          tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                           >  1:36354408/1‑68 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:19465131/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:27332982/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:4570089/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:37062554/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:30692799/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:9116763/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:17462970/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:22269026/67‑1 (MQ=255)
                             cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                                         <  1:6028060/67‑1 (MQ=255)
                              atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt                                        <  1:32120458/65‑1 (MQ=255)
                              atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt                                        <  1:17866889/65‑1 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:32989097/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:36684820/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:34895677/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:29478167/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:19827896/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:18766355/1‑68 (MQ=255)
                              aCCGCGTTGTGCTGTTCAGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                                       >  1:1094593/1‑68 (MQ=255)
                               tcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGtct                                      <  1:9558709/67‑3 (MQ=255)
                                 gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                    <  1:28368179/68‑1 (MQ=255)
                                 gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                    <  1:5066733/68‑1 (MQ=255)
                                 gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                                    <  1:12407490/68‑1 (MQ=255)
                                   gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                   >  1:28661600/1‑67 (MQ=255)
                                   gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                                   >  1:23815000/1‑67 (MQ=255)
                                     tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc                                 <  1:11842288/67‑1 (MQ=255)
                                        gCTGTTCGGCAATCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATAGAAg                             >  1:8749594/1‑68 (MQ=255)
                                        gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                             >  1:34315332/1‑68 (MQ=255)
                                        gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                             >  1:21459388/1‑68 (MQ=255)
                                        gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                              <  1:24362150/67‑1 (MQ=255)
                                          tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                            <  1:19717475/67‑1 (MQ=255)
                                             tCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCa                         <  1:29798768/67‑1 (MQ=255)
                                             tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         >  1:30132108/1‑67 (MQ=255)
                                             tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa                         <  1:3990220/67‑1 (MQ=255)
                                              cGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGc                       <  1:15834300/68‑1 (MQ=255)
                                               ggCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCt                      <  1:14826246/68‑1 (MQ=255)
                                               ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGAGTAAAGGTTCTCAGGTTTATATCGAAGGtaa                         >  1:13956095/1‑63 (MQ=255)
                                                  aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc                    <  1:25049042/67‑1 (MQ=255)
                                                  aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:7165209/1‑68 (MQ=255)
                                                  aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg                   >  1:24808974/1‑68 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATTTCGAAGGTCAGCTGCGTAcc               <  1:29117369/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:16502386/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:16709117/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:32390668/68‑1 (MQ=255)
                                                      tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc               <  1:26514966/68‑1 (MQ=255)
                                                       ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc              <  1:1154618/68‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGt            <  1:12697970/66‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:7558272/67‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:11297001/67‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:35594523/67‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:25680379/67‑1 (MQ=255)
                                                          ggAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           <  1:32596016/67‑1 (MQ=255)
                                                           gAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTaa          >  1:26200852/1‑68 (MQ=255)
                                                           gAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAGGTCAGCTTCGCACCCGTa           >  1:35917869/1‑67 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCTAAGGTCAGCTGCGTACCCGTa           <  1:21750509/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:33380753/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:27697926/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:29358027/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:31760154/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:23013094/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:5580913/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:22514816/65‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt        <  1:6571567/68‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt        <  1:5250017/68‑1 (MQ=255)
                                                             aGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAAt        >  1:21075091/1‑68 (MQ=255)
                                                             aGTGCCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTa           <  1:4367469/65‑1 (MQ=255)
                                                                ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg      <  1:1213752/67‑1 (MQ=255)
                                                                ggCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATgg      <  1:367890/67‑1 (MQ=255)
                                                                   gagcgaATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCg  <  1:33184756/68‑1 (MQ=255)
                                                                   |                                                                   
CCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGTAAATGGACCG  >  CP009273/4264189‑4264323

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: