Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,688 A→G 13.1% 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,6880AG13.1% 338.0 / 43.1 183intergenic (+295/+420)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (44/90);  new base G (5/16);  total (52/131)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 193 total aligned reads displayed.

CAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304540‑2304837
                                                                                                                                                    |                                                                                                                                                     
cAAGTCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCaa                                                                                                                                                      <  2:720324/150‑1 (MQ=25)
cAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCaa                                                                                                                                                      <  2:4627140/150‑1 (MQ=35)
cAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCaa                                                                                                                                                      <  2:6764144/150‑1 (MQ=35)
 aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAg                                                                                                                                                     >  1:5677783/1‑150 (MQ=34)
  aGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGc                                                                                                                                                    >  2:5440747/1‑150 (MQ=12)
    gCGGTCACGCCGCATCCGACATCCATGGGCCGTGGCGTTTGCCTGTTGGGTCTCTGGAGCGTCTTTTGAGGCCTAAACTGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAAATCTAACGCCCAAGCCg                                                                                                                                                  <  1:2617730/150‑1 (MQ=11)
    gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCg                                                                                                                                                  >  2:1190488/1‑150 (MQ=21)
    gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  >  1:5168782/1‑150 (MQ=34)
    gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAACCGCTGTGAAGTGCTCCACCCCGTTGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  1:747974/150‑1 (MQ=21)
    gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGAAGCTGGCGCGTCTTATCAGCACTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAGCCGACGTCTGACGCCCAAGCCg                                                                                                                                                  >  1:5246319/1‑150 (MQ=11)
    gCGGTCACGCCGCATCCGACAGCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  1:6584268/150‑1 (MQ=21)
      ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCTTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGTATAAGGCGTTTACGCCGAATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                                <  1:2016277/150‑1 (MQ=255)
      ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                                >  2:1654710/1‑150 (MQ=34)
      ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                                >  2:3788492/1‑150 (MQ=34)
       gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               <  1:8108516/150‑1 (MQ=34)
       gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                               <  1:6414412/150‑1 (MQ=34)
         cACGCCGCTTCCGACATCCTACGCCCGAGCGGGTGGCGTTAGTCTACTCTGGCTCGTCTTATGTGGCCTACACCGCCGTGAAGTGCTGCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:5343469/150‑1 (MQ=11)
         cACGCCGCCTCCGACATCCCACGCCCGAGCCGGTTGCCTGTTGCGTGGCTGGCGCGTGTTATCAGGCCTTCCCCGCTGTGTTTCGCTCCACGCCGTTGGTCGGATATGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  2:5246319/150‑1 (MQ=11)
         cACGCCGCCACCGACATCCAACGCCCGAGCCGGTTGGATGATGCGCCGGTGGCGCGTCTTATCATGCCTTCCCCTCTGTGTTTTGCTCCACGCCGTTGGTCGGATATGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  2:5168785/150‑1 (MQ=11)
         cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATTCGAAGCTGTCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGc                                                                                                                                             <  2:2831563/150‑1 (MQ=255)
         cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  2:5832469/150‑1 (MQ=32)
           cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  2:6892116/150‑1 (MQ=32)
           cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  2:7377671/150‑1 (MQ=32)
            gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTg                                                                                                                                          <  1:1475563/150‑1 (MQ=32)
            gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGACCTAGCCGGTTGCCTg                                                                                                                                          >  2:6213492/1‑150 (MQ=17)
             cctcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         <  2:2269617/147‑1 (MQ=17)
             ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGGGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         <  2:7827239/150‑1 (MQ=17)
             ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         <  1:4317177/150‑1 (MQ=32)
             ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         >  1:2964261/1‑150 (MQ=32)
             ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCAGATGCGACGCTGGCGCGTATTATAAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTAGGCTAAGGCGGTTACGCCGCATCCGACAGCTAACCCCCACGCCCGTTGCCTGa                                                                                                                                         >  1:6464571/1‑150 (MQ=11)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:1542962/150‑1 (MQ=32)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:731500/150‑1 (MQ=32)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:7939956/150‑1 (MQ=32)
              cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  2:1669427/150‑1 (MQ=32)
               gcATCCGACATCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  2:4595809/150‑1 (MQ=17)
               gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCAACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCCAGCCGGTTGCCGGATg                                                                                                                                       >  1:8077604/1‑150 (MQ=11)
               gcATACTACATACAAAGCCCGAGCCGTTTGCCTGATTCTACTATGTCTAGTCTTATCATGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  2:2381119/150‑1 (MQ=11)
                  tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                                                                    <  2:2988142/150‑1 (MQ=32)
                  tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACCCTGGCGCGTCTTATAAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGTTAAGGCGGTTACGCCGCCTCCGACCTCTAACGCCCAAGCCGGGTGCCTggggcgg                                                                                                                                    >  1:4216402/1‑144 (MQ=11)
                     gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                                                                 >  2:3446533/1‑150 (MQ=32)
                     gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGACTGATGCGActc                                                                                                                                 >  2:6935562/1‑148 (MQ=17)
                     gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCAACACCGCTGTGAAGTGCTCCACCCCCTAGGTCGGATAAGGCGGTTACGCCGCATCCGACAACAAAAGCCCAAGCCGGTTGCCAGAAGCGACGc                                                                                                                                 >  1:2776313/1‑150 (MQ=11)
                     gACATCCAACGCCCGAGCCGGGTGCCTGATGCGACGCGGGCGCGTCTTATAAGGCCTACACCCCTGTGAAGTGGGCCACCCCGGAGGGCGGATAAGGCGGTTAAGCCGCATACGACATCAAACGCCCAAGCCGGTTTCCTGGTGCGgcgc                                                                                                                                 >  1:2153262/1‑150 (MQ=14)
                        aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACCGCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                              >  2:4084425/1‑150 (MQ=12)
                         tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                             >  2:4798215/1‑150 (MQ=32)
                            caCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCAACACCGCTGTGAAGGGCTCCACCCCGTAGGTCGGAGAAGGCGGTTACGCCGCATCGGACAGCTCACGCCGAAGCCGGGTGCCGGATGCGTCGGGGGcgcg                                                                                                                          >  1:4572640/2‑150 (MQ=11)
                            aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                          >  2:3417838/1‑150 (MQ=32)
                              cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  1:4798215/150‑1 (MQ=32)
                               gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                       <  1:1230224/150‑1 (MQ=31)
                                                                                                  aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg                                                    <  2:5955067/150‑1 (MQ=11)
                                                                                                   ccccGTAGGTCCGACAAGGCGGTTACGCCGCATCCGACATCTACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:2125462/150‑1 (MQ=11)
                                                                                                    cccGTAGGTCGGATTAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGAc                                                  <  2:6464571/150‑1 (MQ=17)
                                                                                                         aGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                                             <  1:3417325/150‑1 (MQ=33)
                                                                                                          ggTCTGATAAGGCGGTCACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                            <  1:1590682/150‑1 (MQ=255)
                                                                                                          ggTCCGATAAGGCGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                            <  2:6384159/150‑1 (MQ=21)
                                                                                                           gTCTGATAAGGCGGTTACGCCGCATCCGACATCAAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc                                           <  2:7587091/150‑1 (MQ=21)
                                                                                                           gTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc                                           <  1:1027474/150‑1 (MQ=34)
                                                                                                               gATAAGGCGGTTACGCCGCATCCGACATCTATCGCCCTTGCCGTGTGCCTGTTGCTACGCTGGCGCGTCTTATGAGGCCTAAACCGCTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTAACGCCGCATCCGACATCCAACGTcc                                       <  1:1219065/150‑1 (MQ=11)
                                                                                                                 tAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCGGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                                     <  1:3788492/150‑1 (MQ=37)
                                                                                                                 tAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGTTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                                     <  1:6900819/150‑1 (MQ=37)
                                                                                                                 tAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                                     <  1:1654710/150‑1 (MQ=37)
                                                                                                                   aGGCGGTCACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTTGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGc                                   <  1:7616802/150‑1 (MQ=255)
                                                                                                                    ggCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGcc                                  >  2:1581594/1‑150 (MQ=21)
                                                                                                                     ggtgTTTCGCTGCATCCTACTTCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  1:7421121/147‑1 (MQ=11)
                                                                                                                     gCGGTTACGCCTCATCCTACATATAACGCCCTAGCGGTTTGCCTTATTCTAATCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  1:1437194/150‑1 (MQ=11)
                                                                                                                     gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  2:8077575/150‑1 (MQ=37)
                                                                                                                     gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGACTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  1:7484331/150‑1 (MQ=37)
                                                                                                                     gCGGTCACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  1:7851870/150‑1 (MQ=255)
                                                                                                                     gCGGTAACGCCGCATCCGACATCTAACGCCGCAGCCGGTTTCCTGATTCGTCGCTGGGTCTTCTTATGTGTCCTAGTCCGATGTGAAGTGCTGCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  1:6049666/150‑1 (MQ=11)
                                                                                                                      cGGATACGCCGCTCCCGAAATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTTGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCgg                                <  1:7013121/150‑1 (MQ=11)
                                                                                                                        gtcaCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                              <  1:2464916/147‑1 (MQ=12)
                                                                                                                        gtcaCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                              <  1:193952/147‑1 (MQ=12)
                                                                                                                        gTTACGCCGCATGCGACATCTAACTCCCGAGCGGGTTGCCTGATTCGAAGTTGTCTAGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                              <  1:4374862/150‑1 (MQ=11)
                                                                                                                         ttACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTg                             <  2:2117026/150‑1 (MQ=37)
                                                                                                                           aCGGCGCATCCGACATCTTACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                           <  1:7155910/150‑1 (MQ=255)
                                                                                                                           aCGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                           <  2:4994082/150‑1 (MQ=37)
                                                                                                                             gtcgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                         >  2:1011749/3‑150 (MQ=37)
                                                                                                                             gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                         <  2:3163409/150‑1 (MQ=37)
                                                                                                                             gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                         >  2:2080176/1‑150 (MQ=37)
                                                                                                                              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGa                        <  1:6093400/150‑1 (MQ=255)
                                                                                                                               cgcATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt                       <  1:6220914/150‑1 (MQ=37)
                                                                                                                               cgcATCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt                       <  2:4754410/150‑1 (MQ=37)
                                                                                                                               cgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   >  1:3747921/1‑122 (MQ=31)
                                                                                                                                tcATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACTCCGCTGTGAAGTGCTCCACCGCGTAGGTTTTATTAGGCGGTTACGCCTCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATg                      <  1:4932827/149‑1 (MQ=11)
                                                                                                                                gcATCCGTAATCTAACGCCCAAGCCGGGTGCCTGATGCGACGCTGGCGCGCCTTACCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATg                      <  2:8093025/150‑1 (MQ=11)
                                                                                                                                    ccGACATCTAACGCCCTAGCCGGTTGCCTGTTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGTTATGGCGGTTTCGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGAc                  <  2:2522798/150‑1 (MQ=11)
                                                                                                                                         aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTgg             <  2:3487847/150‑1 (MQ=12)
                                                                                                                                          tCTAACGCCCCCGCGGGTTGCCTGTTGCGACGGTGGCGCGTATTTTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTTGGTTGGATATGGGGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATAAGACGcgcca            <  1:3430821/150‑5 (MQ=11)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          <  2:3342832/150‑1 (MQ=37)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          <  1:904087/150‑1 (MQ=37)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          <  1:3501331/150‑1 (MQ=37)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGGCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  2:453547/1‑150 (MQ=21)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGc                            >  1:8178756/1‑132 (MQ=17)
                                                                                                                                            tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCAAACGACCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  2:4649446/1‑150 (MQ=17)
                                                                                                                                              ccGCCCCAGCCGGGTGCCTGCTGCGCCGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt        <  1:301977/149‑1 (MQ=11)
                                                                                                                                              aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt        <  2:5067698/150‑1 (MQ=37)
                                                                                                                                               cGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGGTACGCCGCATCCGACATCCAACGGCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTc       >  2:3685937/1‑150 (MQ=255)
                                                                                                                                                gCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGGTGGATCAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt      >  2:3031907/1‑150 (MQ=255)
                                                                                                                                                   cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCTCCACCCCGTAGGTGGGAGAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGGTGCCTGAGGCGACGCGGGcgcggcggcg   >  2:7325978/1‑144 (MQ=11)
                                                                                                                                                    aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc  >  2:2118957/1‑150 (MQ=37)
                                                                                                                                                    |                                                                                                                                                     
CAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304540‑2304837

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