Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,187,952 A→G 6.2% intergenic (+68/+130) yqiK → / ← sibD PHB family membrane protein, function unknown/sRNA antisense regulator of toxic IbsD protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,187,9520AG 6.2% 1233.6 / 41.7 393intergenic (+68/+130)yqiK/sibDPHB family membrane protein, function unknown/sRNA antisense regulator of toxic IbsD protein
Reads supporting (aligned to +/- strand):  ref base A (217/146);  new base G (5/19);  total (228/165)
Fisher's exact test for biased strand distribution p-value = 2.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 451 total aligned reads displayed.

ACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGA  >  CP009273/3187885‑3188019
                                                                   |                                                                   
 cccTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAg                                                                    >  1:15368315/1‑68 (MQ=255)
  ccTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                                                                   <  1:27720902/68‑1 (MQ=255)
  ccTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                                                                   <  1:517319/68‑1 (MQ=255)
   ctctTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTg                                                                  <  1:11651247/68‑1 (MQ=255)
   ctctTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                                                                   >  1:2963315/1‑67 (MQ=255)
   ctctTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                                                                   >  1:21970728/1‑67 (MQ=255)
   ctctTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGt                                                                   >  1:18552992/1‑67 (MQ=255)
     ctTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCa                                                                <  1:18002034/68‑1 (MQ=255)
       tCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCaaa                                                              >  1:25007184/1‑68 (MQ=255)
        ccATGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAAt                                                             <  1:26010385/68‑1 (MQ=255)
        ccAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt                                                             >  1:8524564/1‑68 (MQ=255)
         cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg                                                            <  1:3611090/68‑1 (MQ=255)
         cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg                                                            <  1:11485423/68‑1 (MQ=255)
         cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg                                                            <  1:21470393/68‑1 (MQ=255)
         cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATg                                                            <  1:7513059/68‑1 (MQ=255)
          atgAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATgc                                                           <  1:9629304/66‑2 (MQ=255)
          acgAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATgc                                                           <  1:10973070/66‑2 (MQ=255)
          aGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATgc                                                           <  1:1065594/68‑2 (MQ=255)
            gAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATc                                                         >  1:11823488/1‑68 (MQ=255)
            gAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATc                                                         >  1:20328931/1‑68 (MQ=255)
            gAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATc                                                         >  1:19676723/1‑68 (MQ=255)
             aataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAATGTGCAAATGATCa                                                        >  1:14373922/1‑68 (MQ=255)
             aataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATc                                                         <  1:18382304/67‑1 (MQ=255)
             aataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCa                                                        >  1:8359729/1‑68 (MQ=255)
              ataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg                                                       >  1:27987042/1‑68 (MQ=255)
              ataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg                                                       >  1:2729443/1‑68 (MQ=255)
              ataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAg                                                       >  1:30159391/1‑68 (MQ=255)
                  tCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttg                                                   >  1:9083665/1‑68 (MQ=255)
                  tCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtag                                                   >  1:26697044/1‑66 (MQ=255)
                   cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGATCAggtggt                                                  >  1:29047491/1‑68 (MQ=255)
                   cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAggtggt                                                  >  1:25682004/1‑68 (MQ=255)
                   cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgt                                                  >  1:30671536/1‑68 (MQ=255)
                   cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgt                                                  >  1:17470608/1‑68 (MQ=255)
                   cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgt                                                  >  1:24723756/1‑68 (MQ=255)
                   cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgt                                                  >  1:27043626/1‑68 (MQ=255)
                   cccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgatgt                                                  >  1:426856/1‑68 (MQ=255)
                    ccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgtt                                                 <  1:31616603/68‑1 (MQ=255)
                    ccTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgtt                                                 <  1:17006366/68‑1 (MQ=255)
                      tGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAATGTGCAAATGATCAGGATGTTgc                                               >  1:31482860/1‑68 (MQ=255)
                      tGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgc                                               <  1:6199382/68‑1 (MQ=255)
                         cccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcagc                                            >  1:29463875/1‑68 (MQ=255)
                          ccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcagca                                           >  1:9220840/1‑68 (MQ=255)
                          ccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgcagc                                            <  1:21706625/67‑1 (MQ=255)
                           cAAACAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAt                                          <  1:28273893/68‑1 (MQ=255)
                           cAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAt                                          <  1:13212744/68‑1 (MQ=255)
                           cAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAt                                          <  1:4717762/68‑1 (MQ=255)
                              aaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAtttt                                       <  1:22419710/68‑1 (MQ=255)
                              aaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAtttt                                       <  1:5441294/68‑1 (MQ=255)
                              aaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAtttt                                       <  1:7821268/68‑1 (MQ=255)
                               aaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAtttt                                       >  1:5676168/1‑67 (MQ=255)
                                aGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGc                                     >  1:12547066/1‑68 (MQ=255)
                                aGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGc                                     >  1:31738941/1‑68 (MQ=255)
                                 gCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCg                                    >  1:31959729/1‑68 (MQ=255)
                                 gCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCg                                    >  1:1885854/1‑68 (MQ=255)
                                  ccGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGt                                   >  1:22039107/1‑68 (MQ=255)
                                        aTTTTTTATTTCTGCCTTCCGCTATAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACgg                              <  1:17300800/67‑1 (MQ=255)
                                        aTTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCCGGTTGTTGCAGCATTTTGCGTGACgg                              <  1:31893902/67‑1 (MQ=255)
                                         ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                             <  1:4877047/67‑1 (MQ=255)
                                         ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                             <  1:128221/67‑1 (MQ=255)
                                         ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                             <  1:22309013/67‑1 (MQ=255)
                                         ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                             <  1:31245933/67‑1 (MQ=255)
                                         ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                             <  1:3401928/67‑1 (MQ=255)
                                         ttttttATATCTCCGTTCCGCTAAAAAGTGCAAATGATCACGTTGTTGCAGCATTTTGCGTGACGGc                             <  1:29773726/67‑1 (MQ=18)
                                         ttttttATATCTCCGTTCCGCTAAAAAGGGCAAATGATCACGTTGGTGCAGCATTTTGCGTGACGGc                             <  1:26460876/67‑1 (MQ=18)
                                             ttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCggg                        <  1:6821867/68‑1 (MQ=255)
                                             ttATATCTGCGTTCCGCAAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCggg                        <  1:19899857/68‑1 (MQ=255)
                                                tatCTGCGTTCCGCTAAAAAGTGCAAATGCTCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGaa                     >  1:15140703/1‑68 (MQ=255)
                                                tatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGaa                     >  1:18987872/1‑68 (MQ=255)
                                                tatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGaa                     >  1:15848124/1‑68 (MQ=255)
                                                 atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                    <  1:6633081/68‑1 (MQ=255)
                                                 atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                    <  1:19113858/68‑1 (MQ=255)
                                                 atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                    >  1:23672329/1‑68 (MQ=255)
                                                  tCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGTCGGCTCGGGGAAg                    >  1:29274406/1‑67 (MQ=255)
                                                  tCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                    >  1:21019827/1‑67 (MQ=255)
                                                  tCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                    >  1:13474763/1‑67 (MQ=255)
                                                  tCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                    >  1:23339695/1‑67 (MQ=255)
                                                  tCTGCGTCCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAg                    >  1:987400/1‑67 (MQ=255)
                                                   cTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGc                  <  1:30336287/68‑1 (MQ=255)
                                                       gTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaa               >  1:29758457/1‑67 (MQ=255)
                                                       gTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaa               <  1:23042805/67‑1 (MQ=255)
                                                         tCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaaagt            >  1:19191534/1‑66 (MQ=255)
                                                         tCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaa               >  1:13808799/1‑65 (MQ=255)
                                                         tCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaa               >  1:4926070/1‑65 (MQ=255)
                                                         tCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaaa               >  1:24361285/1‑65 (MQ=255)
                                                           cGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGc          <  1:6369080/68‑1 (MQ=255)
                                                           cGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGc          <  1:17838256/68‑1 (MQ=255)
                                                           cGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGc          <  1:28603495/68‑1 (MQ=255)
                                                            gCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCa         >  1:19304034/1‑68 (MQ=255)
                                                                aaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCag     <  1:14717422/68‑1 (MQ=255)
                                                                aaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGGGACGGCTCGGGGAAGGCAAAACTGCATCag     <  1:30500321/68‑1 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGCCGGCTCGGGGAAGGCAAAACTGAATCaga    >  1:2078972/1‑68 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAAGTGCATCaga    >  1:10056310/1‑68 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga    >  1:22769178/1‑68 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga    >  1:7300017/1‑68 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga    >  1:2403184/1‑68 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga    >  1:25828463/1‑68 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga    <  1:16283504/68‑1 (MQ=255)
                                                                 aaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga    >  1:16056900/1‑68 (MQ=255)
                                                                  aaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCagag   <  1:5920697/68‑1 (MQ=255)
                                                                   aGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAgaga  >  1:10345155/1‑68 (MQ=255)
                                                                   |                                                                   
ACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGA  >  CP009273/3187885‑3188019

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: