Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,689,128 A→C 27.6% pseudogene (218/753 nt) bcsQ ← involved in cellulose production, minD superfamily (pseudogene)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,689,1280AC27.6% 1212.4 / inf 615pseudogene (218/753 nt)bcsQinvolved in cellulose production, minD superfamily (pseudogene)
Reads supporting (aligned to +/- strand):  ref base A (105/340);  new base C (26/144);  total (131/484)
Fisher's exact test for biased strand distribution p-value = 2.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 627 total aligned reads displayed.

GGATTTTCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTAAAATC  >  CP009273/3689062‑3689195
                                                                  |                                                                   
ggATTTTCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCa                                                                     >  1:35080493/1‑67 (MQ=255)
  aTTTTCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAAc                                                                   <  1:2696734/67‑1 (MQ=255)
  aTTTTCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAAc                                                                   <  1:2779154/67‑1 (MQ=255)
  aTTTTCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAAc                                                                   <  1:15118419/67‑1 (MQ=255)
  aTTTTCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAAc                                                                   <  1:24277073/67‑1 (MQ=255)
     ttCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCa                                                                >  1:25473686/1‑67 (MQ=255)
      tCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGc                                                              <  1:29308113/68‑1 (MQ=255)
      tCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGc                                                              <  1:7111971/68‑1 (MQ=255)
       cTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCg                                                             >  1:37513677/1‑68 (MQ=255)
            tcttcAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCAc                                                         >  1:46861103/1‑67 (MQ=255)
            tcttcAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCAc                                                         >  1:10095428/1‑67 (MQ=255)
            tcttcAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCAc                                                         >  1:124133/1‑67 (MQ=255)
              ttcAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGcc                                                      <  1:6833966/68‑1 (MQ=255)
               tcAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCa                                                     >  1:37145907/1‑68 (MQ=255)
                cAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAg                                                    <  1:1540607/68‑1 (MQ=255)
                cAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAg                                                    <  1:34707182/68‑1 (MQ=255)
                cAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAg                                                    <  1:39133664/68‑1 (MQ=255)
                cAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAg                                                    <  1:42964093/68‑1 (MQ=255)
                cAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAg                                                    <  1:11902976/68‑1 (MQ=255)
                cAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAg                                                    <  1:26516868/68‑1 (MQ=255)
                 aaTGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGt                                                   <  1:5314209/68‑1 (MQ=255)
                 aaTGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGt                                                   <  1:10877918/68‑1 (MQ=255)
                 aaTGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGt                                                   <  1:12974465/68‑1 (MQ=255)
                 aaTGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGt                                                   <  1:49517706/68‑1 (MQ=255)
                 aaTGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGt                                                   <  1:11146736/68‑1 (MQ=255)
                 aaTGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGGAGCGCAACCCAGCGTCACGCCAgg                                                   <  1:43259543/68‑2 (MQ=255)
                 aaTGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGGGTAGCGCAACCCAGCGTCACGCCAGt                                                   <  1:25758401/68‑1 (MQ=255)
                    ggATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGt                                                   >  1:35810025/1‑65 (MQ=255)
                      aTAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTg                                               >  1:5995949/1‑67 (MQ=255)
                      aTAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTg                                               <  1:367005/67‑1 (MQ=255)
                      aTAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTg                                               <  1:17068397/67‑1 (MQ=255)
                      aTAACTGACCAAAAGGCAGCAAAACGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTg                                               >  1:11602742/1‑67 (MQ=255)
                        aaCTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGcc                                            <  1:46313953/68‑1 (MQ=255)
                         aCTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCa                                           >  1:30903700/1‑68 (MQ=255)
                         aCTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCa                                           >  1:29804100/1‑68 (MQ=255)
                         aCTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCa                                           >  1:44571135/1‑68 (MQ=255)
                          cTGACGAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCAt                                          <  1:14052371/68‑1 (MQ=255)
                          cTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCAt                                          <  1:19963485/68‑1 (MQ=255)
                          cTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCAt                                          <  1:23926124/68‑1 (MQ=255)
                          cTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCAt                                          <  1:18306777/68‑1 (MQ=255)
                           tGACCAAAAGGCAGCAAATCGAGCTGCGAGTTGTAGCGCAACCCAGCGACACGCCAGTCCTGGCCAt                                          >  1:1769851/1‑67 (MQ=255)
                           tGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCAt                                          >  1:24736318/1‑67 (MQ=255)
                              ccAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCa                                       <  1:8399639/67‑1 (MQ=255)
                               cAAAAGGCAGCAAATCGCGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCATTCCTGGCCATCCAGc                                     <  1:23175185/68‑1 (MQ=255)
                               cAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGc                                     <  1:31581387/68‑1 (MQ=255)
                                   aGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATc                                  >  1:17427991/1‑67 (MQ=255)
                                   aGGCACCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCg                                 <  1:18816973/68‑1 (MQ=255)
                                      ctgcaAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGACAGTCCTGGCCATCCAGCATCGc                                <  1:26330284/64‑1 (MQ=255)
                                      cagcaAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGct                               <  1:12606268/67‑1 (MQ=255)
                                      cagcaAATCGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGc                                <  1:50294139/66‑1 (MQ=255)
                                      cagcaAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGct                               <  1:2217097/67‑1 (MQ=255)
                                      cagcaAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGct                               <  1:43695445/67‑1 (MQ=255)
                                      cagcaAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGct                               <  1:28287499/67‑1 (MQ=255)
                                      cagcaAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGc                                <  1:36497076/66‑1 (MQ=255)
                                      cagcaAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGc                                <  1:19376222/66‑1 (MQ=255)
                                      cagcaAATCGAGCGGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGct                               <  1:21539715/67‑1 (MQ=255)
                                        gcaAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTg                            >  1:18961901/1‑68 (MQ=255)
                                          aaaTCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCCTCGCTCTGGc                          <  1:50629286/68‑1 (MQ=255)
                                          aaaTCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTgg                           <  1:28686334/67‑1 (MQ=255)
                                          aaaTCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTgg                           <  1:4516554/67‑1 (MQ=255)
                                          aaaTCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGc                          <  1:39979771/68‑1 (MQ=255)
                                          aaaTCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGc                          <  1:13434686/68‑1 (MQ=255)
                                              cGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAg                      <  1:33831454/68‑1 (MQ=255)
                                              cGAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAg                      <  1:28108094/68‑1 (MQ=255)
                                              cGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAg                      <  1:45009486/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:38764444/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:48315428/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:38537898/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:45511397/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:1457120/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:9306254/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:42013759/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     <  1:36248126/68‑1 (MQ=255)
                                               gAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATAGCTCTGGCCCAGc                     <  1:8733208/68‑1 (MQ=255)
                                                aGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGc                     >  1:10349739/1‑67 (MQ=255)
                                                 gCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:32636970/67‑1 (MQ=255)
                                                 gCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:20942461/67‑1 (MQ=255)
                                                 gCTGCGAGGTGTAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:21857003/67‑1 (MQ=255)
                                                 gCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:16157279/67‑1 (MQ=255)
                                                 gCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:32055267/67‑1 (MQ=255)
                                                 gCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:33146156/67‑1 (MQ=255)
                                                 gCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:4725459/67‑1 (MQ=255)
                                                 gCTCCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGcc                    <  1:47951335/67‑1 (MQ=255)
                                                       aGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACgg             >  1:34281564/1‑68 (MQ=255)
                                                          tgtAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTg           <  1:44098366/67‑1 (MQ=255)
                                                           gtAGCGCCACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTgg          <  1:41372098/67‑1 (MQ=255)
                                                           gtAGCGCAACCCAGCGTCCCGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTgg          <  1:39470286/67‑1 (MQ=255)
                                                           gtAGCGCAACCCAGCGTCACGCCAGTCCTGGCGATCCAGCATCGCTCTGGCCCAGCCCTGACGGTggg         <  1:45848056/68‑1 (MQ=255)
                                                           gtAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTggg         >  1:48605594/1‑68 (MQ=255)
                                                           gtAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTggg         <  1:20702143/68‑1 (MQ=255)
                                                           gtAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTgg          <  1:7543937/67‑1 (MQ=255)
                                                           gtAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTgg          <  1:16689441/67‑1 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    <  1:49229851/68‑1 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    >  1:37714218/1‑68 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    >  1:32490431/1‑68 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    >  1:31040632/1‑68 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    >  1:3044491/1‑68 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    >  1:40684062/1‑68 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    <  1:40915769/68‑1 (MQ=255)
                                                                gcAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTaaaa    >  1:9768632/1‑68 (MQ=255)
                                                                  |                                                                   
GGATTTTCTTGTTCTTCAATGGATAACTGACCAAAAGGCAGCAAATCGAGCTGCGAGGTGTAGCGCAACCCAGCGTCACGCCAGTCCTGGCCATCCAGCATCGCTCTGGCCCAGCCCTGACGGTGGGTAAAATC  >  CP009273/3689062‑3689195

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: