Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,264,186 C→A 5.5% A45E (GCG→GAG) ‡ ssb → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,264,1860CA 5.5% 2128.7 / 71.3 592A45E (GCG→GAG) ‡ssbsingle‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (297/262);  new base A (23/10);  total (320/272)
Fisher's exact test for biased strand distribution p-value = 7.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 613 total aligned reads displayed.

TACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTG  >  CP009273/4264119‑4264253
                                                                   |                                                                   
 aCATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCg                                                                    >  1:11991341/1‑68 (MQ=255)
 aCATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCg                                                                    >  1:17078618/1‑68 (MQ=255)
      ccAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCg                                                                >  1:11631891/1‑67 (MQ=255)
        aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                             >  1:23766016/1‑68 (MQ=255)
        aaaTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCg                                                             >  1:22930889/1‑68 (MQ=255)
          aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                           >  1:31738034/1‑68 (MQ=255)
          aTGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgag                                                           >  1:27596945/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                          >  1:9591661/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                          <  1:3336862/68‑1 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                          >  1:32754415/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                          >  1:32463252/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                          >  1:26706860/1‑68 (MQ=255)
           tggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCgaga                                                          >  1:22473053/1‑68 (MQ=255)
            ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                         >  1:12166247/1‑68 (MQ=255)
            ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                         >  1:13126106/1‑68 (MQ=255)
            ggtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGAt                                                         >  1:14076055/1‑68 (MQ=255)
             gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                        >  1:39314784/1‑68 (MQ=255)
             gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                        >  1:32208888/1‑68 (MQ=255)
             gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                        >  1:28649930/1‑68 (MQ=255)
             gtggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                        >  1:21980547/1‑68 (MQ=255)
              tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                        <  1:5051534/67‑1 (MQ=255)
              tggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATg                                                        <  1:20877014/67‑1 (MQ=255)
               ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                      >  1:27219856/1‑68 (MQ=255)
               ggCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGaa                                                      >  1:21048690/1‑68 (MQ=255)
                  gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                    <  1:8944853/67‑1 (MQ=255)
                  gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                    <  1:25413066/67‑1 (MQ=255)
                  gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAg                                                    <  1:22137225/67‑1 (MQ=255)
                  gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                   <  1:15555756/68‑1 (MQ=255)
                  gcAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGa                                                   >  1:38620603/1‑68 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                               <  1:20724048/67‑1 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                               <  1:19648052/67‑1 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                               <  1:34790226/67‑1 (MQ=255)
                       tGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGa                                              <  1:2500967/68‑1 (MQ=255)
                         ccAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAg                                               >  1:6177100/1‑65 (MQ=255)
                          cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                            <  1:8107808/67‑1 (MQ=255)
                          cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACt                                            <  1:17349223/67‑1 (MQ=255)
                          cAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTg                                           >  1:37222847/1‑68 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:26046017/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:4818816/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:6433443/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:23208137/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:4441417/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:3013834/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:10247858/68‑1 (MQ=255)
                           aaCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGa                                          <  1:14852642/68‑1 (MQ=255)
                            aCATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGaa                                         >  1:29236712/1‑68 (MQ=255)
                                tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                     <  1:16706287/68‑1 (MQ=255)
                                tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                     <  1:5197201/68‑1 (MQ=255)
                                tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                     <  1:19959670/68‑1 (MQ=255)
                                tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGc                                     <  1:19426916/68‑1 (MQ=255)
                                tACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCCAGATGAAAGAACAGACTGAATGGc                                     <  1:28324822/68‑1 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:24434683/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:13754062/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:7121466/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:16203939/1‑68 (MQ=255)
                                   gCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:22284550/1‑68 (MQ=255)
                                    cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:8725799/1‑67 (MQ=255)
                                    cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:10545012/1‑67 (MQ=255)
                                    cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:7750946/1‑67 (MQ=255)
                                    cTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCAcc                                  >  1:29385484/1‑67 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                <  1:30625252/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                <  1:20346903/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                <  1:34409498/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                <  1:18714531/68‑1 (MQ=255)
                                     tGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACcgc                                <  1:11001690/68‑1 (MQ=255)
                                        cTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGtt                             >  1:39518520/1‑68 (MQ=255)
                                          aCTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                           <  1:9377131/68‑1 (MQ=255)
                                             tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                           >  1:31318927/1‑65 (MQ=255)
                                             tCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTtgt                           >  1:35568303/1‑65 (MQ=255)
                                               cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                      <  1:21576722/68‑1 (MQ=255)
                                               cGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGt                      >  1:30290558/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                   >  1:4199158/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                   >  1:33673883/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                   >  1:30912126/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                   >  1:35097760/1‑68 (MQ=255)
                                                  aTCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCg                   >  1:4041357/1‑68 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                 >  1:38864747/2‑68 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                 >  1:10078067/2‑68 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                 >  1:36431113/2‑68 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                 >  1:17632634/2‑68 (MQ=255)
                                                    gcTGGCGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGc                 >  1:29848451/2‑68 (MQ=255)
                                                    ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                 >  1:37439161/1‑68 (MQ=255)
                                                    ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                 >  1:36830940/1‑68 (MQ=255)
                                                    ccTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGc                 >  1:14487468/1‑68 (MQ=255)
                                                         cGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAAc             <  1:13358568/67‑1 (MQ=255)
                                                         cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc             >  1:15704166/1‑65 (MQ=255)
                                                         cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc             >  1:7443300/1‑65 (MQ=255)
                                                         cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc             <  1:12701015/67‑3 (MQ=255)
                                                         cGTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCaagc             <  1:27994499/67‑3 (MQ=255)
                                                          gTGACAAACAGACGGGGGAGATGCGGGAGCAGACAGAATGGCATCGCGTGGTGCTGTTCGGCAAGCt            <  1:35285408/67‑1 (MQ=255)
                                                            gATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTgg          <  1:24584936/67‑1 (MQ=255)
                                                              tAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAg       >  1:18574536/1‑68 (MQ=255)
                                                                 aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg    >  1:25715445/1‑68 (MQ=255)
                                                                 aGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAg    <  1:37887921/68‑1 (MQ=255)
                                                                  gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt   >  1:18990744/1‑68 (MQ=255)
                                                                  gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt   >  1:35951238/1‑68 (MQ=255)
                                                                  gCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGt   >  1:15974354/1‑68 (MQ=255)
                                                                   cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg  >  1:2654307/1‑68 (MQ=255)
                                                                   cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg  >  1:5646989/1‑68 (MQ=255)
                                                                   cGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTg  >  1:17874336/1‑68 (MQ=255)
                                                                   |                                                                   
TACATGCCAAATGGTGGCGCAGTTGCCAACATTACGCTGGCTACTTCCGAATCCTGGCGTGATAAAGCGACCGGCGAGATGAAAGAACAGACTGAATGGCACCGCGTTGTGCTGTTCGGCAAACTGGCAGAAGTG  >  CP009273/4264119‑4264253

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