Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,268,331 G→A 8.6% L123L (TTG→TTA chaC → cation transport regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,268,3310GA 8.6% 1223.3 / 79.3 394L123L (TTG→TTAchaCcation transport regulator
Reads supporting (aligned to +/- strand):  ref base G (164/195);  new base A (20/14);  total (185/209)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 431 total aligned reads displayed.

CGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGAC  >  CP009273/1268265‑1268398
                                                                  |                                                                   
   aCGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTgg                                                                  <  1:12045874/67‑1 (MQ=255)
    cGCACCACAGGCGTCGCCTATCGACTGGCAGAAGAGACGCTGGAGCAGGAACTAACCCGGTTGTGGaa                                                                >  1:13624148/1‑68 (MQ=255)
     gCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAg                                                               <  1:11301008/68‑1 (MQ=255)
     gCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAg                                                               <  1:32166085/68‑1 (MQ=255)
     gCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTATGGaa                                                                >  1:3668043/1‑67 (MQ=255)
     gCACCACAGGCGGCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAg                                                               <  1:23673505/68‑1 (MQ=255)
      caccacAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGTACTAACCCTGTTGTGGAAGc                                                              >  1:31108584/1‑68 (MQ=255)
         cacaGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCgag                                                           <  1:25684815/68‑1 (MQ=255)
         cacaGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCgag                                                           <  1:2687251/68‑1 (MQ=255)
          acaGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCgaga                                                          >  1:17708276/1‑68 (MQ=255)
            aGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTATGGAAGCgagaga                                                        >  1:13369512/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:26195746/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:32882434/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:4875243/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:33539045/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:6033018/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:9081463/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:9549122/1‑68 (MQ=255)
             ggCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:34055801/1‑68 (MQ=255)
             ggCGTCGCCTATCGACAGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgat                                                       >  1:25817713/1‑68 (MQ=255)
               cggcGCCTATCGACTGCCAGACGACACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCgagaga                                                        <  1:25606751/62‑1 (MQ=255)
               cggcGCCTATCGACTGCCAGAAGAGCCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCgagaga                                                        <  1:22904229/62‑1 (MQ=255)
               cGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTACCCCTGTTGTGGAAGCGAGAgatga                                                     >  1:31950786/1‑68 (MQ=255)
               cGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatga                                                     >  1:5558215/1‑68 (MQ=255)
               cGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatga                                                     >  1:21997844/1‑68 (MQ=255)
               cGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatga                                                     >  1:4531414/1‑68 (MQ=255)
               cGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatga                                                     >  1:12293712/1‑68 (MQ=255)
               cGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatga                                                     >  1:11594577/1‑68 (MQ=255)
                gTCGCCTATCGACTGCCAGAAGAGCCGCTGGAGCAGGAACTAACCCTGGTGTGGAAGCGAGAgatgat                                                    <  1:18167582/68‑1 (MQ=255)
                gTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatgat                                                    <  1:10004097/68‑1 (MQ=255)
                gTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatgat                                                    <  1:13960104/68‑1 (MQ=255)
                gTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatgat                                                    <  1:15875856/68‑1 (MQ=255)
                       aTCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGc                                             <  1:17389695/68‑1 (MQ=255)
                       aTCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGc                                             <  1:30121262/68‑1 (MQ=255)
                         cGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCg                                               >  1:3434920/1‑64 (MQ=255)
                           aCTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGtt                                         >  1:16875344/1‑68 (MQ=255)
                           aCTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGtt                                         >  1:11077960/1‑68 (MQ=255)
                           aCTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGCGGAAGCGAGAGATGATTACCGGCTGtt                                         >  1:7426357/1‑68 (MQ=255)
                              gCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTATGGAAGCGAGAGATGATACCGGCTGTTAtat                                     >  1:4366937/1‑66 (MQ=255)
                                cAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCt                                     <  1:20277504/67‑1 (MQ=255)
                                cAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCt                                     <  1:29583200/67‑1 (MQ=255)
                                cAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCt                                     <  1:30327176/67‑1 (MQ=255)
                                cAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCt                                     <  1:12764177/67‑1 (MQ=255)
                                  gaagaGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGcc                                  >  1:23759235/1‑68 (MQ=255)
                                    agagACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                <  1:17969979/68‑1 (MQ=255)
                                    agagACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                <  1:20844930/68‑1 (MQ=255)
                                    agagACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                <  1:33182856/68‑1 (MQ=255)
                                     gagaCGCTGGAGCAGGAACTACCCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                >  1:17044555/1‑67 (MQ=255)
                                     gagaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                >  1:22684506/1‑67 (MQ=255)
                                     gagaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                >  1:25095139/1‑67 (MQ=255)
                                     gagaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                >  1:27984803/1‑67 (MQ=255)
                                     gagaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCAGCaaa                                >  1:29114495/1‑67 (MQ=255)
                                     gagaCGCTGGAGCAGGAACTAACCCTGTTATGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCaa                                >  1:18670156/1‑67 (MQ=255)
                                       gaggCTGGACCAGGAACTAACCCTGTTGTGGATGCGAGAGATGATTACCGGCTGTTATCTGCCAAc                               <  1:5138717/63‑1 (MQ=255)
                                       gaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAAcc                              <  1:4034465/67‑1 (MQ=255)
                                       gaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAAcc                              <  1:9724032/67‑1 (MQ=255)
                                       gaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAAcc                              <  1:1188763/67‑1 (MQ=255)
                                       gaCGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAAc                               <  1:23324769/66‑1 (MQ=255)
                                         cGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTgg                           >  1:2178575/1‑68 (MQ=255)
                                         cGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTgg                           >  1:3059777/1‑68 (MQ=255)
                                         cGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCGGTTATCTGCCAACCTgg                           <  1:31277254/68‑1 (MQ=255)
                                          gctggagctggaACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGgt                          <  1:24631128/68‑1 (MQ=255)
                                          gCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGgt                          <  1:28908938/68‑1 (MQ=255)
                                          gCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGgt                          >  1:12531246/1‑68 (MQ=255)
                                          gCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGgt                          <  1:15345195/68‑1 (MQ=255)
                                                gCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAg                     <  1:7710524/67‑1 (MQ=255)
                                                gCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAg                     <  1:15162197/67‑1 (MQ=255)
                                                gCAGGAACTAACCCTGTTGGGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAg                     <  1:24411910/67‑1 (MQ=255)
                                                  aGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGc                    <  1:18908009/66‑1 (MQ=255)
                                                    gAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATGTGCCAACCTGGTGACAGcatga                >  1:14585593/1‑68 (MQ=255)
                                                     aaCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttga                <  1:6778028/67‑1 (MQ=255)
                                                     aaCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttga                >  1:21155986/1‑67 (MQ=255)
                                                     aaCTAACCCTGTTATGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttga                <  1:30097783/67‑1 (MQ=255)
                                                      aCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgat               <  1:30931248/67‑1 (MQ=255)
                                                      aCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgat               <  1:27176474/67‑1 (MQ=255)
                                                      aCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgat               <  1:32545882/67‑1 (MQ=255)
                                                      aCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgat               <  1:20070340/67‑1 (MQ=255)
                                                      aCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttga                >  1:10287252/1‑66 (MQ=255)
                                                      aCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttga                >  1:8839145/1‑66 (MQ=255)
                                                      aCTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttga                >  1:28437440/1‑66 (MQ=255)
                                                      aCTAACCCTGTTATGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgat               <  1:14146754/67‑1 (MQ=255)
                                                       cTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatct             >  1:15441635/1‑68 (MQ=255)
                                                       cTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatct             <  1:2229785/68‑1 (MQ=255)
                                                          aCCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttga          <  1:15988839/68‑1 (MQ=255)
                                                          aCCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttga          <  1:19008380/68‑1 (MQ=255)
                                                          aCCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttga          >  1:198095/1‑68 (MQ=255)
                                                          aCCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttga          >  1:25353712/1‑68 (MQ=255)
                                                          aCCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttga          >  1:6512891/1‑68 (MQ=255)
                                                           cccTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttgat         >  1:17100805/1‑68 (MQ=255)
                                                           cccTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttgat         >  1:6354028/1‑68 (MQ=255)
                                                           cccTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttgat         >  1:8694394/1‑68 (MQ=255)
                                                           cccAGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatcttgat         >  1:7336766/1‑68 (MQ=255)
                                                            ccTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTtgatg        >  1:13165077/1‑68 (MQ=255)
                                                            ccTGTTATGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTtgatg        >  1:8359729/1‑68 (MQ=255)
                                                             cTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTtgatga       <  1:26441480/68‑1 (MQ=255)
                                                             cTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTtgatga       <  1:21509566/68‑1 (MQ=255)
                                                               gttgtGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGatgatg     <  1:29827831/68‑1 (MQ=255)
                                                               gttgtGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGatgatg     <  1:16193970/68‑1 (MQ=255)
                                                                ttgtGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATgg    >  1:14752357/1‑68 (MQ=255)
                                                                 tgtgGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGa   <  1:20704508/68‑1 (MQ=255)
                                                                  |                                                                   
CGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGAC  >  CP009273/1268265‑1268398

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: