Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,688 A→G 7.3% 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 7.3% 274.9 / 10.3 122intergenic (+295/+420)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (46/53);  new base G (2/6);  total (48/74)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304541‑2304837
                                                                                                                                                   |                                                                                                                                                     
aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAg                                                                                                                                                     >  2:3135197/1‑150 (MQ=34)
aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAg                                                                                                                                                     >  1:3997307/1‑150 (MQ=34)
aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAg                                                                                                                                                     <  1:3612992/150‑1 (MQ=34)
   gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  2:838104/150‑1 (MQ=34)
   gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  2:3441520/150‑1 (MQ=34)
   gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  2:1780716/150‑1 (MQ=34)
   gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  >  2:1388096/1‑150 (MQ=34)
   gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  1:3357623/150‑1 (MQ=34)
   gCGGTAACGCCGCATCCGACATCCAACGCCCGAGCCGTTTGCCTGTTGCGAGGCTGGCGCGTCTTCTGAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  1:3816530/150‑1 (MQ=11)
      gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACTCTGGCGCGTCTTATCAGGCCTACACAGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               <  2:4144330/150‑1 (MQ=21)
      gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               <  2:3724788/150‑1 (MQ=34)
      gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               <  2:1977277/150‑1 (MQ=34)
      gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               >  1:4967741/1‑150 (MQ=34)
      gTCACGCCGCACCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATGTGGCGTACACGGCTGTGATGTGCTCGACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               <  1:5090893/150‑1 (MQ=11)
         aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                            <  1:5064504/150‑1 (MQ=32)
          cgccgcTTCCGACATCCAACGCCCGAGCCGTTTGCGTTATGCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  1:1335897/150‑1 (MQ=11)
          cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  1:1804231/150‑1 (MQ=32)
          cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           >  2:138508/1‑150 (MQ=32)
          cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  1:2950695/150‑1 (MQ=32)
          cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  2:583848/150‑1 (MQ=32)
          cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGAAGGTCGGATAAGGCGGTTACGCCGCATCCTACATCTAACGCCCAAGCCGGATGCCt                                                                                                                                           >  2:138507/1‑150 (MQ=12)
          cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAGAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCggaggccg                                                                                                                                           >  2:5144433/1‑144 (MQ=255)
           gccgcATCCGACATCAAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGAGAAGGCGGTTACGCCGCATCCGACATCCAACGACCAAGACggaagcgag                                                                                                                                          >  2:3147158/1‑143 (MQ=11)
            ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGAGGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         <  1:772713/150‑1 (MQ=17)
             cgcATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  2:1479930/150‑1 (MQ=17)
             cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:5457135/150‑1 (MQ=32)
             cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:4875351/150‑1 (MQ=32)
             cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  2:4632510/150‑1 (MQ=32)
              gcATCCGACATCCAACGCCCGAGCCGGTTGGCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  1:2778641/150‑1 (MQ=17)
              gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCTGGTTGCCTGATg                                                                                                                                       <  1:3892696/150‑1 (MQ=17)
              gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGACAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  2:625114/150‑1 (MQ=17)
                aTCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCaa                                                                                                                                                      >  1:4207182/1‑133 (MQ=32)
                aTCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCaa                                                                                                                                                      <  2:4207182/133‑1 (MQ=32)
                  ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGAc                                                                                                                                   <  1:5092089/150‑1 (MQ=32)
                   cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGGTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                                                  >  2:4398287/1‑150 (MQ=17)
                      cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAACAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTg                                                                                                                               <  2:910546/150‑1 (MQ=17)
                        tCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAGGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                             <  1:5250830/150‑1 (MQ=17)
                          cAACGCCCGAGCCGGTTGCCTTATGCGACGCTGGCGCGTATTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  1:4621102/150‑1 (MQ=17)
                          cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  1:4621361/150‑1 (MQ=32)
                          cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAACAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTAACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  1:4621135/150‑1 (MQ=17)
                           aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                          >  2:5180348/1‑150 (MQ=32)
                            ccGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                         <  1:2804833/149‑1 (MQ=17)
                            aCGCCCGAGGCGGTGGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCGGGCGCGt                                                                                                                         <  1:4754209/150‑1 (MQ=255)
                            aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGTTTGCCTGATGCGACGCTGGCGCGt                                                                                                                         <  1:4624284/150‑1 (MQ=17)
                             cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  1:3167093/150‑1 (MQ=32)
                             cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  1:1814664/150‑1 (MQ=32)
                              gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                       <  2:4931918/150‑1 (MQ=31)
                              gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                       <  1:1768431/150‑1 (MQ=31)
                              gCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                       <  2:5360908/150‑1 (MQ=14)
                                                                                                 aCCCCGTAGGTCTGAGAAGGCGGTGACGCCGCATCCGACATCGAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg                                                    <  2:4647248/150‑1 (MQ=11)
                                                                                                  ccccGTAGGTCTGACTAGGCGGTTACGCCGCATCCGACATCTAAAGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:3463192/150‑1 (MQ=11)
                                                                                                  ccccGTAGGTCTGACAAGGCGGTGACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:1875554/150‑1 (MQ=9)
                                                                                                  ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  1:5434974/150‑1 (MQ=31)
                                                                                                  ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:14164/150‑1 (MQ=31)
                                                                                                     cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                                >  2:5103043/1‑150 (MQ=32)
                                                                                                     cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                                >  1:4004657/1‑150 (MQ=32)
                                                                                                      ttAGGTCTGAGAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATc                                               <  2:1927353/149‑1 (MQ=17)
                                                                                                      gTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATc                                               >  2:1044942/1‑150 (MQ=32)
                                                                                                        acgTCCGATAAGGCGGTAACGCCGCAACCGACATCTAACGCCCAAGCCGCTAGCCTGATGCGACGCAGGCGCGACATATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                                             <  2:3666814/148‑1 (MQ=11)
                                                                                                         ggTCTGATAAGGCGGTGACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                            <  1:3545663/150‑1 (MQ=255)
                                                                                                          gTCGGATAAGGCGGTTTCGACGGATCCGCCATCTAACTCCCGAGCCGGTCGGCTGATTCGAATCTGCCTCGTCTTATGTGGCCTACACCTCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc                                           <  1:482821/150‑1 (MQ=11)
                                                                                                           tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCATACACCGATGTGAAGTGCTCCACCCCGTCGGTTTGATAACGCGGATACGCAGAAAACGACATCAAACg                                          >  2:225517/1‑150 (MQ=11)
                                                                                                            cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                         <  2:4304546/150‑1 (MQ=25)
                                                                                                             tgATATGGCGTTCACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATACGGCGGTTACGCCGCATCCGACATCCAACGTc                                        <  1:3147158/149‑1 (MQ=11)
                                                                                                                tAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                                     >  2:4435285/1‑150 (MQ=37)
                                                                                                                  aGGCGGTCACGCCGCATCCGACATCTCACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTTGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGc                                   <  1:5073379/150‑1 (MQ=11)
                                                                                                                    gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 >  2:945434/1‑150 (MQ=37)
                                                                                                                    gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGCCg                                 >  2:896533/1‑150 (MQ=21)
                                                                                                                    gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 >  2:945446/1‑150 (MQ=37)
                                                                                                                      ggTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGt                               >  2:4529055/1‑150 (MQ=37)
                                                                                                                      ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGt                               <  2:173689/150‑1 (MQ=255)
                                                                                                                        tcaCGCCGCATCCGCCTTCGTACGCCCCAGCGTGTTGCCTGTTGCGTCGCTGGCGGGCCTTATCAGGCGCTCTGCGGTGTGAAGTGCTGCACCCCGTAGGTGTGAGAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTg                             <  2:504006/148‑1 (MQ=2)
                                                                                                                         tACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:4371693/150‑1 (MQ=37)
                                                                                                                         caCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:317216/149‑1 (MQ=18)
                                                                                                                          aCGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                           <  1:1263369/150‑1 (MQ=37)
                                                                                                                            tccgcATCCGACATCTTAAGCCCCAGCGGTTTGCCTGATGCTACTCTGGCGCGTCTTATCTGGCCTACACGGCTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                         <  1:1670195/149‑1 (MQ=11)
                                                                                                                            gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGCCGGTTGCCTg                         >  2:2508721/1‑150 (MQ=21)
                                                                                                                            gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                         >  1:4337914/1‑150 (MQ=37)
                                                                                                                            gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                         >  2:595862/1‑150 (MQ=37)
                                                                                                                             ccgcATCCGACATCTAACGCCCTCGCCGGTTGCCTGTTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGa                        <  1:590841/150‑1 (MQ=255)
                                                                                                                             ccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGa                        >  1:4359825/1‑150 (MQ=37)
                                                                                                                             ccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCgg                                >  1:1358973/1‑142 (MQ=17)
                                                                                                                              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt                       >  2:3703260/1‑150 (MQ=21)
                                                                                                                                 aTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCAAACGTCCGAGCCGGTTGCCTGATGCg                    >  2:4538769/1‑150 (MQ=21)
                                                                                                                                    cGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGAGTAGTGCTCCACCCCGTAGGTAGGATAAGGCGGTTGCGCAGCAACCGACATCGAACGGCCGAGCCGGTTGCCTGATGCGACg                 >  2:3172960/1‑150 (MQ=11)
                                                                                                                                     gACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGc                >  2:3691327/1‑150 (MQ=37)
                                                                                                                                     gACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTAGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGc                >  2:3339085/1‑150 (MQ=37)
                                                                                                                                          ccaaCGCCCTAGCCGGTTGCCTGAGGCGACCCTGGCGCGGCGGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcg           <  1:5407567/148‑1 (MQ=11)
                                                                                                                                          caaaCGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTTGACAAGGCGGTAACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcg           <  2:2883568/148‑1 (MQ=11)
                                                                                                                                           tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  2:1463038/1‑150 (MQ=37)
                                                                                                                                           tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGGTGGATAAGGCGGTTACGCCGCATCCGACATCCAACTTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  1:13499/1‑150 (MQ=255)
                                                                                                                                            aaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCGGGcgcg         >  1:2763569/1‑150 (MQ=255)
                                                                                                                                             aCGCCCATGGCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt        <  2:3023826/150‑1 (MQ=255)
                                                                                                                                               gCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGAAACCCCCCATCAGGAGTGCGCCACCCCGTAGGTAGGAAAAGGCGGTTGCGCCGCATCCGACATCAAACGACCGATCAGGTAGCCTGATGCGACGCTGGAGCGTCt      >  2:2597327/1‑150 (MQ=2)
                                                                                                                                                cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCtt     <  1:3211554/150‑1 (MQ=37)
                                                                                                                                                cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCtt     <  2:3485890/150‑1 (MQ=17)
                                                                                                                                                 ccACGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTACCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACACCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa    <  2:1162326/150‑1 (MQ=255)
                                                                                                                                                 ccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATAAGGCCTACACCCCTGTGAAGTGCTACACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAACCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa    >  1:4303507/1‑150 (MQ=11)
                                                                                                                                                  cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAACCAACGACCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt   >  1:3872413/1‑150 (MQ=21)
                                                                                                                                                   aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc  <  1:3735241/150‑1 (MQ=37)
                                                                                                                                                   |                                                                                                                                                     
AAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304541‑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: