Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,680 T→C 11.3% intergenic (+287/+428) 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,6800TC11.3% 385.0 / 26.3 172intergenic (+287/+428)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (49/96);  new base C (5/13);  total (57/115)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

TAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGC  >  NC_000913/2304531‑2304829
                                                                                                                                                     |                                                                                                                                                     
tAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCt                                                                                                                                                       <  1:1705934/150‑1 (MQ=255)
 aGGTCAGACAAGGCGGTCCCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                                                                                                                                      <  1:2564246/150‑1 (MQ=255)
   gTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                                                                                                                                     >  1:1141083/1‑149 (MQ=255)
   gTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                                                                                                                                     <  2:1141083/149‑1 (MQ=255)
   gTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                                                                                                                                    <  2:133059/150‑1 (MQ=40)
   gTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                                                                                                                                    >  1:2096889/1‑150 (MQ=40)
     cAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTTACGc                                                                                                                                                  <  1:6125788/150‑1 (MQ=38)
     cAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                                                                                                                                  >  1:2358166/1‑150 (MQ=38)
      aGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGcc                                                                                                                                                 >  1:161789/1‑150 (MQ=38)
       gACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGccc                                                                                                                                                >  2:6400535/1‑150 (MQ=37)
       gACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGccc                                                                                                                                                >  1:4947022/1‑150 (MQ=37)
        aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                               <  1:3601099/150‑1 (MQ=35)
        aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                                                               >  2:5495014/1‑150 (MQ=35)
          aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGGGAAGTGCTCCACCCCCTAGGTCGGAGCAGGGGGTGACGCCGCATCCGACATCTAACGCCCagg                                                                                                                                             >  1:2560295/1‑148 (MQ=11)
           tggCTGTCTCGCCGCATCCTACATACAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTTGGTCGTATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGc                                                                                                                                            <  1:1082620/149‑1 (MQ=11)
           aGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGc                                                                                                                                            <  2:1882413/150‑1 (MQ=34)
           aGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGc                                                                                                                                            <  2:3857742/150‑1 (MQ=34)
            ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGcc                                                                                                                                           <  2:1868334/150‑1 (MQ=34)
             gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                                                          >  2:1391808/1‑150 (MQ=34)
               ggTCACGCCTCATCCGACATCCAACGCCCGAGCCGTTTGCCTTATTCGAATCTGGCTCGTCTTATGATGGCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                        <  2:628764/150‑1 (MQ=11)
               ggTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGt                                                                                                                                        <  2:59858/150‑1 (MQ=34)
                gTCTCGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCGTGATTGGACTCTGGCGCGTCTTTTGTTGGCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                       <  2:851420/150‑1 (MQ=11)
                gTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                       <  2:1331504/150‑1 (MQ=34)
                 tCACGCCGCATCCGAGATCCATCGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                                                      <  2:695483/150‑1 (MQ=16)
                  cACGCCGCCTCCGACATCCAACGCCCGAGCCGGTTGCGTGATGCGTCGCTGGCTCGTCTTATCTGTCCTACACAGCTTTGAAGTGCTCGACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                     <  2:6182223/150‑1 (MQ=11)
                   aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                                                                    >  1:2723501/1‑150 (MQ=32)
                      ccgccgccGACCTCCCACGCCCGCGCCGGGTGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACAGCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATACGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                 <  2:4431155/144‑1 (MQ=11)
                      ccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGGTGCCTGa                                                                                                                                 >  1:3414661/1‑150 (MQ=17)
                      ccgcACCCGACATCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCTTCTTATCCGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGa                                                                                                                                 <  1:44758/150‑1 (MQ=255)
                       cgcATCCGACATCCAACGCCCGTGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                <  1:2270308/150‑1 (MQ=17)
                       cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGTAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                <  1:5839623/150‑1 (MQ=17)
                       cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                                                                                                                                    >  1:3725249/1‑130 (MQ=32)
                       cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                <  2:2056307/150‑1 (MQ=32)
                       cgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                <  1:5417416/150‑1 (MQ=32)
                        gcttCCGACATCCTTCGCCCGAGCCGTTTGCCTGATTCTACGCTGGCTCGTCTTATCAGGCCTACACCGCCGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                               <  1:6086342/147‑1 (MQ=11)
                        gcATCCGACATCCAACGCCCGAGCCGTTTGCCTGATGCTACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                               <  1:2203717/150‑1 (MQ=17)
                        gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                               <  1:4867904/150‑1 (MQ=32)
                         cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGGAGGTCGGATAAGGCGGTTACGCCGCATCCGCCAACTAACGCCCAACCCCGTTTCCAGAGGc                                                                                                                              >  1:1583980/1‑150 (MQ=11)
                           tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                     >  1:4579601/1‑141 (MQ=32)
                           tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                                                     <  2:4579601/141‑1 (MQ=32)
                             cGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCTCCACCCCGTATGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                                          <  1:4104882/150‑1 (MQ=17)
                             cGACAGCCCACGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                                          <  2:2412280/150‑1 (MQ=12)
                              gACATCCATCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATGAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                                                         <  2:3194885/150‑1 (MQ=17)
                                 aTCCCACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCGGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                      <  2:2667326/150‑1 (MQ=17)
                                  tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                     <  2:3360119/150‑1 (MQ=17)
                                    cAACTCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTAACAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTAGGATAAGGCGGTTACGCCGCCTCCGACCTCTAACGCCCGAGCCGGTTGCGTGCGGCGGCGCTggggg                                                                                                                   >  2:3749409/1‑147 (MQ=11)
                                    cAACGCCCGAGCCGTTTGCCTGATGCTACGCTGGCGCGACAGAACAGGCCTACACCTACTTGAAGTGCTCCACCCCGTAGGTCGGACAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                   <  2:3048712/150‑1 (MQ=11)
                                    cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                                                   <  1:464631/150‑1 (MQ=32)
                                     aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                  >  2:1297307/1‑150 (MQ=32)
                                     aaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                  >  2:1297288/1‑150 (MQ=32)
                                     aaCGCCCGAGCCGGTTGCCAGATGCGACGCTGGCGCGTCTTATCAGGCCTACAACGCAGTGAAGTGCTCCACCCCGTAGGTCGGTTAAGGCGGTTACGCCACATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                  >  1:4693264/1‑150 (MQ=11)
                                      aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                 <  1:1991980/150‑1 (MQ=32)
                                      aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                 <  1:6446296/150‑1 (MQ=32)
                                      aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                 <  2:2723501/150‑1 (MQ=32)
                                      aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                 <  1:5600761/150‑1 (MQ=32)
                                      aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                 <  2:5453729/150‑1 (MQ=32)
                                      aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                 <  1:1297307/150‑1 (MQ=32)
                                      aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                 <  1:515049/150‑1 (MQ=32)
                                        gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                                                               <  1:2045145/150‑1 (MQ=31)
                                                                                                               cGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                        <  1:3545072/150‑1 (MQ=17)
                                                                                                               cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACAt                                        <  2:6093988/150‑1 (MQ=32)
                                                                                                                gTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATc                                       <  2:3414661/150‑1 (MQ=17)
                                                                                                                   ggTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                    <  1:2998921/150‑1 (MQ=21)
                                                                                                                   ggTCGGAGAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCaa                                    <  2:3584557/150‑1 (MQ=21)
                                                                                                                    gTCGGATATGGCGGTTACGCAGCCTCCGACATCCCTCTGCCCTTCCGTTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc                                   <  1:6172281/150‑1 (MQ=11)
                                                                                                                    gTCGGAGAAGGCGGTGACGCCGCCTCCACCCACCACCGCCCCCGCCGGGTGCCGTCTGCGGCTCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc                                   <  2:4326364/150‑1 (MQ=11)
                                                                                                                    gTCAGACAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc                                   >  1:3011379/1‑150 (MQ=21)
                                                                                                                     tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCGAACg                                  >  1:2956940/1‑150 (MQ=17)
                                                                                                                     tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCCACg                                  >  2:5551204/1‑150 (MQ=18)
                                                                                                                     tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACg                                  >  1:2868127/1‑150 (MQ=34)
                                                                                                                      cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                 <  1:4515675/150‑1 (MQ=35)
                                                                                                                       ggATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTc                                <  2:2235515/150‑1 (MQ=35)
                                                                                                                         aTAAGGCGTTTACGCCGTATCGGAAATCCAACGCCCAAGCCGGTTGCCTGCTGCGACGGTGGCGCGTCTTATCTGGCCTACACCGCTGTTAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCg                              <  1:6090166/150‑1 (MQ=11)
                                                                                                                         aTAAGGCGGTCACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCg                              <  2:4767807/150‑1 (MQ=37)
                                                                                                                          tAAGGCTGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGTTGCGACGCTGGGGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                             <  1:5355954/150‑1 (MQ=255)
                                                                                                                          tAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTTAAGTGCTCCACCCCGTTGGTTGTATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                             <  1:5356237/150‑1 (MQ=255)
                                                                                                                          tAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                             <  2:5245684/150‑1 (MQ=37)
                                                                                                                          tAAGGCGGTTACGCCGCATCCGACATCGAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGa                             <  2:2849767/150‑1 (MQ=37)
                                                                                                                             ggCGGTTACGCCGCATCCGACATCTAACGCCCACGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCTGGCCTACTCCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGcc                          <  1:6615518/150‑1 (MQ=255)
                                                                                                                             ggCGGGCACGCCGCATCCGACATCCAACGCCCCAGCCGGTTGCCCGATGCGACGCTGGCGCGTATTATCAGGCCAACACCGCTGTGAAGTGCGCCACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGcc                          <  1:2097544/150‑1 (MQ=11)
                                                                                                                              gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCg                         <  1:1391808/150‑1 (MQ=37)
                                                                                                                                  ttACGCCACATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTg                     <  2:4693264/150‑1 (MQ=37)
                                                                                                                                   gaCGCCGCACCCGACATCTAACGCCCTTGCCGGTTGGCTGATGCGACGCTGGCTCGTCTTATCAGGCCTACTGCGCTGTGAAGTGGGCCTCCCGGTAGGTAGGATAAGGCGGTAACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                    <  2:4999357/149‑1 (MQ=14)
                                                                                                                                   caCGCCGCAACCGACATATAACTCCCAATCCTGTTGACTGATGTGACGCTGGCTCGTCTTATCAGCCCTAAACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGc                    <  1:618775/149‑1 (MQ=11)
                                                                                                                                    aCGCCGCATCCGACATCCAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTTTGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGcc                   <  2:4948855/150‑1 (MQ=11)
                                                                                                                                      gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCGGGGACGGGCTCCACCCCGGAGGGTGGAGAAGGCGGGTACGCCGCATCCGACATCCCACGCCCGAGCCGGGGGCCTg                 >  2:5981114/1‑150 (MQ=11)
                                                                                                                                      gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                 >  2:6410591/1‑150 (MQ=37)
                                                                                                                                      gccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGACCGAGCCGGTTGCCTg                 >  2:3960666/1‑150 (MQ=21)
                                                                                                                                      gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                 <  2:907066/150‑1 (MQ=21)
                                                                                                                                      gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                 <  2:3091179/150‑1 (MQ=21)
                                                                                                                                      gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                 <  1:4340340/150‑1 (MQ=21)
                                                                                                                                      gccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                 <  1:1920606/150‑1 (MQ=21)
                                                                                                                                         gcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATg              <  1:6254764/150‑1 (MQ=37)
                                                                                                                                            tCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGa           <  1:1294232/150‑1 (MQ=37)
                                                                                                                                              cGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACg         >  1:4373563/1‑150 (MQ=37)
                                                                                                                                                 cATCGACCGCCCCCGCCGGGTGCCTGATGCGCCGCTGGCGCGTCGGTTGAGGGCTACACCGCTGTGAAGTGCTGCACCCCGTAGGTGGGAGACGGCGGTGACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTg      <  2:5357562/150‑1 (MQ=11)
                                                                                                                                                     tAACGCCCAAGCCGGTTGCCTGATGTGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc  >  1:2964786/1‑150 (MQ=37)
                                                                                                                                                     tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgcgc  <  1:6293116/150‑1 (MQ=37)
                                                                                                                                                     tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGAGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGATGCCTGATGCGACGCTGgcgc  >  2:550242/1‑150 (MQ=255)
                                                                                                                                                     tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGGAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTACGAGCCGGGTGCCTGATGCGACGCTGgcgc  >  1:5284269/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGC  >  NC_000913/2304531‑2304829

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: