Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,472 C→A 11.5% 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,4720CA11.5% 194.8 / 4.3 87intergenic (‑374/+224)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (28/47);  new base A (5/7);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATC  >  NC_000913/4153326‑4153619
                                                                                                                                                  |                                                                                                                                                   
gCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                                                                                                                  <  1:790740/150‑1 (MQ=255)
gCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                                                                                                                  >  2:2015133/1‑150 (MQ=255)
  cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTTGGAATCGTATTAGTCACTgtaggtcg                                                                                                                                                >  2:1048933/1‑150 (MQ=14)
  cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcgttcg                                                                                                                                                <  1:54788/150‑1 (MQ=255)
  cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCAGCGTTTGAGTCACTgtcggtcg                                                                                                                                                >  1:1785587/1‑150 (MQ=255)
   ccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGtcggtcgg                                                                                                                                               <  2:869773/150‑1 (MQ=255)
   ccAAGTAACACCAAAGGTGTAGGGCGGATAAGATGCGCACGTATCGCATCCGACATTATTGCGGCACTGGCGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGtcggtcgg                                                                                                                                               <  1:1210077/150‑1 (MQ=255)
    cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                                                                                                              <  1:548884/150‑1 (MQ=255)
      aGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTAACGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATa                                                                                                                                            <  1:104640/150‑1 (MQ=255)
        tAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                                                                          >  1:1724093/1‑150 (MQ=255)
           cacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATg                                                                                                                                       >  2:1788757/1‑150 (MQ=255)
            acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACGGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                                                                      <  1:587183/150‑1 (MQ=255)
            acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgg                                                                                                                                      >  1:1933327/1‑149 (MQ=255)
            acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                                                                      <  1:136363/150‑1 (MQ=255)
               aaaGGTGTAGGTCGGATACGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCa                                                                                                                                   <  1:857877/150‑1 (MQ=255)
               aaaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCa                                                                                                                                   <  2:1721799/150‑1 (MQ=40)
                 tggTGTATGTCGGATAAGATGCGCAAGTCTCTCATCCTACAGTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCTGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAg                                                                                                                                 <  2:156195/149‑1 (MQ=11)
                 agttgtAGGTCTGATAATATGCTCAAGTATCGCATACGACATTATTGCGGCACTGGAGTTGGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAg                                                                                                                                 <  1:1762450/147‑1 (MQ=11)
                       aGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGc                                                                                                                           >  2:1508022/1‑150 (MQ=17)
                           cGGATACGATGCGCAATTATCGAATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGATCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATcc                                                                                                                       <  2:164775/150‑1 (MQ=255)
                           cGGATAAGATGCGAAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCACATcc                                                                                                                       <  2:1480060/150‑1 (MQ=255)
                            gtaTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCg                                                                                                                      <  1:740865/148‑1 (MQ=37)
                              aTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGAc                                                                                                                    <  1:1063767/150‑1 (MQ=37)
                                      cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGc                                                                                                            <  2:1583571/150‑1 (MQ=37)
                                       gcAAGTATCGCATCAGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCg                                                                                                           <  1:725940/150‑1 (MQ=37)
                                        cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCgg                                                                                                          <  1:79674/150‑1 (MQ=37)
                                               cGCATCCGACATTATTGATTCACTGGAGTTTGGCAACAGTGCCGTATGCGTCTCTAGCGACTTATCCTGCCTACAGTTGTGCATAGTTTGATTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                                                                   <  2:1933327/150‑1 (MQ=11)
                                                ccATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAg                                                                                                  <  1:1065837/149‑1 (MQ=37)
                                                    cagtcATTATTGCGGCACTGTATTTTGGCAACAGTGCCGGATTCGGCGCGAGCGCCTTATCCGGACTACTGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTg                                                                                              <  1:763499/146‑1 (MQ=11)
                                                     ggACATTATTGCGGCACTGGAGTTTGTCAACTGTGCCGGATTCGGCGCGAGCGCCTTATCCTGCATACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                                             <  2:1700012/149‑1 (MQ=11)
                                                     cGACATTATTGCGGCACTTGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCACGGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                                            <  2:1154241/150‑1 (MQ=12)
                                                     cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                                             <  1:1641736/150‑1 (MQ=37)
                                                      gACATTATTGCGGCACTGTATTTTGTCAACTGTGTCGTATGCGGATCGAGCGCCTTATCCGTCCTACAGGTGGGCAGCGTTTGATTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                                            <  2:982176/150‑1 (MQ=11)
                                                         attattGCGGCACTGGAGTTTGGCAACAGTGCATTATGCGGCTCGATCGCCTTATCCTGCATACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAAc                                                                                         <  2:1390753/150‑1 (MQ=11)
                                                              tGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATg                                                                                                                                       >  1:1883683/1‑99 (MQ=37)
                                                                cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATGCGGCCTACAGTTGGGCATCGTGTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                                                  <  1:1061547/150‑1 (MQ=255)
                                                                cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACATTTTTGAATCGTTTGAGTCACTGTAGGTCGGATAAGATGAGAAAGTATCGCATCCGCCATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                                                  >  2:1370594/1‑150 (MQ=11)
                                                                cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                                                  <  2:1372500/150‑1 (MQ=39)
                                                                cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                                                  <  1:1629818/150‑1 (MQ=39)
                                                                   cACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGAATCGCATCCGCCATTATTGCGGCACTGTAGTTTGGCAACAGAGCCGGAt                                                                               >  2:1781746/1‑150 (MQ=255)
                                                                    aCTTGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATg                                                                              <  1:1099331/150‑1 (MQ=255)
                                                                    aCTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATg                                                                              <  1:1788757/150‑1 (MQ=39)
                                                                       ggAGTTTGGCAACAGTGACGGATGCGGCGAGAGCGCGTTATCCGGACTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcgg                                                                           <  2:1570452/150‑1 (MQ=255)
                                                                          gTTTGGCAACAGTGCCGGATGCGTCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTATAGTCACTGTAGGTCTGATAAGATTCGCAATTATAGAATCCGACATTATTGCTGCACTGTATTTTGGAAACAGTGCCGGATGCGgcgc                                                                        >  2:1867744/1‑150 (MQ=11)
                                                                          gTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGgcgc                                                                        >  1:1043397/1‑150 (MQ=39)
                                                                          gTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATAAGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGgcgc                                                                        >  2:294026/1‑150 (MQ=255)
                                                                              ggCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgc                                                                    <  2:465731/150‑1 (MQ=255)
                                                                                cAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCCTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAtcgc                                                                  >  1:1442087/1‑150 (MQ=255)
                                                                                cAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                                                                  >  2:7145/1‑150 (MQ=39)
                                                                                 aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                                                                 <  1:1536450/150‑1 (MQ=39)
                                                                                  accgTGCCGTATCCGGCTCGAGCTCCTTATCCTTCCTAAAGTTGGGCATATTTTGCGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCt                                                                <  1:1335318/147‑1 (MQ=14)
                                                                                  aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCt                                                                <  1:1881113/150‑1 (MQ=39)
                                                                                    aGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTa                                                              <  1:2015133/150‑1 (MQ=39)
                                                                                     gTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                             >  1:1675075/1‑150 (MQ=39)
                                                                                       gtcGGATGGGGGGCGAGCGGCTTATCCGTCATACTGTTGGGCAGCGTTTGAGTCTGTGTAGGTAGTATAAGATGAGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                                           <  2:751664/148‑1 (MQ=11)
                                                                                        tcgtaTGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCg                                                          <  2:1274974/146‑1 (MQ=255)
                                                                                        ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGTGCGAGCGCCTTATCCg                                                          >  2:6573/1‑150 (MQ=255)
                                                                                          ggATGCGGCTCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                        <  2:1713133/150‑1 (MQ=255)
                                                                                           gATGCGGCGCGAGCGCCTTATCCGGCCTACAGGTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                       <  1:1955424/150‑1 (MQ=255)
                                                                                                 gcgcGAGCGCCTTATCCGGCCTACAGTTTGGCATCGTTTTATAAACTGTATTTCGGATAATATGCGCAAGAATGGCATCCTACATTATTGCGTCACTGGAGTTTTGCAACAGAGCCGGATGCGGCGCGAGCGACTTATCCGGCCTACGGt                                                 >  2:1884187/1‑150 (MQ=11)
                                                                                                 gcgcGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTgg                                                                                                    >  2:470154/1‑99 (MQ=18)
                                                                                                   gcgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTg                                               <  2:525550/150‑1 (MQ=255)
                                                                                                    cgATCGCCTTATCCGTCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                              <  2:1637924/150‑1 (MQ=255)
                                                                                                    cgAGCGCGTTATCCGGCCTACAGTTGGGAATAGATTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                              <  2:958771/150‑1 (MQ=255)
                                                                                                    cgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATAGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                                              >  1:489707/1‑150 (MQ=255)
                                                                                                     gAGCGCCTTATCCGGCCTACAGGTGGGCAGCGTTTGACTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTggg                                             <  2:1000970/150‑1 (MQ=255)
                                                                                                      aGCGCCTTATCCTGCCTACAGTTGTGCAGCGATTGAGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTTCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGc                                            <  1:865094/150‑1 (MQ=14)
                                                                                                        cgcCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGACTACGGTTGGTCAt                                          >  1:990702/1‑150 (MQ=255)
                                                                                                          cgttATGCGGCCTACAGGTGGGCATCGTTTGAGTCAGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCg                                        <  2:1963534/148‑1 (MQ=255)
                                                                                                          ccTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGGCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCg                                        <  2:584839/150‑1 (MQ=255)
                                                                                                             tATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                     <  2:1129051/150‑1 (MQ=255)
                                                                                                              aTCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGtttt                                    >  2:1471032/1‑149 (MQ=255)
                                                                                                               tCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGCAACGTTTGa                                   >  2:1006000/1‑150 (MQ=255)
                                                                                                                ccGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTATGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGCGCAACGTTagag                                  >  1:934735/1‑150 (MQ=14)
                                                                                                                ccGGCATACAGTTGGGCATCGTTTGCGTCTCTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAg                                  <  1:1992590/150‑1 (MQ=255)
                                                                                                                    aatACAGGTGGGCAGATTTTGAGTATCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCAc                              <  2:393627/148‑1 (MQ=14)
                                                                                                                     cTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACt                             >  2:18370/1‑150 (MQ=255)
                                                                                                                     ataaaGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACt                             <  2:994511/146‑1 (MQ=255)
                                                                                                                      tACAGTTGGTCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTg                            <  1:449291/150‑1 (MQ=255)
                                                                                                                      tACAGTTGGGCAGCGTTTGCGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTg                            <  1:24929/150‑1 (MQ=255)
                                                                                                                        cAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGAGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGAATCGTATGAGTGACTGTa                          >  1:553613/1‑150 (MQ=255)
                                                                                                                        cAGTTGGGCAACGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATACGGCCTACGGTTGGGCATCGTTTGAGACACTgaa                          >  1:1813885/1‑148 (MQ=255)
                                                                                                                            tGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTc                      <  2:1912978/150‑1 (MQ=255)
                                                                                                                            tGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGGCTCGTTTAATTAACTGAAGTTc                      >  1:379110/1‑150 (MQ=14)
                                                                                                                               gCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGCAACGTTTGATTCACTGTAGTTCCGAt                  >  2:1234112/1‑150 (MQ=17)
                                                                                                                                   cGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTGGGAATCGTTTGAGTCACTGTAGTTCGGATaat               >  2:792140/1‑149 (MQ=255)
                                                                                                                                   cGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTATTAGTCACTGTAGTTCGGATAAg               >  2:372166/1‑150 (MQ=255)
                                                                                                                                    gTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGa              >  1:1412671/1‑150 (MQ=255)
                                                                                                                                    atttGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTGGCATCGTTTGAGTCACTGTAGGTCGGATAAGa              >  1:547067/2‑150 (MQ=255)
                                                                                                                                      ttgtttCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCAGCGTTTGAGTCACTGTAGGTCGGATAAGATg            <  2:890234/145‑1 (MQ=255)
                                                                                                                                        gAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGATACGTTAGTCTCAGAGTAGTACCGATAAGATgag          >  1:1982545/1‑148 (MQ=14)
                                                                                                                                         ctgcccTGTCGGTCGGATAAGATGCGCAAGTATCGCATCGGACATTATTGCGGCAGTGGAGGTTGGGTACTGTGCCGGATGCGGCGCGAGCGGCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCTCTGTAGGTCGGATAAGATgcgc         <  2:886758/145‑1 (MQ=11)
                                                                                                                                         aGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTAGGGAATCGTTTGAGTCACTGTAGGTCGGTTAAGATccgc         >  1:1262257/1‑150 (MQ=255)
                                                                                                                                               tGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGAGCATCGTATGATTCACTGTAGGTCGGATAAGATGCGCAAgcat   >  1:174536/1‑150 (MQ=255)
                                                                                                                                                gtcggtcgGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGTAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc  <  1:551730/150‑1 (MQ=255)
                                                                                                                                                gtcggtcgGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCACgcctc  >  2:966752/1‑150 (MQ=255)
                                                                                                                                                gtcggtcgGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcatc  <  1:1179852/150‑1 (MQ=39)
                                                                                                                                                  |                                                                                                                                                   
GCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATC  >  NC_000913/4153326‑4153619

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: