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

TAGCCTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCG  >  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)
                    aggaggCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCCAGCTGGCCCAGCGCTTCGACCAGGGAGCCCCGACCAAAGgcg                                                                                                                 >  8:357387/1‑135 (MQ=255)
                        ggCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCCGGGAGACCCGACCAAAGGCGCCGa                                                                                                             >  3:447950/1‑135 (MQ=255)
                        ggCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCGGGCTACAGGGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGa                                                                                                             <  2:431132/135‑1 (MQ=255)
                         gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCCCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCCCCACGTGGCCCAGCCCATGGACCAGGGAGCCCCGCCCAAAGGCGGCGAg                                                                                                            >  1:103828/1‑135 (MQ=255)
                         gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGCTTCCCGCAAGCGGGCCCCGCGCGTCGACCAGGTAGACCCCACCAAAGGCCCCgc                                                                                                             >  3:299767/1‑133 (MQ=255)
                         gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCAAGCGGGCCCAGCGCGTCGCCCAGGTAGCCCCGCCCACAGGCCCCgcg                                                                                                            >  3:294037/1‑133 (MQ=255)
                         gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGACCAGGTAGACCCGCCCCAAGGCGCCGAg                                                                                                            >  8:411390/1‑135 (MQ=255)
                         gCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGCCCAGGGAGCCCCGACCAAAGGCGCCGAg                                                                                                            >  3:16951/1‑135 (MQ=255)
                            aGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGGCGGTCGGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTg                                                                                                         <  2:350317/135‑1 (MQ=255)
                             gCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCCTTCTCCCCCCGGTCGGGCTGTTGGCTGGCTACACCGCCGCTGCCAGCAAGCTGCCCCCGCGCTTCGCCCAGGTAGCCCCGACCAAACGCGCCGCGCAgg                                                                                                        >  7:219961/1‑135 (MQ=255)
                             gCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCCCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCCGGGAGGCCCGCCCAAAGGCGCCGCGCTgg                                                                                                        >  7:17818/1‑135 (MQ=255)
                               cGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCAAGCTGGCCCAGCGCGTCGCCCCGGGAGCCCCGACCAAAGGCGCCGAGCTGGTg                                                                                                      >  7:135517/1‑135 (MQ=255)
                               cGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCTTCGCCCAGGTAGACCCGCCCAAACGCGCCGAGCTGGTg                                                                                                      >  6:166937/1‑135 (MQ=255)
                                      tctACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCCAGCTGGCCCCGCGCGTCGCCCAGGTAGACCCGCCCAAAGGCGCCGCGCTGGTGCGCAAGc                                                                                               >  5:398831/1‑135 (MQ=255)
                                           gAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGCCCAGGTAGCCCCGACCAAAGGCGCCGCGCTGGTGCGCAAGCTGGcc                                                                                          >  6:224183/1‑135 (MQ=255)
                                           gAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGcc                                                                                          >  8:441966/1‑135 (MQ=255)
                                              ggCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGAt                                                                                       <  3:409443/135‑1 (MQ=255)
                                              ggCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCGGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGAt                                                                                       <  4:410739/135‑1 (MQ=255)
                                                 cGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCGTCGACCAGGGAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATgg                                                                                     >  3:294103/1‑134 (MQ=255)
                                                   caccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTTGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCaa                                                                                  >  3:126139/1‑135 (MQ=255)
                                                   caccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCaa                                                                                  <  7:441966/135‑1 (MQ=255)
                                                    accacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAg                                                                                 >  7:296861/1‑135 (MQ=255)
                                                     ccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGCCCAGGTAGCCCCGCCCAAAGGCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGc                                                                                >  8:42613/1‑135 (MQ=255)
                                                     ccacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAACCTGGCCCCGCGCTTCGCCCAGGTAGCCCCGACCCAAGGCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGc                                                                                >  6:166651/1‑135 (MQ=255)
                                                     ccacGGTCGGCATTCTCACCCCCGTCGGTCTGTTGGCTGGCTACAGCGCCGAGGCCAGCAAACTGGGCCAGCCCGTGCTCCCCGTATACCCCCCCCAAGGCGGCCCGCATGGGCGCAATCAGGCCGATGGCCAGc                                                                                >  6:5712/1‑135 (MQ=255)
                                                      cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCCGCCAGCTGGCCCAGCGCTTAGCCCCGGTAGACCCGCCACAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGCTGGCCAGCt                                                                               >  8:30773/1‑135 (MQ=255)
                                                      cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCAGGTAGACCCGACCCAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGCTGGcccgct                                                                               >  8:23554/1‑130 (MQ=255)
                                                      cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                               <  6:368508/135‑1 (MQ=255)
                                                      cacGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                               <  3:350529/135‑1 (MQ=255)
                                                      cacGGTCGGCATTCTCACCCCGGTCGGTCTGTGGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                               <  2:379464/135‑1 (MQ=255)
                                                      cacGGTCGGCATTCTCACCCCGGTCGGTCGGTGGGCGGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCt                                                                               <  7:411390/135‑1 (MQ=255)
                                                        cGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGa                                                                             <  5:204679/135‑1 (MQ=255)
                                                          gTCGGCTTTCTCACCCCGGGCGGTCTGTTGGCGGGCTACTGCGCCGATGCCAGCAAGGTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGacc                                                                           <  7:362800/135‑1 (MQ=255)
                                                           tCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCAGGTAGCCCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGacca                                                                          >  8:259691/1‑135 (MQ=255)
                                                                   ctcACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTcccc                                                                  >  8:467280/1‑133 (MQ=255)
                                                                      aCCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTGCGCCCAGGGAGACCCGCCCAAAGGCGCCCAGCTGGTGCGCAAGCTGGCCGATGGCCAGCTGACCCCCCTCCACGAc                                                               >  8:335543/1‑135 (MQ=255)
                                                                               ggtcTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAgcgc                                                      >  7:262397/1‑135 (MQ=255)
                                                                                             ctACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTg                                          <  8:325322/133‑1 (MQ=255)
                                                                                              tACAGCGCCGATGCCCGCAAGCTGGCCCAGCGCTTCGCCCAGGTTGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCCAGCGGACCACCCTCCACGCCCGGCCGCGCGTGAAGGGGCTggcg                                       >  3:139421/1‑135 (MQ=255)
                                                                                              tACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTggcg                                       >  8:141102/1‑135 (MQ=255)
                                                                                              tACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCGTCGCCCAGGTAGACCCGACCAAAGGCGCCGCGCTGGTGCGCAAGCTGGCCGCTGGCCAGCTGACCACCCGCCCCGCCCCGCAGCGCGTGAAGGTGCTggcg                                       >  4:18055/1‑135 (MQ=255)
                                                                                                 cgcgcGGATGCCTGCAAGCTGGACCAGGGCTTCGACCAGGTACACCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCGGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTggcggc                                     <  6:190614/133‑1 (MQ=255)
                                                                                                    gcCGATGCCAGCAAGCTGGCCCCGCGCGTCGCCCAGGGTGACCCGCCCAAAGGCGCCGCGCAGGGGCGCCAGCTGGCCGCTGGCAAGCGGACCACCCACCACGCCCGGCAGCGCGTGAACGTGCTGGCGGCATTc                                 >  2:458968/1‑135 (MQ=255)
                                                                                                     cggttGCCAGCACGCGGCCCAAGCTCTGCGCCCCGGTAGACCCCACCTACGGCTCCGAGAGGGTGCGCAAGCTGCCCGAGGGCAAGCTGCCCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTcc                                <  4:249965/131‑1 (MQ=255)
                                                                                                           ccAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGCCCAAAGGCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCCCCACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCaa                          >  4:426813/1‑135 (MQ=255)
                                                                                                            cAGCAAGCGGGCCCAGCGCTTCGAACAGGTAGACCCGACCAACTGCGCCGAGCTGTTGCGCAAGCTGGCCTATGGAAAGCTGACCACACTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAAt                         <  7:357387/135‑1 (MQ=255)
                                                                                                               cacgCTGGCCCTGCGCTTCGACGAGGTCTACCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGATGGTGCTGGCGGCATTCCGGCCAATccc                      <  2:88674/132‑1 (MQ=255)
                                                                                                                   cgggCCCAGGGCTACTCCCCGGTACCCCCTCCAAACGGCGCCGAGCGGGTGCGCATGGTGGCCGATGCCAAGCTGACCACCCGCCACGACCGGCAGCGCTTGACGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  2:430250/133‑1 (MQ=255)
                                                                                                                   cgggCCCAGCGCTTCGACCAGGCAGACCCGACCAAAGGCGCCGAGCTGGTGCGCATGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  8:245028/133‑1 (MQ=255)
                                                                                                                   cgggCCCAGCGCGTCGACCCGGTCGCCCCGCCCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  8:128725/133‑1 (MQ=255)
                                                                                                                   cggccccAGCGCTCCCACCTGGACCCCCCCACCACCGGCGCCGAGGTGGTGCGCAAGCTGGCCGAGGGAAAGCTCACCCCCCGCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  3:18055/131‑1 (MQ=255)
                                                                                                                   cTGGCCCTGCGCTTCTCCCCGGAAGCCCCGCCCAACGGCGCCGAGCTGGGGCGCAAGCGGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  5:447150/135‑1 (MQ=255)
                                                                                                                   cTGGCCCAGGGTTTCGACCAGCAATACCCGACAAATGGCGCCGAGGTGGTGGGAAAGCTGGCCTATGGCAAGGTTACCACCCCCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  5:23663/135‑1 (MQ=255)
                                                                                                                   cTGGCCCAGCGCTTCGACCAGGTCGACCCGACCAAAGGCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  8:241131/135‑1 (MQ=255)
                                                                                                                   cTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  8:262397/135‑1 (MQ=255)
                                                                                                                   cTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  7:318029/135‑1 (MQ=255)
                                                                                                                   cTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGACCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  7:393027/135‑1 (MQ=255)
                                                                                                                   cTGGCCCAGCGCTTCGACCAGGTACACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCTATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  8:14297/135‑1 (MQ=255)
                                                                                                                   cTGGCCCAGCGCTTCCACCAGGTAGACCCGACCAAAGCCGCCGCGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  2:33671/135‑1 (MQ=255)
                                                                                                                   cTGCCCCTGCGCTTCGCCCAGGTAGACCCCACCAAGCGCCCCGCGCTGGTGCGAAAGCTGGCCTATTGAAAGCTGACCCCCCACAAAGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  8:402484/135‑1 (MQ=255)
                                                                                                                    tGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGCGCGCCGAGCTGGGGCGCACGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGAc                  <  2:167167/135‑1 (MQ=255)
                                                                                                                             gcTTCGACCAGGTAGACCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCt        <  8:296861/135‑1 (MQ=255)
                                                                                                                               ttCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTgg      >  4:363388/1‑135 (MQ=255)
                                                                                                                                tCGACCCGGTCGACCCGCCACACGGGGCGGGGGGGGTGCGCAAGCGGGCCGAGGGCATGCTGCCCACCCTCCCCGACCGGCCGCGCGTGAAGGGGCGGGCGGCATGCCGGCCAATCCCGGACGCCCAGCCCTggg     <  4:299767/135‑1 (MQ=255)
                                                                                                                                  gACCAGGTAGACCCGACCAAAGGCGCCGAGCGGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTggg     <  8:17818/133‑1 (MQ=255)
                                                                                                                                  caCCAGGAACACCCGACCAACGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTgggg    <  7:41128/133‑1 (MQ=255)
                                                                                                                                   aCCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCg  >  4:459451/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TAGCCTGCAGGCACTGAGCCAGGAGGCCAGCCGCAGCCTCTACGAAGGCCGCACCACGGTCGGCATTCTCACCCCGGTCGGTCTGTTGGCTGGCTACAGCGCCGATGCCAGCAAGCTGGCCCAGCGCTTCGACCAGGTAGACCCGACCAAAGGCGCCGAGCTGGTGCGCAAGCTGGCCGATGGCAAGCTGACCACCCTCCACGACCGGCAGCGCTTGAAGGTGCTGGCGGCATTCCGGCCAATCCCGGACGCCCAGCCCTGGGGCG  >  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: