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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,2470AG11.3% 2078.1 / 179.1 665A65A (GCA→GCGssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (245/344);  new base G (35/40);  total (280/385)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 671 total aligned reads displayed.

ATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGT  >  CP009273/4264180‑4264313
                                                                   |                                                                  
aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                                                    >  1:9041397/1‑68 (MQ=255)
 tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa                                                                    <  1:225408/67‑1 (MQ=255)
     gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                               >  1:17979388/1‑68 (MQ=255)
     gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt                                                               >  1:8752722/1‑68 (MQ=255)
      cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg                                                              >  1:24761700/1‑68 (MQ=255)
      agACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg                                                              >  1:4172381/2‑68 (MQ=255)
       gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg                                                             <  1:23418080/68‑1 (MQ=255)
        aCCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgc                                                            >  1:27105916/1‑68 (MQ=255)
             cGAGATGAAAGAACAGACTGACTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                        <  1:2213873/67‑1 (MQ=255)
             cGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagc                                                        <  1:28124259/67‑1 (MQ=255)
              gagaTGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                      <  1:14659905/68‑1 (MQ=255)
              gagaTGAAAGAACAGACTGAATGGCACCGCGGTGGGCGGTTCGGCAAACTGGCAGAAGTGgcgagcga                                                      <  1:3849406/68‑1 (MQ=255)
                  tGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAata                                                   <  1:3042345/67‑1 (MQ=255)
                      acaaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGACGTGGCGAGCGAATATCTGc                                              <  1:15096320/66‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTTTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                              <  1:33553493/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                              <  1:8144889/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                              <  1:14028823/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                              <  1:18272950/68‑1 (MQ=255)
                      aGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGc                                              <  1:6481126/68‑1 (MQ=255)
                       gAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCg                                             >  1:19304692/1‑68 (MQ=255)
                        agcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgc                                            >  1:34068704/3‑67 (MQ=255)
                        aaCAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGt                                            >  1:25101994/1‑68 (MQ=255)
                         gcagacagaATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCgca                                           >  1:15391437/2‑66 (MQ=255)
                         aCAGACTAAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTa                                           >  1:8583624/1‑68 (MQ=255)
                            gACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAg                                        >  1:25489784/1‑68 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTaaa                                         <  1:4642682/65‑1 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                      <  1:17312729/68‑1 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                      >  1:2444476/1‑68 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                      >  1:14524333/1‑68 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                      <  1:9734360/68‑1 (MQ=255)
                              cTGAATGGCACCGCGTTGTGCTGTTCGACAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGt                                      <  1:19824220/68‑1 (MQ=255)
                                  aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                    <  1:5506094/66‑4 (MQ=255)
                                  aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAgggcgc                                    <  1:20119311/66‑4 (MQ=255)
                                  aTGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAggg                                       <  1:13439537/63‑1 (MQ=255)
                                  aTGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTtctc                                  >  1:17062920/1‑68 (MQ=255)
                                   tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                 >  1:18823486/1‑68 (MQ=255)
                                   tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                 >  1:34476589/1‑68 (MQ=255)
                                   tGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCa                                 >  1:32665129/1‑68 (MQ=255)
                                     gCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               >  1:28734029/1‑68 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               <  1:25718566/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               <  1:10216490/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               <  1:29037630/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               <  1:33060493/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               <  1:23841957/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               <  1:7874789/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAgg                               <  1:18506735/67‑1 (MQ=255)
                                      cACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATATGCGTAAAGGTTCTCAgg                               <  1:20837301/67‑1 (MQ=255)
                                       atcgcgTGGTGCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGt                              <  1:21880582/65‑1 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:31081387/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:11888535/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:29478700/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:13078421/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:33868490/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:1607042/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:17604489/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:19631412/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:20474405/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:5273090/1‑68 (MQ=255)
                                       aCCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGtt                             >  1:32823354/1‑68 (MQ=255)
                                          gcgTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtat                          <  1:30717998/68‑1 (MQ=255)
                                            gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTtata                         <  1:31771666/67‑1 (MQ=255)
                                            gTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTatat                        <  1:7187158/68‑1 (MQ=255)
                                                 gCTGTTCGGCAAGCTCGCGGAAGTGGCAGGTGAATGTCTGCGCAAGGGCGCGCAGGTCTACATCGAAg                   >  1:26438741/1‑68 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTTCGTAAAGGTGCTCAGGTTTATATCGaa                    <  1:23240929/67‑1 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGaa                    <  1:6691948/67‑1 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                   >  1:29874904/1‑68 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                   >  1:7602072/1‑68 (MQ=255)
                                                 gCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAg                   >  1:4848550/1‑68 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGTGCGATTATCTGCGTATTGGTTCTCAGGTTTATATCGAAgg                  <  1:9861489/67‑1 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                  <  1:31402465/67‑1 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                  <  1:10002020/67‑1 (MQ=255)
                                                   tGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAgg                  >  1:11466251/1‑67 (MQ=255)
                                                      tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCTAAGGTCa               >  1:5026362/1‑67 (MQ=255)
                                                      tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa               >  1:21275835/1‑67 (MQ=255)
                                                      tCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCa               >  1:14284475/1‑67 (MQ=255)
                                                       cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc             <  1:10433656/68‑1 (MQ=255)
                                                       cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc             <  1:4051412/68‑1 (MQ=255)
                                                       cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc             <  1:1065215/68‑1 (MQ=255)
                                                       cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc             <  1:28458305/68‑1 (MQ=255)
                                                       cGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGc             <  1:6958867/68‑1 (MQ=255)
                                                           aaaCTGGCAGAAGTTGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg         >  1:6145679/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGCAGGTCAGCTGCg         >  1:1645009/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc          <  1:27795841/67‑1 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGc          <  1:4386907/67‑1 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg         >  1:32009673/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg         >  1:9489318/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg         >  1:9299837/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg         >  1:12179967/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg         >  1:16666486/1‑68 (MQ=255)
                                                           aaaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCg         >  1:26141305/1‑68 (MQ=255)
                                                            aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt        <  1:21057718/68‑1 (MQ=255)
                                                            aaCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt        <  1:26827032/68‑1 (MQ=255)
                                                             aCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGt        >  1:12703285/1‑67 (MQ=255)
                                                             aCTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGT‑‑TCAGGTTTATATCGAAGGTCAGCTGCGTAcc     <  1:15583786/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc     <  1:5865285/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc     <  1:11696305/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc     <  1:11004879/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc     <  1:21644225/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc     <  1:22942002/68‑1 (MQ=255)
                                                               tGGCAGAAGTGGCGAGCGAATATCCGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTAcc     <  1:27575517/68‑1 (MQ=255)
                                                                   |                                                                  
ATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGAGCGAATATCTGCGTAAAGGTTCTCAGGTTTATATCGAAGGTCAGCTGCGTACCCGT  >  CP009273/4264180‑4264313

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: