Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,688 A→G 8.4% 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.4% 303.8 / 14.8 135intergenic (+295/+420)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (38/70);  new base G (2/8);  total (42/93)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

CAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304540‑2304837
                                                                                                                                                    |                                                                                                                                                     
 aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAg                                                                                                                                                     >  2:45042/1‑150 (MQ=34)
  aGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTATTATCAGGCCTACACCGCTTTGAAGTGCTCCACCCCGTTGGTGGTATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGc                                                                                                                                                    <  1:718617/150‑1 (MQ=11)
    gCGGTCACGCCGCATCCGACATCCAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTTTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  1:4625066/150‑1 (MQ=21)
    gCGGTCACGCCGCATCCGACATACAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAACCGCTGTGAAGTGCTCCACCCCGTAGGTCGTATAAGGCGGTTACGCCGCATCCGAAATCTAACGCCCAAGCCg                                                                                                                                                  <  1:2809064/150‑1 (MQ=255)
      ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                                       >  2:4624530/1‑143 (MQ=34)
       gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               <  1:1126095/150‑1 (MQ=34)
        tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                              <  1:4227435/150‑1 (MQ=33)
        tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                              <  1:4906789/150‑1 (MQ=33)
        tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                              <  1:3772249/150‑1 (MQ=33)
        tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                              <  2:1720095/150‑1 (MQ=33)
        tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                              <  2:1720169/150‑1 (MQ=33)
        tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCATATCAGGCCTACACCGCGGTGAAGAGCTCCCCCCCGTAGGTCGGCTAAGGCGGGTACGCCGCATCCGACATCGAACGCCCCCGCCggggg                                                                                                                                              >  1:3547975/1‑147 (MQ=11)
         cACGCCGCATCCGACATCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAGAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  2:4134229/150‑1 (MQ=17)
         cACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:3546231/150‑1 (MQ=12)
         cACGCCGCATCCGACATCCAACGCCCGAGCGGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  2:3580794/150‑1 (MQ=17)
         cACGCCGCATCCGACATCCAACGCCCGAGCCGTTTTCCTGATGCGACGCTGGCGCGTCTTATCTGGGCTACTGCTGTTTTAAGTGGTGCTCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGc                                                                                                                                             <  1:3048575/150‑1 (MQ=11)
          aCGCCGCATCCGACATCCAACGCGCGAGCGGTTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                            <  2:5093786/150‑1 (MQ=12)
           cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCt                                                                                                                                           >  2:1300347/1‑150 (MQ=17)
           cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  1:36922/150‑1 (MQ=32)
             ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCGAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         >  1:2991543/1‑150 (MQ=17)
              cgccgccGACATCCCCCGCCCGAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCGTCTCAGGCCTACACCGCTGTGACGTGCTCCACCCCGTAGGTCGGATACGGCGGTTACGCCGCCTCCGACATCTAACGCCCCCGCCGGTTGCCTGAt                                                                                                                                        <  2:4016834/145‑1 (MQ=11)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGCTGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  2:1061778/150‑1 (MQ=17)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:4460933/150‑1 (MQ=32)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  2:2490141/150‑1 (MQ=32)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:1358179/150‑1 (MQ=32)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:928754/150‑1 (MQ=32)
              cgcAGCCGACAGCCCCCGCCCGAGCCGGGGGCCTGCTGCGCCGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGACGGGCGCCACCCCGTAGGTCGGAGAAGGCGGTGACGCCGCAGCCGACATCTAACGCCCCCGCCGGGTGCCTGAt                                                                                                                                        <  2:1062009/150‑1 (MQ=11)
               gcctCCGACCTCCAACGCCCGAGCCGGTTGCCTGCTGCGACGCTGGCGCGACCTATCAGGCCTACACCGCGGTGAAGCGCTCCACGCCGTAGGTCTGACAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  2:1081569/147‑1 (MQ=11)
               gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  2:3570544/150‑1 (MQ=32)
                cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGc                                                                                                                                      <  1:2216065/150‑1 (MQ=32)
                cATCAGACATCCAAAGCCGGAGCCGGTTGCCTGCTGAGTCGTTGTCGAGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGc                                                                                                                                      <  1:2216134/150‑1 (MQ=11)
                  tCCGACATCCCACGCCCGAGCGGGTTGCCTGATGCGACGCTGGCTCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                                                                    <  2:4883231/150‑1 (MQ=255)
                  tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                                                                    <  2:4883238/150‑1 (MQ=32)
                    cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                                                  >  2:4275454/1‑150 (MQ=32)
                    cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                                                  <  1:5263688/150‑1 (MQ=32)
                     gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                                                                 >  2:4848850/1‑150 (MQ=32)
                     gACATCCAACGCCCGAGCCGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                                                                 >  2:5468526/1‑150 (MQ=17)
                       cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTg                                                                                                                               <  1:3908600/150‑1 (MQ=32)
                        ctCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGACTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCTTAGGTCGGATAAGGCGGTTACGCCGGATCCGACTTCTAACGCCCAATCCGGTTGCCTGATGCGACGCTgg                                                                                                                              <  1:4659528/149‑1 (MQ=11)
                        aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                              <  1:2271784/150‑1 (MQ=32)
                        aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                              <  1:5690349/150‑1 (MQ=32)
                           cAACGCGCGAGCCGGTTGCCTTATGCGAAGCTGGCTCGTCTTATCATGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  2:1687822/150‑1 (MQ=11)
                            aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                          <  1:4275483/150‑1 (MQ=32)
                            aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                          <  1:4275454/150‑1 (MQ=32)
                             aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                         <  1:2964850/150‑1 (MQ=32)
                             aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                         <  2:2858771/150‑1 (MQ=32)
                              cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  1:2373429/150‑1 (MQ=32)
                              cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  1:2099035/150‑1 (MQ=32)
                              cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  1:4167889/150‑1 (MQ=32)
                              cGCCCGAGCCGCTTGCCTGATGCGAAGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  1:2099882/150‑1 (MQ=17)
                                                                                                   ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   >  1:1145357/1‑150 (MQ=31)
                                                                                                     ccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACa                                                 <  1:5622335/150‑1 (MQ=32)
                                                                                                     ccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACa                                                 <  1:2332626/150‑1 (MQ=32)
                                                                                                      cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCATGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTGTGAAGTGCTCCACCCCGTGGGGTGGATAAGGCGGTTACGCCGCATCCGACag                                                >  1:5420354/1‑149 (MQ=255)
                                                                                                      cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAGACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGGTAAGGCGGTTACGCCGCATCCGACag                                                >  1:655156/1‑149 (MQ=17)
                                                                                                      cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                                >  1:5382487/1‑150 (MQ=32)
                                                                                                       gTAGGTCGGATAAGGCGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATc                                               <  2:2361420/150‑1 (MQ=17)
                                                                                                       gTAGGTCCGATAAGGCGGTCACGCCGCCTCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATc                                               <  2:2713850/150‑1 (MQ=255)
                                                                                                           gTCGGATTAGGCGGTGACGCCGCATCCGACATCTCACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATGAGGCGGTTACGCCGCATCCGACATCCAAc                                           <  2:4203276/150‑1 (MQ=11)
                                                                                                            tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  1:3687678/1‑150 (MQ=34)
                                                                                                            tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  2:21196/1‑150 (MQ=34)
                                                                                                            tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  2:21220/1‑150 (MQ=34)
                                                                                                            tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCGGTGAAGTGCTCCACCCCGGAGGTTGGATAAGGCGGTGACGCCGCATCCGACATCGAACg                                          >  1:2360438/1‑150 (MQ=255)
                                                                                                            tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGACCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  2:2582524/1‑150 (MQ=21)
                                                                                                            tCGGATAAGGCGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  2:5245847/1‑150 (MQ=21)
                                                                                                             cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                         >  2:3904769/1‑150 (MQ=35)
                                                                                                              ggAGAAGGCGGTTACGCCGCATCCGACATCCTACGCCCACGCCGGTTGCGTGAGTCGACGCTGGCGCGTCTTATCTGGGCTACACCGCCGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTc                                        <  2:1180021/150‑1 (MQ=11)
                                                                                                               gATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTcc                                       >  1:2718728/1‑150 (MQ=37)
                                                                                                               gATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCACAACACCCATCAGACCTGCTCAACCCCCAAGGCTGGATAACGACGTTACGCCGCATCCGGCATCAAACGGcc                                       >  2:3227721/1‑150 (MQ=11)
                                                                                                                 tAAGGCGTTTACGCCGCATCCTACATCTAACGCCCCAGCGGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                                     <  1:3092153/150‑1 (MQ=255)
                                                                                                                     gCGGTTACGCCGCATCCGACATCTAACGCGCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATGAGGGCTACACCGGTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  2:1823544/150‑1 (MQ=11)
                                                                                                                        gTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGCCGGtt                              >  1:5620309/1‑150 (MQ=21)
                                                                                                                         ttACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTg                             <  1:1224610/150‑1 (MQ=37)
                                                                                                                         tcaCGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTg                             <  1:5623113/148‑1 (MQ=37)
                                                                                                                          tACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTTATGCGACTCTGGCGCTTCTTATCAGGAATACACCGCTGTGATTTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  1:5571405/150‑1 (MQ=11)
                                                                                                                          tACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  1:4281878/150‑1 (MQ=37)
                                                                                                                            cgccgcCTCCGACAGCGAACGGCCTCGCCGGTTGCCTGCTGCGGCGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCt                          <  1:2787531/150‑1 (MQ=11)
                                                                                                                            cgccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCt                          >  1:3373155/1‑150 (MQ=37)
                                                                                                                              ccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACGTCCGAGCCg                                 >  1:4040744/1‑141 (MQ=17)
                                                                                                                               cgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt                       >  1:1810827/1‑150 (MQ=37)
                                                                                                                                  aTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCg                    <  2:2376062/150‑1 (MQ=37)
                                                                                                                                   tgcGACACCGAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGa                   <  1:793755/148‑1 (MQ=255)
                                                                                                                                   tCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGa                   <  1:3509763/150‑1 (MQ=37)
                                                                                                                                     cGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCAAACGACCGAGCCGGTTGCCTGATGCGACg                 >  2:3704050/1‑150 (MQ=17)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  2:691345/1‑150 (MQ=37)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  1:1829152/1‑150 (MQ=37)
                                                                                                                                             aaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg         >  2:4987901/1‑150 (MQ=37)
                                                                                                                                             aaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGCTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg         >  2:3845142/1‑150 (MQ=255)
                                                                                                                                             aaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGAGAAGTGCGACACCCCGTAGGTAGGATAAGGCGGTTACGCCGCATACGACAACCAACGACCGAGCCGGATGCCTGATGCGACGCTGGcgcg         >  1:3924420/1‑150 (MQ=11)
                                                                                                                                              aCGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt        <  2:1966230/150‑1 (MQ=37)
                                                                                                                                               cGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTc       <  1:3451407/150‑1 (MQ=37)
                                                                                                                                                gCCCTATCCGGGTGGCTGATGAGACGCTGGCGCGTCTTATCAGGCCTACTCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt      <  2:785049/150‑1 (MQ=11)
                                                                                                                                                gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt      <  1:4041256/150‑1 (MQ=37)
                                                                                                                                                 cccTAGCCGGTTGCCTGATGCGACGCTGGCGCCTCGGATGAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTAGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCtt     <  2:1217330/150‑1 (MQ=11)
                                                                                                                                                   cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt   <  1:435518/150‑1 (MQ=37)
                                                                                                                                                   cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt   >  2:4003253/1‑150 (MQ=37)
                                                                                                                                                    aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCGAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATc  >  2:3268430/1‑150 (MQ=21)
                                                                                                                                                    aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCGAACGACCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc  >  2:3268416/1‑150 (MQ=17)
                                                                                                                                                    aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc  >  2:1898105/1‑150 (MQ=17)
                                                                                                                                                    aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAAGCCGCAACCCACATACAACGACCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc  >  1:4895046/1‑150 (MQ=11)
                                                                                                                                                    |                                                                                                                                                     
CAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304540‑2304837

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: