Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,029,468 A→C 17.8% T292P (ACC→CCC)  mcpU_3 → Methyl‑accepting chemotaxis protein McpU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,029,4680AC17.8% 113.8 / 3.2 73T292P (ACC→CCC) mcpU_3Methyl‑accepting chemotaxis protein McpU
Reads supporting (aligned to +/- strand):  ref base A (23/34);  new base C (12/4);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 2.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGTG  >  ACB122/4029338‑4029590
                                                                                                                                  |                                                                                                                           
gAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCCAGCTGGCCCCGCGCGTCGCCCAGGTAGCCCCGAccca                                                                                                                         >  7:347556/1‑133 (MQ=255)
  gCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAg                                                                                                                       <  5:328812/135‑1 (MQ=255)
     aggaggCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCCAGCTGGCCCAGCGCTTCGACCAGGGAGCCCCGACCAAAG‑gcg                                                                                                                   >  8:357387/1‑135 (MQ=255)
         ggCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCCGGGAGACCCGACCAAAG‑GCGCCGa                                                                                                               >  3:447950/1‑135 (MQ=255)
         ggCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCGGGCTACAGGGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGa                                                                                                               <  2:431132/135‑1 (MQ=255)
          gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCCCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCCCCACGTGGCCCAGCCCATGGACCAGGGAGCCCCGCCCAAAG‑GCGGCGAg                                                                                                              >  1:103828/1‑135 (MQ=255)
          gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGCTTCCCGCAAGCGGGCCCCGCGCGTCGACCAGGTAGACCCCACCAAAG‑GCCCCgc                                                                                                               >  3:299767/1‑133 (MQ=255)
          gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCAAGCGGGCCCAGCGCGTCGCCCAGGTAGCCCCGCCCACAG‑GCCCCgcg                                                                                                              >  3:294037/1‑133 (MQ=255)
          gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGACCAGGTAGACCCGCCCCAAG‑GCGCCGAg                                                                                                              >  8:411390/1‑135 (MQ=255)
          gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGCCCAGGGAGCCCCGACCAAAG‑GCGCCGAg                                                                                                              >  3:16951/1‑135 (MQ=255)
             aGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGGCGGTCGGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTg                                                                                                           <  2:350317/135‑1 (MQ=255)
              gCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCCTTCTCCCCCCGGTCGGGCTGTTGGCTGGCTACACCGCCGCTGCCAGCAAGCTGCCCCCGCGCTTCGCCCAGGTAGCCCCGACCAAAC‑GCGCCGCGCAgg                                                                                                          >  7:219961/1‑135 (MQ=255)
              gCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCCCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCCGGGAGGCCCGCCCAAAG‑GCGCCGCGCTgg                                                                                                          >  7:17818/1‑135 (MQ=255)
                cGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCAAGCTGGCCCAGCGCGTCGCCCCGGGAGCCCCGACCAAAG‑GCGCCGAGCTGGTg                                                                                                        >  7:135517/1‑135 (MQ=255)
                cGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGCCCAGGTAGACCCGCCCAAAC‑GCGCCGAGCTGGTg                                                                                                        >  6:166937/1‑135 (MQ=255)
                       tctACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCCAGCTGGCCCCGCGCGTCGCCCAGGTAGACCCGCCCAAAG‑GCGCCGCGCTGGTGCGCAAGc                                                                                                 >  5:398831/1‑135 (MQ=255)
                            gAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGCCCAGGTAGCCCCGACCAAAG‑GCGCCGCGCTGGTGCGCAAGCTGGcc                                                                                            >  6:224183/1‑135 (MQ=255)
                            gAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGcc                                                                                            >  8:441966/1‑135 (MQ=255)
                               ggCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGAt                                                                                         <  3:409443/135‑1 (MQ=255)
                               ggCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCGGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGAt                                                                                         <  4:410739/135‑1 (MQ=255)
                                  cGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCGTCGACCAGGGAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATgg                                                                                       >  3:294103/1‑134 (MQ=255)
                                    caccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTTGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCaa                                                                                    >  3:126139/1‑135 (MQ=255)
                                    caccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCaa                                                                                    <  7:441966/135‑1 (MQ=255)
                                     accacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAg                                                                                   >  7:296861/1‑135 (MQ=255)
                                      ccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGCCCAGGTAGCCCCGCCCAAAG‑GCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGc                                                                                  >  8:42613/1‑135 (MQ=255)
                                      ccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAACCTGGCCCCGCGCTTCGCCCAGGTAGCCCCGACCCAAG‑GCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGc                                                                                  >  6:166651/1‑135 (MQ=255)
                                      ccacGGTCGGCATTCTCACCCCCGTCGGTCTGTTGGCTGGCTACAGCGCCGAGGCCAGCAAACTGGGCCAGCCCGTGCTCCCCGTATACCCCCCCCAAG‑GCGGCCCGCATGGGCGCAATCAGGCCGATGGCCAGc                                                                                  >  6:5712/1‑135 (MQ=255)
                                       cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCCAGCTGGCCCAGCGCTTAGCCCCGGTAGACCCGCCACAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGCTGGCCAGCt                                                                                 >  8:30773/1‑135 (MQ=255)
                                       cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCAGGTAGACCCGACCCAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGCTGGcccgct                                                                                 >  8:23554/1‑130 (MQ=255)
                                       cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                                 <  6:368508/135‑1 (MQ=255)
                                       cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                                 <  3:350529/135‑1 (MQ=255)
                                       cacGGTCGGCATTCTCACCCCGGTCGGTCTGTGGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                                 <  2:379464/135‑1 (MQ=255)
                                       cacGGTCGGCATTCTCACCCCGGTCGGTCGGTGGGCGGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                                 <  7:411390/135‑1 (MQ=255)
                                         cGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGa                                                                               <  5:204679/135‑1 (MQ=255)
                                           gTCGGCTTTCTCACCCCGGGCGGTCTGTTGGCGGGCTACTGCGCCGATGCCAGCAAGGTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGacc                                                                             <  7:362800/135‑1 (MQ=255)
                                            tCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCAGGTAGCCCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGacca                                                                            >  8:259691/1‑135 (MQ=255)
                                                    ctcACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTcccc                                                                    >  8:467280/1‑133 (MQ=255)
                                                       aCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTGCGCCCAGGGAGACCCGCCCAAAG‑GCGCCCAGCTGGTGCGCAAGCTGGCCGATGGCCAGCTGACCCCCCTCCACGAc                                                                 >  8:335543/1‑135 (MQ=255)
                                                                ggtcTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAgcgc                                                        >  7:262397/1‑135 (MQ=255)
                                                                              ctACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTg                                            <  8:325322/133‑1 (MQ=255)
                                                                               tACAGCGCCGATGCCCGCAAGCTGGCCCAGCGCTTCGCCCAGGTTGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCCAGCGGACCACCCTCCACGCCCGGCCGCGCGTGAAGGGGCTggcg                                         >  3:139421/1‑135 (MQ=255)
                                                                               tACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTggcg                                         >  8:141102/1‑135 (MQ=255)
                                                                               tACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCGTCGCCCAGGTAGACCCGACCAAAG‑GCGCCGCGCTGGTGCGCAAGCTGGCCGCTGGCCAGCTGACCACCCGCCCCGCCCCGCAGCGCGTGAAGGTGCTggcg                                         >  4:18055/1‑135 (MQ=255)
                                                                                  cgcgcGGATGCCTGCAAGCTGGACCAGGGCTTCGACCAGGTACACCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCGGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTggcggc                                       <  6:190614/133‑1 (MQ=255)
                                                                                     gcCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCAGGGTGACCCGCCCAAAG‑GCGCCGCGCAGGGGCGCCAGCTGGCCGCTGGCAAGCGGACCACCCACCACGCCCGGCAGCGCGTGAACGTGCTGGCGGCATTc                                   >  2:458968/1‑135 (MQ=255)
                                                                                      cggttGCCAGCACGCGGCCCAAGCTCTGCGCCCCGGTAGACCCCACCTACG‑GCTCCGAGAGGGTGCGCAAGCTGCCCGAGGGCAAGCTGCCCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTcc                                  <  4:249965/131‑1 (MQ=255)
                                                                                            ccAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGCCCAAAG‑GCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCCCCACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCaa                            >  4:426813/1‑135 (MQ=255)
                                                                                             cAGCAAGCGGGCCCAGCGCTTCGAACAGGTAGACCCGACCAACT‑GCGCCGAGCTGTTGCGCAAGCTGGCCTATGGAAAGCTGACCACACTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAAt                           <  7:357387/135‑1 (MQ=255)
                                                                                                cacgCTGGCCCTGCGCTTCGACGAGGTCTACCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGATGGTGCTGGCGGCATTCCGGCCAATccc                        <  2:88674/132‑1 (MQ=255)
                                                                                                    cgggCCCAGGGCTACTCCCCGGTACCCCCTCCAAACG‑GCGCCGAGCGGGTGCGCATGGTGGCCGATGCCAAGCTGACCACCCGCCACGACCGGCAGCGCTTGACGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  2:430250/133‑1 (MQ=255)
                                                                                                    cgggCCCAGCGCTTCGACCAGGCAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCATGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  8:245028/133‑1 (MQ=255)
                                                                                                    cgggCCCAGCGCGTCGACCCGGTCGCCCCGCCCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  8:128725/133‑1 (MQ=255)
                                                                                                    cggccccAGCGCTCCCACCTGGACCCCCCCACCACCG‑GCGCCGAGGTGGTGCGCAAGCTGGCCGAGGGAAAGCTCACCCCCCGCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  3:18055/131‑1 (MQ=255)
                                                                                                    cTGGCCCTGCGCTTCTCCCCGGAAGCCCCGCCCAACG‑GCGCCGAGCTGGGGCGCAAGCGGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  5:447150/135‑1 (MQ=255)
                                                                                                    cTGGCCCAGGGTTTCGACCAGCAATACCCGACAAATG‑GCGCCGAGGTGGTGGGAAAGCTGGCCTATGGCAAGGTTACCACCCCCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  5:23663/135‑1 (MQ=255)
                                                                                                    cTGGCCCAGCGCTTCGACCAGGTCGACCCGACCAAAG‑GCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  8:241131/135‑1 (MQ=255)
                                                                                                    cTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  8:262397/135‑1 (MQ=255)
                                                                                                    cTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  7:318029/135‑1 (MQ=255)
                                                                                                    cTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGACCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  7:393027/135‑1 (MQ=255)
                                                                                                    cTGGCCCAGCGCTTCGACCAGGTACACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCTATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  8:14297/135‑1 (MQ=255)
                                                                                                    cTGGCCCAGCGCTTCCACCAGGTAGACCCGACCAAAG‑CCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  2:33671/135‑1 (MQ=255)
                                                                                                    cTGCCCCTGCGCTTCGCCCAGGTAGACCCCACCAAGC‑GCCCCGCGCTGGTGCGAAAGCTGGCCTATTGAAAGCTGACCCCCCACAAAGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  8:402484/135‑1 (MQ=255)
                                                                                                     tGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGCGCGCCGAGCTGGGGCGCACGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                    <  2:167167/135‑1 (MQ=255)
                                                                                                              gcTTCGACCAGGTAGACCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCt          <  8:296861/135‑1 (MQ=255)
                                                                                                                ttCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTgg        >  4:363388/1‑135 (MQ=255)
                                                                                                                 tCGACCCGGTCGACCCGCCACACG‑GGGCGGGGGGGGTGCGCAAGCGGGCCGAGGGCATGCTGCCCACCCTCCCCGACCGGCCGCGCGTGAAGGGGCGGGCGGCATGCCGGCCAATCCCGGACGCCCAGCCCTggg       <  4:299767/135‑1 (MQ=255)
                                                                                                                   gcccAGGTAGACCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTgggg      <  4:294103/132‑1 (MQ=255)
                                                                                                                   gACCAGGTAGACCCGACCAAAG‑GCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTggg       <  8:17818/133‑1 (MQ=255)
                                                                                                                   caCCAGGAACACCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTgggg      <  7:41128/133‑1 (MQ=255)
                                                                                                                    aCCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCg    >  4:459451/1‑135 (MQ=255)
                                                                                                                     ccAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGt   >  6:432589/1‑135 (MQ=255)
                                                                                                                     ccAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGAGGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGt   <  7:295757/135‑1 (MQ=255)
                                                                                                                      cAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGtg  <  7:361485/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                           
GAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCGTG  >  ACB122/4029338‑4029590

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