Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,029,454 A→C 27.3% D287A (GAC→GCC)  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,4540AC27.3% 98.2 / 6.9 74D287A (GAC→GCC) mcpU_3Methyl‑accepting chemotaxis protein McpU
Reads supporting (aligned to +/- strand):  ref base A (18/32);  new base C (19/4);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 3.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCCTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCG  >  ACB122/4029323‑4029588
                                                                                                                                   |                                                                                                                                       
tAGCCTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCa                                                                                                                                      <  3:420441/135‑1 (MQ=255)
 aGCCTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAg                                                                                                                                     >  7:14297/1‑135 (MQ=255)
 aGCCTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAg                                                                                                                                     >  4:114308/1‑135 (MQ=255)
    cTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTa                                                                                                                                  <  3:448001/135‑1 (MQ=255)
        aGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCCCCCCGCTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCACGCTCGCCCAGCGCGTCGACCAGGTATAcc                                                                                                                              >  1:36481/1‑135 (MQ=255)
        aGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGACAGCAAGCTGGCCCAGCGCGTCGCCCAGGTAGccc                                                                                                                              >  7:466690/1‑135 (MQ=255)
          gCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGACCAGGGAGACCCg                                                                                                                            >  3:289641/1‑135 (MQ=255)
               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)
                                                                                                                                  gACCAGGTAGACCCGACCAAAG‑GCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTggg     <  8:17818/133‑1 (MQ=255)
                                                                                                                                  caCCAGGAACACCCGACCAACG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTgggg    <  7:41128/133‑1 (MQ=255)
                                                                                                                                   aCCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCg  >  4:459451/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TAGCCTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAG‑GCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCG  >  ACB122/4029323‑4029588

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: ‑