Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 692,043 T→A 5.2% noncoding (30/75 nt) glnV ← tRNA‑Gln

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273692,0430TA 5.2% 1606.7 / 20.1 524noncoding (30/75 nt)glnVtRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base T (316/173);  new base A (18/16);  total (335/189)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 546 total aligned reads displayed.

AAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACA  >  CP009273/691976‑692110
                                                                   |                                                                   
aaGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAt                                                                     <  1:17616686/68‑1 (MQ=255)
  ggTAAGCGTCTTGGATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:49958355/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:26972891/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:40778464/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:4246214/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:39506706/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:38797224/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:45135680/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:3707152/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:46544986/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:35953411/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:33308307/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:1476332/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:50110971/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:6251678/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:35808953/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:35339381/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:9112293/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACATCGGAATGCCGGAATCAGAATcc                                                                   >  1:2830119/1‑68 (MQ=255)
  ggTAAGCGTCTTGAATAAATTGGCTCTTGGACGAGGAGTCGAACCTCGGATTGACTGAATGAGAATcc                                                                   >  1:20518044/1‑68 (MQ=16)
   gTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCg                                                                  >  1:31584328/1‑68 (MQ=255)
   gTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAACCg                                                                  >  1:26652813/1‑68 (MQ=255)
    tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                 >  1:44385597/1‑68 (MQ=255)
    tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                 >  1:30190719/1‑68 (MQ=255)
    tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGAATCAGAATACgg                                                                 >  1:214746/1‑68 (MQ=255)
    tAAGCGTCTTGAATAAATTGGCTGGGGTACCAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                 >  1:7168480/1‑68 (MQ=255)
     aaGGGTCTTTAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                                <  1:22914993/68‑1 (MQ=255)
      aGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                                >  1:15319549/1‑67 (MQ=255)
        cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCATAATCCGGTg                                                               >  1:34361684/1‑66 (MQ=255)
        cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                               >  1:23697028/1‑66 (MQ=255)
        cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                               >  1:32421516/1‑66 (MQ=255)
        cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                               >  1:12142900/1‑66 (MQ=255)
        cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCCAACCTCGGAAAGCCGGAATCAGAATCCGCTg                                                               >  1:45732429/1‑66 (MQ=17)
         ttCTTTAATAAATTGGCTGGGGTACGAGGATTCGAACCTGGGAATGCCGGTATCAAAAACCGGTGcc                                                             <  1:48702812/66‑1 (MQ=14)
         gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCt                                                            >  1:18328217/1‑68 (MQ=255)
         gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCt                                                            >  1:24898878/1‑68 (MQ=255)
         gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAACCGGTGCCt                                                            >  1:40125026/1‑68 (MQ=255)
          tCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCtt                                                           >  1:1007817/1‑68 (MQ=255)
             ttaaTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                        <  1:6741008/66‑1 (MQ=255)
             ttaaTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                        <  1:10529613/66‑1 (MQ=255)
               aaTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                                      <  1:24434714/68‑1 (MQ=255)
                aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                     >  1:966839/1‑68 (MQ=40)
                aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                     <  1:29029699/68‑1 (MQ=40)
                aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                     >  1:17210573/1‑68 (MQ=40)
                  aaaTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTg                                                   >  1:22175641/1‑68 (MQ=37)
                   atttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTgg                                                  <  1:32847812/66‑1 (MQ=14)
                    aTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 <  1:50706362/68‑1 (MQ=33)
                    aTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 <  1:11812043/68‑1 (MQ=33)
                     ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                >  1:16444614/1‑68 (MQ=31)
                     ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                >  1:19456910/1‑68 (MQ=31)
                     ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                >  1:17748261/1‑68 (MQ=31)
                     ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCATGTCg                                                >  1:33619640/1‑68 (MQ=255)
                                aCGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaaca                                     <  1:7852466/68‑1 (MQ=31)
                                   aGGATTCGAAGCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:50611322/1‑68 (MQ=255)
                                   aGGATTCGAACCTCGGAATGCGGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACaaaa                                  >  1:5818165/1‑67 (MQ=255)
                                   aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacattt                                  >  1:33920344/1‑65 (MQ=14)
                                   aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacattt                                  >  1:47904107/1‑65 (MQ=14)
                                   aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacagtg                                  >  1:49143796/1‑65 (MQ=14)
                                   aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:2530809/1‑68 (MQ=37)
                                   aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:36462984/1‑68 (MQ=37)
                                   aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:42260917/1‑68 (MQ=37)
                                     gATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACATAtt                                 >  1:4469444/1‑67 (MQ=18)
                                       ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGt                              <  1:43638755/68‑1 (MQ=255)
                                       ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGt                              <  1:29792438/68‑1 (MQ=255)
                                       ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGt                              <  1:38165117/68‑1 (MQ=255)
                                       ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGt                              <  1:35100734/68‑1 (MQ=255)
                                            aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTAATTGGtatttc                         >  1:32043687/1‑66 (MQ=11)
                                            aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGa                         >  1:47133439/1‑68 (MQ=255)
                                                cGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAAttt                     <  1:13585662/68‑1 (MQ=255)
                                                   aaTGCCGGAATCCCAATCCGGAGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGc                   <  1:18811705/67‑1 (MQ=18)
                                                   aaTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGc                   <  1:22561455/67‑1 (MQ=255)
                                                   aaTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGc                   <  1:19997410/67‑1 (MQ=255)
                                                    aTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAATTTGGTTTTGAATTTGcc                  >  1:5142386/1‑67 (MQ=255)
                                                    aTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGATTTGAATTTGCCg                 >  1:2118644/1‑68 (MQ=255)
                                                     tGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa                <  1:4321741/68‑1 (MQ=255)
                                                      gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa               <  1:14463705/68‑1 (MQ=255)
                                                      gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa               <  1:43128332/68‑1 (MQ=255)
                                                      gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa               <  1:19621601/68‑1 (MQ=255)
                                                      gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa               <  1:12924270/68‑1 (MQ=255)
                                                      gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa                >  1:2328170/1‑67 (MQ=255)
                                                      gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAATGGATGTGAATTTGccta                >  1:28620505/1‑65 (MQ=17)
                                                       ccGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTGTTGAATTTTACGAAc              >  1:30961362/1‑68 (MQ=18)
                                                        cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAAc              <  1:29282332/67‑1 (MQ=255)
                                                        cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACa             <  1:12661808/68‑1 (MQ=255)
                                                         ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat            <  1:41468458/68‑1 (MQ=255)
                                                         ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat            <  1:13341520/68‑1 (MQ=255)
                                                         ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat            <  1:3480985/68‑1 (MQ=255)
                                                           aaTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACata           <  1:11425361/67‑1 (MQ=255)
                                                           aaTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACata           <  1:4107477/67‑1 (MQ=255)
                                                           aaTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACata           <  1:26386666/67‑1 (MQ=255)
                                                            aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat          <  1:48300694/67‑1 (MQ=255)
                                                            aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat          <  1:8412465/67‑1 (MQ=255)
                                                            aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat          <  1:3405383/67‑1 (MQ=255)
                                                            aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat          <  1:44111034/67‑1 (MQ=255)
                                                            aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat          <  1:46256392/67‑1 (MQ=255)
                                                             tCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTc        <  1:25665200/68‑1 (MQ=255)
                                                              cAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGATCATATTCg       >  1:37805695/1‑68 (MQ=255)
                                                                  aTCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGa      <  1:39898489/65‑1 (MQ=255)
                                                                  aTCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATa    <  1:15890489/67‑1 (MQ=255)
                                                                   tCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACa  <  1:27379761/68‑1 (MQ=255)
                                                                   |                                                                   
AAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACA  >  CP009273/691976‑692110

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: