Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,244 G→C 11.2% L64L (CTG→CTC ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2440GC11.2% 2140.6 / 206.8 664L64L (CTG→CTCssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (248/341);  new base C (34/40);  total (283/381)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 678 total aligned reads displayed.

GTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCC  >  CP009273/4264177‑4264311
                                                                   |                                                                   
  gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg                                                                    <  1:16281940/67‑1 (MQ=255)
   aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                                                  >  1:9041397/1‑68 (MQ=255)
   aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                                                  >  1:5199707/1‑68 (MQ=255)
    tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                                                  <  1:225408/67‑1 (MQ=255)
       aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg                                                              >  1:24685439/1‑68 (MQ=255)
        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                             >  1:20104843/1‑68 (MQ=255)
        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                             >  1:8752722/1‑68 (MQ=255)
        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                             >  1:33012975/1‑68 (MQ=255)
        gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                             >  1:17979388/1‑68 (MQ=255)
         cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                                            >  1:27653759/1‑68 (MQ=255)
          gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                             <  1:3680420/66‑1 (MQ=255)
          gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                                                           <  1:814164/68‑1 (MQ=255)
                cGAGATGAAAGGACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:11276634/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGACTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:2213873/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:32723880/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:4836533/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:34027599/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:13088647/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:12737091/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:6463746/67‑1 (MQ=255)
                cGAGATGAAAGAACAGACTGAAGGGCACCGCGGTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                      <  1:33469518/67‑1 (MQ=255)
                 gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                    <  1:14659905/68‑1 (MQ=255)
                     tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata                                                 >  1:4170057/1‑67 (MQ=255)
                        aaGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTg                                             <  1:17237329/68‑1 (MQ=255)
                         acaaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGACGTGGCGAGCGAATATCTGc                                            <  1:15096320/66‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                            <  1:30622576/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                            <  1:23356389/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                            <  1:18272950/68‑1 (MQ=255)
                         aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                            <  1:31275561/68‑1 (MQ=255)
                          gagcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGc                                            <  1:21248249/64‑1 (MQ=255)
                          gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                           >  1:19304692/1‑68 (MQ=255)
                            gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                                         >  1:15391437/2‑66 (MQ=255)
                            gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                                         >  1:25392652/2‑66 (MQ=255)
                            aCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                         >  1:6061983/1‑68 (MQ=255)
                                aCTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAgg                                     >  1:717901/1‑68 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                    <  1:3086799/68‑1 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                    >  1:26398747/1‑68 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                    <  1:31610454/68‑1 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                    <  1:14080602/68‑1 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                    >  1:7834346/1‑68 (MQ=255)
                                 cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                    <  1:9734360/68‑1 (MQ=255)
                                   gAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtc                                  <  1:15615807/68‑1 (MQ=255)
                                     aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtctc                                >  1:17062920/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGATTATCTGCGTAAAGGTTCTCa                               >  1:26084057/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                               >  1:207905/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                               >  1:18823486/1‑68 (MQ=255)
                                      tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                               >  1:289937/1‑68 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCTTACTGGGAGAAGTGGGGATGGGTTATCTGCGTAAAGGTTCTCAgg                             <  1:27455807/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:4005821/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:5697212/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:2706097/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:10216490/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:2677001/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:24296382/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:18506735/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:19619721/67‑1 (MQ=255)
                                         cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                             <  1:23225749/67‑1 (MQ=255)
                                          atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt                            <  1:21880582/65‑1 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                           >  1:8519465/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                           >  1:6873799/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                           >  1:4600325/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                           >  1:1607042/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                           >  1:13510253/1‑68 (MQ=255)
                                          aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                           >  1:12249558/1‑68 (MQ=255)
                                           ccGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGttt                          <  1:28328194/68‑1 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                       >  1:17561142/1‑67 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                       >  1:32509916/1‑67 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                       <  1:3195315/67‑1 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                       >  1:22893013/1‑67 (MQ=255)
                                               gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                      <  1:28051481/68‑1 (MQ=255)
                                                 tgtgCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATc                     <  1:29383648/67‑1 (MQ=255)
                                                 ggtgCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCg                    >  1:6657137/2‑68 (MQ=255)
                                                    gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                 >  1:14350297/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                 >  1:29644293/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                 >  1:7602072/1‑68 (MQ=255)
                                                    gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                 >  1:4367330/1‑68 (MQ=255)
                                                      tGTTCGGCAAACTGGCAGAAGTGGCGTGCGATTATCTGCGTATTGGTTCTCAGGTTTATATCGAAgg                <  1:9861489/67‑1 (MQ=255)
                                                      tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                <  1:10002020/67‑1 (MQ=255)
                                                      tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                <  1:24977862/67‑1 (MQ=255)
                                                       gTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTc              >  1:5468097/1‑68 (MQ=255)
                                                         tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa             >  1:21591240/1‑67 (MQ=255)
                                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc           <  1:25709625/68‑1 (MQ=255)
                                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc           <  1:10433656/68‑1 (MQ=255)
                                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc           <  1:14556394/68‑1 (MQ=255)
                                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc           <  1:1065215/68‑1 (MQ=255)
                                                          cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc           <  1:11948463/68‑1 (MQ=255)
                                                           ggCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt          <  1:21052589/68‑1 (MQ=255)
                                                            gCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCt          >  1:9277764/1‑67 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc        <  1:3433574/67‑1 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg       >  1:30275677/1‑68 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg       >  1:9489318/1‑68 (MQ=255)
                                                              aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg       >  1:16666486/1‑68 (MQ=255)
                                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc   <  1:23809156/68‑1 (MQ=255)
                                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc   <  1:2865075/68‑1 (MQ=255)
                                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc   <  1:11702657/68‑1 (MQ=255)
                                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc   <  1:11483203/68‑1 (MQ=255)
                                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc   <  1:20472181/68‑1 (MQ=255)
                                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc   <  1:4980722/68‑1 (MQ=255)
                                                                  tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc   <  1:33805640/68‑1 (MQ=255)
                                                                   ggCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAccc  <  1:19818614/68‑1 (MQ=255)
                                                                   |                                                                   
GTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCC  >  CP009273/4264177‑4264311

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: