Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,472 C→A 14.4% intergenic (‑374/+224) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,4720CA14.4% 187.9 / 8.4 89intergenic (‑374/+224)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base A (2/12);  total (38/51)
Fisher's exact test for biased strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGC  >  NC_000913/4153326‑4153621
                                                                                                                                                  |                                                                                                                                                     
gCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                                                                                                                    <  2:608127/150‑1 (MQ=255)
gCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                                                                                                                    <  2:1776891/150‑1 (MQ=255)
 ccccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggtc                                                                                                                                                   >  1:1831479/1‑150 (MQ=255)
  cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggtcg                                                                                                                                                  <  2:735410/150‑1 (MQ=255)
  cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTGGGGAATCGTGTGAGTCACTgtcggtcg                                                                                                                                                  >  1:988867/1‑150 (MQ=255)
   ccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGtcggtcgg                                                                                                                                                 <  2:1294076/150‑1 (MQ=255)
    cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTAACTGTCGGTCGGa                                                                                                                                                >  1:1596899/1‑150 (MQ=255)
    cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggaccga                                                                                                                                                >  1:166019/1‑145 (MQ=255)
    cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                                                                                                                <  1:2370769/150‑1 (MQ=255)
       gTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATaa                                                                                                                                             <  1:1479521/150‑1 (MQ=255)
        tAACACCAAAGGTGTAGGTCGGATAATATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTTTTGTATGAGGCTCGTGCGACTTATGCGTCCTACAGTTGGGCATCGTTTTAGTCAATGTGTGTCGGATAAg                                                                                                                                            <  2:2302830/150‑1 (MQ=14)
         aaCACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGGCGGATAAGa                                                                                                                                           >  1:1237755/1‑150 (MQ=255)
          acacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGAt                                                                                                                                          <  1:2021160/150‑1 (MQ=255)
            acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                                                                        <  1:904365/150‑1 (MQ=255)
            acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                                                                        <  2:2027683/150‑1 (MQ=255)
             cccaacgtgtAGGTCTTATAATATTCGCAATTATCGCATCCGACATTATTGCGGCACTTTATTTTGTCAACATTTCCTGATGCTTCGCTAGCGCCTTATCCTTAATACTGTTGGGCATCTTTTGATTAAGTGTAGGTCGGATAAGATgcg                                                                                                                                       <  2:450078/144‑1 (MQ=2)
             ccATAGGTTTAGGTCGGATAAGATGCTCAAGTATCGGATCCGACATTATTTCGGGACTGGAGTTTGGCAACAGTGCCGGATTCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGTGTCACTGTCGGGCGGATAAGTTgcg                                                                                                                                       <  2:1037860/150‑1 (MQ=11)
             ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAGAGGTGGGCATCGTGTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                                                                       >  1:1977990/1‑150 (MQ=255)
             ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                                                                       <  1:2376112/150‑1 (MQ=255)
              cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcgc                                                                                                                                      <  2:1608041/150‑1 (MQ=255)
                 aGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAg                                                                                                                                   <  2:416352/150‑1 (MQ=39)
                  ggTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGt                                                                                                                                  >  2:1911768/1‑150 (MQ=39)
                   gtgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTa                                                                                                                                 <  2:1441219/150‑1 (MQ=39)
                     gtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATc                                                                                                                               <  1:917975/150‑1 (MQ=39)
                        ggTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCa                                                                                                                            >  1:1806297/1‑150 (MQ=37)
                        ggTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGGCGGATAAGATGCGCAAGTATCGCa                                                                                                                            >  1:1949411/1‑150 (MQ=37)
                         gTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCAt                                                                                                                           <  1:1424563/150‑1 (MQ=37)
                           cGGATAAGATGCGCAAGTATCGCATCCAACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATcc                                                                                                                         <  1:1734702/150‑1 (MQ=37)
                                aaGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACat                                                                                                                    >  1:538095/1‑150 (MQ=37)
                                  gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACatta                                                                                                                  >  1:879900/1‑150 (MQ=37)
                                  gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCATGTGAGTCACTGTCTTTCGGATAAGATGAGCAAGGATCGCATCCGACatta                                                                                                                  >  2:1323349/1‑150 (MQ=11)
                                     gcccAACTATCGCATCCGACATTATTGCGGCACTTGATTTTGTCAACAGTGCCGTATTCGGCTCGTTCGCGTTATCCGTTATACAGTTTGGCATCGTTTGAGTATCTTTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTg                                                                                                               <  2:1560449/147‑1 (MQ=11)
                                      cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGc                                                                                                              <  1:2014519/150‑1 (MQ=37)
                                        cAAGTATCGCATCCGAAATTATTGCGGGACTGGAGGTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTGGGGCAGCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCgg                                                                                                            <  2:1904538/150‑1 (MQ=11)
                                         aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCTGATGCGGCTCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGATTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGc                                                                                                           <  2:1618196/150‑1 (MQ=255)
                                          aGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTTGGAATCATGTAAAAAACTTTCGGTCAGATAAGATGAGCAAGAAACGCATAAGACATTATTgagtca                                                                                                          >  2:516584/1‑145 (MQ=11)
                                          aGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCa                                                                                                          <  1:603890/150‑1 (MQ=37)
                                              tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCTAGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTgg                                                                                                      <  2:598268/150‑1 (MQ=12)
                                                gCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAg                                                                                                    >  1:1890940/1‑150 (MQ=37)
                                                  aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAg                                                                                                                                   <  2:2130371/117‑1 (MQ=35)
                                                  aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAg                                                                                                                                   >  1:2130371/1‑117 (MQ=35)
                                                         attattGCGGCACTGGATTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGTCATACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAAc                                                                                           <  2:1167365/150‑1 (MQ=255)
                                                           tattGCGGCACTGGAGTTTTGCAACAGTGCCGTATGCGGCGCGATAGCCTTATCCGTCATACAGTTGGGCATCGTTTGATTCTCTGTAGGTCGGATAAGATTCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAg                                                                                         <  2:2011775/150‑1 (MQ=11)
                                                               gctgCACTGGAGTTTGGCAACAGTTCCGGATGAGGCGATAGCTCCTTATCCGTCCTAGAGTTGGTCATCTTTTGATTCTGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGcc                                                                                     <  2:1977990/147‑1 (MQ=11)
                                                                cGGCACTGTAGTTTTTCAACCGTGCCGGATGCGGCGCGTGCGCCTTTTCCGTAATACAGTTGGGAATCTTTTCATTCTCTGTCGGTCTGATAATATGCGCAAGTATCGCATCCTACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                                                    <  2:585231/150‑1 (MQ=14)
                                                                cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                                                    <  2:2057120/150‑1 (MQ=255)
                                                                 ggCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAgc                                                                                                                                  >  1:179365/1‑102 (MQ=25)
                                                                     cTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgc                                                                               <  2:307303/150‑1 (MQ=39)
                                                                     cTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGAGGCAATGGAGTTTTGCAACAGTGCCGGATgc                                                                               >  2:338340/1‑150 (MQ=255)
                                                                     cTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTAACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGAAgc                                                                               >  2:232688/1‑150 (MQ=255)
                                                                      tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGTCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcg                                                                              <  2:1806297/150‑1 (MQ=255)
                                                                       ggAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTGTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcgg                                                                             >  1:785277/1‑150 (MQ=255)
                                                                          gTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGgcgc                                                                          >  2:1941665/1‑150 (MQ=39)
                                                                            tttgCAACAGTGCCGTATGCGGCTCGAGCTCCTTATCCGGTCTACAGTTGGGAATCTTTTGATTCAGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                                                        <  2:988867/147‑1 (MQ=11)
                                                                            ttcgCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                                                        <  1:2206562/147‑1 (MQ=255)
                                                                             tGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAg                                                                       <  2:879900/150‑1 (MQ=39)
                                                                                cAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                                                                    >  1:1204681/1‑150 (MQ=39)
                                                                                    aGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCCTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTa                                                                >  2:1501574/1‑150 (MQ=21)
                                                                                    aGTGCCGCATGCCGCGCGAGCGCCTTATCCGCCCTACAGTTGGGCCGCGTCAGAGTCCGTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTa                                                                <  1:989003/150‑1 (MQ=14)
                                                                                       gCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTTAGAATATATTGAGACACTCTAGGTCGGATAAGATGCGAAAGTATCGCATACGACATTATTGCGGAACTAGAGTTTGGAAACAGAGCCGGATGCGGCGCGAGCGCCTTATcc                                                             >  2:2433324/1‑150 (MQ=11)
                                                                                       gCCGCATGCGGCGCGAGCTCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                                             <  1:168278/150‑1 (MQ=255)
                                                                                        ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCg                                                            >  2:1973335/1‑150 (MQ=39)
                                                                                         cGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCCTCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGGGCGCGCGCCTTATCCgg                                                           >  1:1606430/1‑150 (MQ=255)
                                                                                           gATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                         <  1:1229654/150‑1 (MQ=39)
                                                                                                ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGCTAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACgg                                                    >  1:2485474/1‑150 (MQ=255)
                                                                                                   gcgAGCGCCTTATCCGTCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTg                                                 <  1:1321330/150‑1 (MQ=255)
                                                                                                    cgATCGCCTTATCCGGCATACAGTAGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                <  2:2272925/150‑1 (MQ=255)
                                                                                                    cgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                                <  1:1126554/150‑1 (MQ=255)
                                                                                                      aGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGc                                              <  2:539552/150‑1 (MQ=255)
                                                                                                       gcgcCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCa                                             >  1:863993/1‑150 (MQ=255)
                                                                                                        cgcCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCAt                                            >  1:226241/1‑150 (MQ=255)
                                                                                                             tATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGCCATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGAATCCttt                                       >  2:472063/1‑150 (MQ=255)
                                                                                                              aTCCTGCCTACAGTTGGGAATTTTTTGATTCTCTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTg                                      <  2:1237755/150‑1 (MQ=14)
                                                                                                              aTCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGCCATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGCGCGCCTTATCCGGCCTACGGTTGGGCCTCGTTTg                                      >  2:1227928/1‑150 (MQ=255)
                                                                                                              aTCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGAGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCAGGGTTTg                                      >  2:564359/1‑150 (MQ=255)
                                                                                                                ccttccTACAGTTGGGCATCCTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAg                                    <  2:2250056/146‑1 (MQ=255)
                                                                                                                cctgCCTACAGTTGGGCATCTTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATTGTTTGAg                                    <  2:699055/147‑1 (MQ=255)
                                                                                                                ccGTCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAg                                    <  1:2290786/150‑1 (MQ=255)
                                                                                                                    ccTACAGTTGGGCATCGTTTGATTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTTGTAATATTTCCGTATGCTTCGAGAGCGCGTTATACGGACTAAGGTTGGTCATCTTTTgtttgac                                <  2:781155/150‑7 (MQ=14)
                                                                                                                    ccTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCAc                                <  2:1949411/150‑1 (MQ=255)
                                                                                                                            tGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTc                        >  1:1875777/1‑150 (MQ=255)
                                                                                                                            tGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGCAATCGGTTTAGTCACTGTAGGTc                        >  2:2371855/1‑150 (MQ=255)
                                                                                                                             gggCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATGCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCGTCGTTGGAGTGCCTGTAGGTCg                       >  2:624869/1‑150 (MQ=17)
                                                                                                                               gCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGcc                                                                   >  2:2448501/1‑104 (MQ=17)
                                                                                                                                    gTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGa                >  1:563882/1‑150 (MQ=255)
                                                                                                                                    cactGACTCCCTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGa                <  2:1831479/147‑1 (MQ=255)
                                                                                                                                      ttGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGAAACTTTAAAATAACTGAAGATCGGATAAGATg              >  2:941534/1‑150 (MQ=14)
                                                                                                                                       tGAGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgc             <  1:151223/150‑1 (MQ=255)
                                                                                                                                       tGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAg                                                                                                    >  2:1666954/1‑63 (MQ=255)
                                                                                                                                       tGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATcc                                                             >  2:2318141/1‑102 (MQ=17)
                                                                                                                                        gAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg            >  2:84917/1‑150 (MQ=255)
                                                                                                                                        cattccgtGTCGGGCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg            <  1:2455301/143‑1 (MQ=255)
                                                                                                                                          gTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCCGCGCGAGAGCCTTATCCGGCCTACGTTTGTGCAGAGTTTGAATCACTGTAGGTCGGATAAGAtcagca          >  2:537912/1‑145 (MQ=14)
                                                                                                                                            cACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACTGGTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAg        <  1:1523260/150‑1 (MQ=255)
                                                                                                                                            agcTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGgatgcgatgcttgcg        <  1:192340/148‑6 (MQ=11)
                                                                                                                                             aCTGTCGGTCGGATAAGATGCGCAAGTATCGCAGCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgc       <  1:338340/150‑1 (MQ=255)
                                                                                                                                                 tcggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcg   <  1:1719395/150‑1 (MQ=38)
                                                                                                                                                 gcggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcg   <  1:1911768/149‑1 (MQ=38)
                                                                                                                                                  cggtcggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatcgc  >  1:474982/1‑150 (MQ=38)
                                                                                                                                                  |                                                                                                                                                     
GCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGC  >  NC_000913/4153326‑4153621

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: