Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,688 A→G 13.0% 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.0% 194.1 / 22.7 101intergenic (+295/+420)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (29/50);  new base G (3/9);  total (32/69)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

ACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304539‑2304837
                                                                                                                                                     |                                                                                                                                                     
acacGGCGGTAACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACACGTCTTTGATGTGCTCCACCCCTTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                                       <  1:2790821/150‑1 (MQ=11)
aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCAGATGCGACGCTGGCGCGTCTAATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGAAGGGAGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                                       >  1:3257956/1‑150 (MQ=11)
  aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAg                                                                                                                                                     >  2:4075979/1‑150 (MQ=34)
     gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGATGTGCTCCACCCCGTTGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  1:3572190/150‑1 (MQ=21)
     gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  2:1923776/150‑1 (MQ=34)
     gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                                  <  1:583834/150‑1 (MQ=34)
         tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                              <  1:3709034/150‑1 (MQ=33)
          cACGCCGCCTCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGACTACAGTGCTGTGTTGTGCTCTACCGCGTTGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:3807710/150‑1 (MQ=11)
          cACGCCGCATCCGACATCCTACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:2544027/150‑1 (MQ=17)
          cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:2543998/150‑1 (MQ=32)
          cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:1894832/150‑1 (MQ=32)
          cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:427371/150‑1 (MQ=32)
          cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                             <  1:427427/150‑1 (MQ=32)
           aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCAACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                            >  1:2924965/1‑150 (MQ=17)
           aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGAGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                            <  1:4094442/150‑1 (MQ=17)
            cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                           <  1:3695294/150‑1 (MQ=32)
             gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTg                                                                                                                                          >  2:488543/1‑150 (MQ=32)
             gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCAGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTg                                                                                                                                          >  1:1633877/1‑150 (MQ=17)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGATGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         <  1:2007228/150‑1 (MQ=17)
              ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                         >  1:2465556/1‑150 (MQ=17)
               cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                        <  1:4021981/150‑1 (MQ=32)
                gcttCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  2:3194341/147‑1 (MQ=17)
                gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  1:1496945/150‑1 (MQ=32)
                gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                       <  2:3194306/150‑1 (MQ=32)
                 cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCTGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGc                                                                                                                                      <  1:971591/150‑1 (MQ=17)
                  aTCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCg                                                                                                                                     <  1:1397336/150‑1 (MQ=17)
                   tCCGACCTCCCCCGCCCGTGCGGGTGGCCTGATGCGTCGGGGGCGCGTCTTATCAGGCCTACACCGCTGTGATGTGCTCCACCCGGTAGGGCGGATACGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                                                                    <  1:2474385/150‑1 (MQ=11)
                      gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                            <  1:2739344/139‑1 (MQ=32)
                      gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                            >  2:2739344/1‑139 (MQ=32)
                      gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTAATCAGTCCTACAGCGCTGTGAAGTGCTCCACCCCGTAGGTAGGAAAAGGCGGTTACGCCGAATCAGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                                                                 >  2:4024397/1‑150 (MQ=11)
                         aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                              <  2:2146477/150‑1 (MQ=32)
                          tacAACGCCCGAGCCGGTTGCCTTATTCTACTATGGCCCGTCTTATCAGGCCTACACCGCCGTGAAGTGCTCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCAAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                             <  1:492296/148‑1 (MQ=11)
                          tCCAACGCCCGAGCCGTTTGCCTGATCCTACTATGGCTCGTCTTATCAGGCCTAAACCGCTGTGAAGTGCTCCACCCCGTAGGTCCGATAAGGCGGTAACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                             <  1:492291/150‑1 (MQ=11)
                          tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                             <  1:3314712/150‑1 (MQ=32)
                          tCCAACGCCCGAGCCGGGTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                             <  1:3860620/150‑1 (MQ=17)
                          tCCAACGCCCGAGCCGGGTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                             <  1:2457344/150‑1 (MQ=17)
                           ccttcGCCCGAGCCGGTTGCCTGATGCGACGCGGGCGCGTCGTATCCGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcg                                                                                                                            <  1:2586407/146‑1 (MQ=11)
                           ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcg                                                                                                                            <  1:2586837/150‑1 (MQ=32)
                           ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTATTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcg                                                                                                                            <  1:1242205/150‑1 (MQ=17)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  1:3171634/150‑1 (MQ=32)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  1:1531424/150‑1 (MQ=32)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  1:1531535/150‑1 (MQ=32)
                            cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                           <  1:1531557/150‑1 (MQ=32)
                              aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                         <  1:3280493/150‑1 (MQ=32)
                               cGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                                                        <  2:1306521/150‑1 (MQ=32)
                                                                                                   aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg                                                    <  1:1753004/150‑1 (MQ=11)
                                                                                                   aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg                                                    <  2:2401039/150‑1 (MQ=11)
                                                                                                   aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCg                                                    >  2:2414565/1‑150 (MQ=11)
                                                                                                    ccccGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:2815688/150‑1 (MQ=9)
                                                                                                    ccccGTAGGTCGGCTAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:1633877/150‑1 (MQ=14)
                                                                                                    ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:1993046/150‑1 (MQ=14)
                                                                                                    ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGGAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:2924965/150‑1 (MQ=14)
                                                                                                    ccccGTAGGTCGGAGAAGGCGGTGACGCCGCATCCGACATCCAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa                                                   <  2:3755279/150‑1 (MQ=9)
                                                                                                      ccGTAGGTCCGATAAGGCGGTGACGCCGCATCCGACATCGAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACa                                                 <  2:2159185/150‑1 (MQ=255)
                                                                                                        gTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATc                                               <  2:3997308/150‑1 (MQ=17)
                                                                                                          aGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                                             <  1:3563607/150‑1 (MQ=33)
                                                                                                             tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  1:1754451/1‑150 (MQ=34)
                                                                                                             tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  1:1289082/1‑150 (MQ=34)
                                                                                                              cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                          >  1:3062874/1‑149 (MQ=34)
                                                                                                                gATTAGGCGGTCACGCCTCATCCGTCATCTAACGCCCTAGCGGGTTGCCTGATTCTACGCTGTCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATATGGCGGTTACGCCGAATCCGACATCCAACGTcc                                       <  1:2829835/150‑1 (MQ=11)
                                                                                                                gATATGGCGGTTACGCCGCATCCGAAATCTACCGCCCACGCCGGTTGCCTGTTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTTGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTcc                                       <  1:1684853/150‑1 (MQ=11)
                                                                                                                   aaGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAg                                    >  1:1115285/1‑150 (MQ=37)
                                                                                                                    aGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGc                                   <  2:3955694/150‑1 (MQ=37)
                                                                                                                     ggCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGcc                                  >  2:3569501/1‑150 (MQ=21)
                                                                                                                     ggCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGAAGGTTGGATAAGGCGGTTACGCCGCTTCCGACATCCAACGTCCGAGcc                                  >  1:335875/1‑150 (MQ=255)
                                                                                                                      gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGCCg                                 >  2:3166935/1‑150 (MQ=21)
                                                                                                                         gTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                              >  1:1897333/1‑150 (MQ=37)
                                                                                                                         gTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                              <  1:2414565/150‑1 (MQ=37)
                                                                                                                         gTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                              >  2:2644119/1‑150 (MQ=37)
                                                                                                                         gTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGAAAAGGCGGTTACGCCGCATCCGACATCCAACCGCCGAGCCGGtt                              >  1:874671/1‑150 (MQ=255)
                                                                                                                           tACGCCGCATCCGTCATCTACCGCCCTCGCCGGTTGCCTTATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:2543787/150‑1 (MQ=11)
                                                                                                                           tACGCCGCATCCGCCATCTCTCGCGCTAGCCGGTTGCGGGATTCTACGCTGGCGCGTCTTTAGATGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:1017865/150‑1 (MQ=11)
                                                                                                                           tACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:2837759/150‑1 (MQ=37)
                                                                                                                           gaCGCCGCATCCGACTTCGATCGCCCGAGCCGGTTGCCTTATGCGACGCTGGCGCGTCTTTACATGGCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:3552686/149‑1 (MQ=11)
                                                                                                                           caCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:1148158/149‑1 (MQ=18)
                                                                                                                           caCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                            <  2:3552661/149‑1 (MQ=18)
                                                                                                                            aCGCAGCATCCGACATCAAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTTTGATTAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                           <  2:1055587/150‑1 (MQ=11)
                                                                                                                              gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGGCCGAGCCGGTTGCCTg                         >  2:3714592/1‑150 (MQ=21)
                                                                                                                                cgcATCCGTCATCTAACGCCCTCGACGGTTGCCTGATGCTACGCTGTCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt                       <  1:2552926/150‑1 (MQ=11)
                                                                                                                                cgcATCCGACATCCCCCGCCCCAGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCGGGAGGTGGGAGAAGGCGGTTACCCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGAt                       <  2:236938/150‑1 (MQ=11)
                                                                                                                                    tCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGa                   <  1:3506198/150‑1 (MQ=37)
                                                                                                                                    tCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGa                   <  1:3506241/150‑1 (MQ=37)
                                                                                                                                        aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGACTACACCGCTCGGGATTGCTCCACCCCGTAGGTTGGATAGGGCGGTTACGCCGCATCCGACATCGAACGGCCGAGCCGGTTGCCTGTTGCGACGCt               >  2:1994963/1‑150 (MQ=11)
                                                                                                                                           tCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgc            >  2:3280662/1‑150 (MQ=37)
                                                                                                                                            ccaaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcg           <  2:3053293/148‑1 (MQ=25)
                                                                                                                                            cTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcg           <  1:1597277/150‑1 (MQ=37)
                                                                                                                                            cTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTACACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCGAACGTCCGAGCCGGTTGCCTGATGCGTCGCTGgcg           >  2:1305760/1‑150 (MQ=255)
                                                                                                                                             tAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          <  2:1145307/150‑1 (MQ=37)
                                                                                                                                             tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  1:284935/1‑150 (MQ=37)
                                                                                                                                             tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCt                          >  1:3692029/1‑134 (MQ=17)
                                                                                                                                             caaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          >  2:2608363/2‑150 (MQ=14)
                                                                                                                                             caaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGAGAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc          <  2:3158304/149‑1 (MQ=255)
                                                                                                                                               aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 >  2:1251861/1‑125 (MQ=38)
                                                                                                                                               aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                                 <  1:1251861/125‑1 (MQ=38)
                                                                                                                                               aCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt        >  1:2485145/1‑150 (MQ=37)
                                                                                                                                                   ccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCAACACCGCTGTTAAGTGCTCCAACCCGTAGGTTGGATAAGGCGGTTACGCCGCATACGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa    >  1:2274972/1‑150 (MQ=255)
                                                                                                                                                    cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATAAGGCCAACACCCCAGTGAAGTGCTCCACCCCCTAGGTTGGTTAAGGCGGTTACGCCGCAACCCACAACCAACGACCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt   >  1:771992/1‑150 (MQ=11)
                                                                                                                                                     aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc  >  2:2022388/1‑150 (MQ=37)
                                                                                                                                                     aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc  >  2:809978/1‑150 (MQ=37)
                                                                                                                                                     aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAAAAGGCCAACACCCCAGTGAAGAGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCAACCGACAACCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTaac  >  1:2453847/1‑148 (MQ=11)
                                                                                                                                                     |                                                                                                                                                     
ACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATC  >  NC_000913/2304539‑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: