Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,190 C→G 6.0% T46T (ACC→ACG)  ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,1900CG 6.0% 2151.1 / 83.1 607T46T (ACC→ACG) ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (263/306);  new base G (25/12);  total (288/319)
Fisher's exact test for biased strand distribution p-value = 1.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 614 total aligned reads displayed.

GCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGA  >  CP009273/4264124‑4264257
                                                                  |                                                                   
 ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGCCCg                                                                    >  1:30221369/1‑67 (MQ=255)
  cAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCgg                                                                   <  1:1129875/67‑1 (MQ=255)
  cAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCgg                                                                   <  1:27397684/67‑1 (MQ=255)
      tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                              >  1:25776807/1‑68 (MQ=255)
       ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                             >  1:18467769/1‑68 (MQ=255)
       ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                             <  1:19564788/68‑1 (MQ=255)
       ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                             >  1:21285851/1‑68 (MQ=255)
       ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                             >  1:21721234/1‑68 (MQ=255)
       ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                             >  1:7366409/1‑68 (MQ=255)
        gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            >  1:28034193/1‑68 (MQ=255)
         tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            <  1:15693713/67‑1 (MQ=255)
         tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            <  1:22634870/67‑1 (MQ=255)
         tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            <  1:10785982/67‑1 (MQ=255)
         tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                            <  1:31241881/67‑1 (MQ=255)
          ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                          <  1:585870/68‑1 (MQ=255)
             gcAGTTGCCAACATTACGCTGGCTACTTCCGACTCCTGGCGTGATAAAGCGACCGGGGAGATGAAAg                                                        >  1:13211234/1‑67 (MQ=255)
             gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCTGCGAGATGAAAg                                                        <  1:16431450/67‑1 (MQ=255)
             gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                        <  1:13429066/67‑1 (MQ=255)
             gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                        <  1:6289302/67‑1 (MQ=255)
             gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                        <  1:8104208/67‑1 (MQ=255)
             gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                        <  1:8382636/67‑1 (MQ=255)
             gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                       <  1:26249741/68‑1 (MQ=255)
               aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                     >  1:31051987/1‑68 (MQ=255)
               aGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAAc                                                     >  1:31924911/1‑68 (MQ=255)
                  tGCCAACATTACGCTGGCTACTTCCGAATCCTGGGGTGATAAAGGTACCGGCGAGATGAAAGAACAGa                                                  <  1:28675600/68‑1 (MQ=255)
                  tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                                   <  1:32474696/67‑1 (MQ=255)
                  tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                  <  1:4090337/68‑1 (MQ=255)
                  tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                                  <  1:16932537/68‑1 (MQ=255)
                     cAAGATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACAGGCTAGATGAAAGATCAGACt                                                <  1:20610080/67‑1 (MQ=255)
                     cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                <  1:25296371/67‑1 (MQ=255)
                     cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                <  1:27651239/67‑1 (MQ=255)
                     cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                                <  1:4862241/67‑1 (MQ=255)
                      aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAGCAGACTGa                                              <  1:24479179/68‑1 (MQ=255)
                      aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:11737707/68‑1 (MQ=255)
                      aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:28442823/68‑1 (MQ=255)
                      aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:21005310/68‑1 (MQ=255)
                      aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                              <  1:30519216/68‑1 (MQ=255)
                       aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGGGAGATGAAAGAACAGACTGaa                                             >  1:17503371/1‑68 (MQ=255)
                       aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:653236/1‑68 (MQ=255)
                       aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:32775172/1‑68 (MQ=255)
                       aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:22383902/1‑68 (MQ=255)
                       aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                             >  1:28920261/1‑68 (MQ=255)
                         aTTACGCTGGCTACTTCCGAATCCTGGCGTGGTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                           >  1:34281471/1‑68 (MQ=255)
                         aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                           >  1:1884194/1‑68 (MQ=255)
                         aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATg                                           >  1:19068085/1‑68 (MQ=255)
                         aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACAGAATg                                           >  1:25472802/1‑68 (MQ=255)
                         aTTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGACATGAAAGAACAGACTGAATg                                           >  1:31775637/1‑68 (MQ=255)
                              gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:1046036/1‑68 (MQ=255)
                              gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:8984102/1‑68 (MQ=255)
                              gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:14896950/1‑68 (MQ=255)
                              gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:36869333/1‑68 (MQ=255)
                               cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:35400772/1‑67 (MQ=255)
                               cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:23239558/1‑67 (MQ=255)
                               cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:35667475/1‑67 (MQ=255)
                               cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:36451273/1‑67 (MQ=255)
                               cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                      >  1:32998456/1‑67 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:10021534/68‑1 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:5156641/68‑1 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:32899657/68‑1 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:34001122/68‑1 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:18178612/68‑1 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:18003373/68‑1 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACCGAATGGCACcgc                                    <  1:2682081/68‑1 (MQ=255)
                                tGGCTACTTCCGAATCCTGGCGGGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                    <  1:13801628/68‑1 (MQ=255)
                                   cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                                 >  1:26605960/1‑68 (MQ=255)
                                   cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                  <  1:3484995/67‑1 (MQ=255)
                                   cTACATCCGAATCCTGGCGTGATAAAGCGAGCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGt                                  <  1:4594532/67‑1 (MQ=255)
                                     aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                               <  1:10230131/68‑1 (MQ=255)
                                        tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                               >  1:12271434/1‑65 (MQ=255)
                                        tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCt                            <  1:6817644/68‑1 (MQ=255)
                                          cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                          <  1:11499729/68‑1 (MQ=255)
                                          cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                          <  1:15914262/68‑1 (MQ=255)
                                          cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                          <  1:12637222/68‑1 (MQ=255)
                                          cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                          <  1:19852452/68‑1 (MQ=255)
                                          cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                          <  1:31467502/68‑1 (MQ=255)
                                             aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                       >  1:23967926/1‑68 (MQ=255)
                                               gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                     >  1:28236892/2‑68 (MQ=255)
                                               gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                     >  1:3795342/2‑68 (MQ=255)
                                               gcTGGCGTGACAAACAGACGGGGGAGATGCGGCAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                     >  1:34615631/2‑68 (MQ=255)
                                               ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                     >  1:9454697/1‑68 (MQ=255)
                                                  ggCGTGATAAATCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                   >  1:22907186/1‑67 (MQ=255)
                                                  ggCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaag                  <  1:14493171/68‑2 (MQ=255)
                                                    cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCaa                   >  1:33245251/1‑65 (MQ=255)
                                                    cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                 <  1:22100162/67‑3 (MQ=255)
                                                    cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc                 <  1:721689/67‑3 (MQ=255)
                                                     gTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTg               <  1:29655184/68‑1 (MQ=255)
                                                       gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg              <  1:25071957/67‑1 (MQ=255)
                                                       gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg              <  1:5452401/67‑1 (MQ=255)
                                                       gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg              <  1:15413461/67‑1 (MQ=255)
                                                        aTAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCa            >  1:2928947/1‑68 (MQ=255)
                                                          aaaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGTAGa          >  1:17286433/1‑68 (MQ=255)
                                                           aaGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGaa         >  1:29968255/1‑68 (MQ=255)
                                                            acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg        <  1:20328025/65‑1 (MQ=255)
                                                            acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg        >  1:33557242/4‑68 (MQ=255)
                                                            acagACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAg        >  1:929101/4‑68 (MQ=255)
                                                              cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg      >  1:14118499/1‑68 (MQ=255)
                                                               gACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCTCGCGGAAGTg      >  1:36312999/1‑67 (MQ=255)
                                                               gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTgg     <  1:16679022/68‑1 (MQ=255)
                                                               gACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg      >  1:11057965/1‑67 (MQ=255)
                                                                  cGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGgcga  <  1:7917733/68‑1 (MQ=255)
                                                                  |                                                                   
GCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTGGCGA  >  CP009273/4264124‑4264257

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: ‑