Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,688 A→G 8.8% intergenic (+295/+420) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,6880AG 8.8% 134.9 / 7.4 70intergenic (+295/+420)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (19/40);  new base G (2/4);  total (21/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

ACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTA  >  NC_000913/2304539‑2304835
                                                                                                                                                     |                                                                                                                                                   
aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCAGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGAAGGTCGGATAAGGCGGTTACGCCGCATCCGACAGCTAACGCCCa                                                                                                                                                     >  1:514039/1‑150 (MQ=255)
     gCGGTCTCGCCGCTTCAGACATCCATCGACCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                <  2:625424/150‑1 (MQ=11)
     gCGGTCACTCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                <  2:1120413/150‑1 (MQ=21)
     gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                >  2:2126594/1‑150 (MQ=34)
     gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                >  2:2590830/1‑150 (MQ=34)
     gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                <  2:1120470/150‑1 (MQ=34)
       ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                              >  2:1816843/1‑150 (MQ=34)
       ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                              <  2:2809289/150‑1 (MQ=34)
       ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                              <  2:2809251/150‑1 (MQ=34)
        gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                             <  2:2773736/150‑1 (MQ=34)
         tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                            <  1:1686225/150‑1 (MQ=33)
          cACGCCGCATCCGACATCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCACGCCGGTTGc                                                                                                                                           <  2:225866/150‑1 (MQ=255)
          cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                           <  2:2729892/150‑1 (MQ=32)
          cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                           <  2:2496171/150‑1 (MQ=32)
          cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                           <  2:874123/150‑1 (MQ=32)
           aCGCCGCATCCGACATCCAACGCCCGAGCCGGGTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATTAGGCGGGTGCGCCGCATCCGACAGCTAACGCCCAAGCCGGTTGcc                                                                                                                                          >  1:706664/1‑150 (MQ=11)
            cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                         >  1:353054/1‑150 (MQ=32)
              cctcctCCGACATCCTACGCGCGAGCCGGTTGCCTGATTCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       <  2:662118/145‑1 (MQ=11)
              ccggctCCGTCCTCCACCGCCCGGGGCGGTTGCGTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGCCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       <  1:2746554/145‑1 (MQ=11)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       <  2:745248/150‑1 (MQ=17)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGAGAAGGCGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       <  2:670627/150‑1 (MQ=12)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       <  1:761536/150‑1 (MQ=32)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       >  2:2389954/1‑150 (MQ=32)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       <  1:2229479/150‑1 (MQ=32)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                       <  1:2229435/150‑1 (MQ=32)
               cgcATCCGACATCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                      <  2:2467781/150‑1 (MQ=17)
               agcATCCGACATCCATCGCCACAGCCGTTCTCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                      <  1:313060/149‑1 (MQ=11)
                gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                     <  2:1773563/150‑1 (MQ=32)
                 cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTACTCCGCTGTGAAGTGCTACACCCCGTAGGGAGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTggggc                                                                                                                                    >  2:2312612/1‑146 (MQ=11)
                   tCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                                                                  >  2:1131888/1‑150 (MQ=17)
                   tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                                                                  >  2:1629587/1‑150 (MQ=32)
                   tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                                                                  >  2:1337932/1‑150 (MQ=32)
                    ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGAc                                                                                                                                 <  2:445759/150‑1 (MQ=32)
                    ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGAc                                                                                                                                 <  1:1629587/150‑1 (MQ=32)
                     cGACATCCAACGCCCGAGCCGGTTTCCTGATGCGACGCTGGCGCGTCTTATAAGGCCTACACCGCTGTGAAGTGCTCCAACCCGTAGGTCTGATAAGGCGGTTACGCCGCATACGACATCTAACGCCCAAGCCGGTTGCCTGGTGCGACg                                                                                                                                >  1:714180/1‑150 (MQ=11)
                          tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTTGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                           <  1:23174/150‑1 (MQ=17)
                          tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                           <  1:630139/150‑1 (MQ=32)
                           ccAACGCCCCAGCCGGTTGCCTGATGCGCCGCTGTCGCGACTTATAAGTCGTTCACCGGTGTTAAGTGGTGCCCACCGTAGTTCGGATAAGGCGGTTACGCCGCATCCTACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGgcg                                                                                                                          <  1:1057399/150‑1 (MQ=11)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                         <  1:399228/150‑1 (MQ=32)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                         <  1:2266886/150‑1 (MQ=32)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                         <  1:2906388/150‑1 (MQ=32)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                         <  1:2059582/150‑1 (MQ=32)
                              aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                       <  1:1013800/150‑1 (MQ=32)
                              aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                       <  1:1212774/150‑1 (MQ=32)
                                                                                                   aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg                                                  >  1:1117344/1‑150 (MQ=11)
                                                                                                   aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg                                                  >  1:2398808/1‑150 (MQ=11)
                                                                                                    ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCCGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                 <  1:165245/150‑1 (MQ=14)
                                                                                                     cccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCTCCCCCCCGTAGGGTGGATAAGGCGGTTACGCCGCATCCGAc                                                >  1:2483780/1‑150 (MQ=12)
                                                                                                       cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                              <  1:2838926/150‑1 (MQ=17)
                                                                                                         tAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATcc                                            <  1:583199/150‑1 (MQ=32)
                                                                                                           ggTCGGATTAGGCGGTTTCGCCGCATACTTCATCTAACGCCCCATCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTTCTCCGCTGTGAAGTGCTGCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                          <  2:481347/150‑1 (MQ=11)
                                                                                                           ggTCGGATAAGGCGGTTACGCCGCATCCGACATCGACCGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTTCACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                          <  2:2884419/150‑1 (MQ=255)
                                                                                                              cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGGCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                       >  1:1950707/1‑150 (MQ=25)
                                                                                                               ggATAAGGCGGTGACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTc                                      <  2:1878189/150‑1 (MQ=25)
                                                                                                                  tAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                                   <  1:1580312/150‑1 (MQ=37)
                                                                                                                      gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                               <  2:2147391/150‑1 (MQ=12)
                                                                                                                             cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCt                        >  2:293105/1‑150 (MQ=21)
                                                                                                                             cgccgaATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCt                        <  1:1482026/150‑1 (MQ=37)
                                                                                                                               ccgcATCCGTCATCTAACGCCCAAGCCGGTTGCCTGAGGCGTCGCTGGCGCGTCTTTTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGa                      <  1:737722/150‑1 (MQ=255)
                                                                                                                               ccgcATCCGACTTCTATCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGGCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGTCATCCAACGTCCGAGCCGGTTGCCTGa                      <  1:2111157/150‑1 (MQ=11)
                                                                                                                               ccgcATCCGACATCTACCGCCCAAGCCGGTTGCCTGAAGCGACGCTGGCGCGTCTTATGAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTAACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGa                      <  1:2014417/150‑1 (MQ=255)
                                                                                                                               ccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGa                      <  1:2014398/150‑1 (MQ=37)
                                                                                                                               ccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGa                      >  2:2295003/1‑150 (MQ=17)
                                                                                                                                 gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATg                    <  1:293105/150‑1 (MQ=21)
                                                                                                                                    tCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                              >  2:2793964/1‑121 (MQ=32)
                                                                                                                                     ccGACATCCAACGGCCCTGCCGGTTGCCTGATTCTACTCTGGCTCGTCTTATCTTGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGAc                <  2:922285/150‑1 (MQ=11)
                                                                                                                                      cGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGGAGGGTGGCTAAGGCGGTTACGCCGCATCCGACATCCAACGGCCGAGCCGGGTGCCTGAGGCGACg               >  2:1681015/1‑150 (MQ=11)
                                                                                                                                       gACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTCTGAAGTGCACAACCCCGTAGGTTGGATAAGGCCCTTACGACGCAACCGACATCAAACTTCCGAGCCGGATGCCTGATGCGACGc              >  2:2892627/1‑150 (MQ=11)
                                                                                                                                         cATCTAACGCCCAAGCCGGTTGCCTGATGCGGCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTg            <  2:1950707/150‑1 (MQ=37)
                                                                                                                                             caaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc        <  2:2223973/149‑1 (MQ=255)
                                                                                                                                                   ccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCAAAACACAGATGAGAAGTGCTCCACCCCATATGTTGGATAAGTAGGAAACGCCCCATCCGACATCCAACGTACGAGCCTGATGACTGAAGCGACGCTTGCGCGTCTTa  >  2:1611917/1‑150 (MQ=2)
                                                                                                                                                     |                                                                                                                                                   
ACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTA  >  NC_000913/2304539‑2304835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: