Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,447 A→G 7.4% intergenic (‑349/+249) 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,4470AG 7.4% 335.3 / 7.7 119intergenic (‑349/+249)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (47/61);  new base G (7/2);  total (56/63)
Fisher's exact test for biased strand distribution p-value = 7.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

TAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAG  >  NC_000913/4153301‑4153596
                                                                                                                                                  |                                                                                                                                                     
tAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGtt                                                                                                                                                    >  1:539022/1‑150 (MQ=255)
  aCTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                                                  <  2:2524292/150‑1 (MQ=255)
    ttATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGc                                                                                                                                                >  2:2094088/1‑150 (MQ=255)
       tGGAAATGTTAAAAAATCGCCCCAAGTACCACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATc                                                                                                                                             <  2:256227/150‑1 (MQ=255)
          aaaTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGtt                                                                                                                                          <  1:1365621/150‑1 (MQ=255)
             tGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGa                                                                                                                                       <  1:650140/150‑1 (MQ=255)
             tGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGa                                                                                                                                       <  1:1759454/150‑1 (MQ=255)
               ttAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGAGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGAAACGTATGAGt                                                                                                                                     >  2:2211127/1‑150 (MQ=17)
                         gCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGGATCGCTTCCGACATGATTGCGGCACTGGAGGGTGGCAACAGTGCCGGCTGCGGCGCGAGCGCGGTGACCGGCCGACGGTTGGGCTTCGGTGGTGTCTCTgtcggg                                                                                                                           >  1:961320/1‑149 (MQ=17)
                          ccccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTGACTgtaggtc                                                                                                                          >  2:1941186/1‑150 (MQ=255)
                           cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGCATCGTAAGAGTCACTgtaggtcg                                                                                                                         >  1:1124453/1‑150 (MQ=14)
                           cccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggtcg                                                                                                                         >  1:707404/1‑150 (MQ=255)
                            ccAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGtcggtcgg                                                                                                                        >  2:1372291/1‑150 (MQ=255)
                             cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGtaggtcggc                                                                                                                       >  2:1312351/1‑149 (MQ=255)
                             cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                                                                                       >  2:2284580/1‑150 (MQ=255)
                                 tAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAg                                                                                                                   >  2:657052/1‑150 (MQ=255)
                                  aaCACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGa                                                                                                                  <  2:1112178/150‑1 (MQ=255)
                                    cacCAAAGGTGTAGGTCGGATTAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATg                                                                                                                <  2:1316975/150‑1 (MQ=255)
                                     acCACAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                                               <  2:2252112/150‑1 (MQ=255)
                                     acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                                                               <  2:159030/150‑1 (MQ=255)
                                      ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg                                                                                                              <  2:621842/150‑1 (MQ=17)
                                      ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                                              <  2:735591/150‑1 (MQ=255)
                                      ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                                              <  2:136564/150‑1 (MQ=255)
                                      ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                                              >  2:1390638/1‑150 (MQ=255)
                                         aaGCTGTAGGACGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTATGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGATTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                                                           <  2:185316/150‑1 (MQ=255)
                                            gtgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACac                                                                                                                                                      <  2:1495023/104‑2 (MQ=25)
                                            gtgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAg                                                                                                                                                      >  1:1495023/1‑104 (MQ=35)
                                              gtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATc                                                                                                      <  2:917399/150‑1 (MQ=39)
                                                   tCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                                                                                 <  2:2585634/150‑1 (MQ=37)
                                                   tCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                                                                                 <  2:1248146/150‑1 (MQ=37)
                                                   tCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                                                                                 <  2:2456377/150‑1 (MQ=37)
                                                      gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGa                                                                                              <  2:1591740/150‑1 (MQ=37)
                                                      gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGa                                                                                              <  2:946251/150‑1 (MQ=37)
                                                           gATGCGCAAGTATCGCATCCGACATTATTGTGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACatta                                                                                         <  1:107118/150‑1 (MQ=37)
                                                                 cAAGTATCGCATCCGGCATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCgg                                                                                   <  1:772283/150‑1 (MQ=37)
                                                                  aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGc                                                                                  <  1:2545505/150‑1 (MQ=37)
                                                                  aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCATTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGc                                                                                  <  1:1027784/150‑1 (MQ=37)
                                                                  aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCATTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGc                                                                                  <  1:1027803/150‑1 (MQ=37)
                                                                  aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCATTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGc                                                                                  <  1:2314521/150‑1 (MQ=37)
                                                                        cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                                            <  2:1304573/150‑1 (MQ=37)
                                                                        cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                                            <  2:944100/150‑1 (MQ=37)
                                                                        cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                                            <  1:1323427/150‑1 (MQ=37)
                                                                        cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGATTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                                            <  2:944083/150‑1 (MQ=37)
                                                                          aaTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGAATCGTTTGAGACACTGTCTGTCGTATAAGATGAGCACGAATCGCATCCGACATTATTGCGGCACTGGAGt                                                                          >  2:2245281/2‑150 (MQ=11)
                                                                            tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCACGAATCGCATCCGACATTATTGCGGCACTGGAGttt                                                                        >  1:1967118/1‑150 (MQ=21)
                                                                              cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                      <  1:62008/150‑1 (MQ=37)
                                                                               gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGTTGCGCAAGTATCGCTTCCGACATTATTGCGGCACTGGAGTTTGGc                                                                     >  1:1031771/1‑150 (MQ=255)
                                                                               gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                     >  2:43368/1‑150 (MQ=37)
                                                                               gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGAATCGCATCCGACATTATTGCGGCACTGGAGTTTGGc                                                                     >  2:43373/1‑150 (MQ=21)
                                                                               gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTAACTGTCGGTCGGATAAGATGCGCAAGAATCGCATCCGACATTATTGCGGCACAGGAGTTTGGc                                                                     >  2:1612296/1‑150 (MQ=255)
                                                                                 cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                   <  1:1877628/150‑1 (MQ=37)
                                                                                 cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                   <  1:601664/150‑1 (MQ=37)
                                                                                 cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                   <  1:2120686/150‑1 (MQ=37)
                                                                                   ttattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACa                                                                 <  1:1328204/150‑1 (MQ=37)
                                                                                     attGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGt                                                               <  1:594954/150‑1 (MQ=38)
                                                                                         cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGCGCGCCTGATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                           <  2:2458054/150‑1 (MQ=255)
                                                                                         cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                           >  1:127350/1‑150 (MQ=39)
                                                                                           gctcTGGAGTTTGGCGACCGTGCCGTAAGGGGCTCGTGCGCGTAATGCGGCCTACAGGGGGGCATCTTCTGAGTCCGTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGa                                                         <  2:2501768/147‑1 (MQ=14)
                                                                                               tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACCGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCa                                                                                 >  1:1465606/1‑122 (MQ=25)
                                                                                               tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcg                                                     <  1:147504/150‑1 (MQ=39)
                                                                                               tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATgcg                                                     >  1:2215925/1‑150 (MQ=37)
                                                                                               tGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCa                                                                                 >  1:1430230/1‑122 (MQ=17)
                                                                                                    tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGGTTGGCAACAGTGCCGGATGCGGcgcg                                                >  1:2444520/1‑150 (MQ=255)
                                                                                                    tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACAGGAGTTTGGCAACAGTGCCGGATGCGGcgcg                                                >  1:2497414/1‑150 (MQ=255)
                                                                                                         cAACAGTTCCGGATTCGGCGCGAGCGCCTTATCCGTCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                                           <  1:167905/150‑1 (MQ=255)
                                                                                                         cAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgcgc                                           <  2:636077/150‑1 (MQ=39)
                                                                                                          caCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                                          <  2:384211/149‑1 (MQ=39)
                                                                                                          aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                                          >  2:312107/1‑150 (MQ=39)
                                                                                                             aGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGGGTGAGTCACTGTCGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTa                                       >  2:2161341/1‑150 (MQ=255)
                                                                                                               tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                                     >  2:1353267/1‑150 (MQ=21)
                                                                                                                 ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCg                                   >  2:234262/1‑150 (MQ=37)
                                                                                                                   ggataaggcgCTCGCGCCCCATCCGGCCTACAGTTGGGCATCGTTTGAGTCACGGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACAGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                 >  2:2070394/7‑150 (MQ=14)
                                                                                                                   ggataaggcgCGCGCGCCGCAACCCGCCTACAGTTGGGCCGCGGTTGAGGCCCTGTCGGTCGGTTAAGATTCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACACTGCCGGCTGCGGGGCGGGCGCCGTATCCggg                                 >  2:1746272/7‑149 (MQ=2)
                                                                                                                   ggATGCGGCGCGAGCGCCTTATCCGGCCTACCGTTGGGCAACGTGTGAGTCGCTGTCGGTCGGCGAAGATGCGCTAGGATCGCATCAGTCATTAGAGAGGGTCGGGAGATTGGCCACTGTGACGGAGGCGGCGCGAGCGCCTTGTACGGc                                 >  2:1296313/1‑150 (MQ=2)
                                                                                                                       gcggcgCGAGCGCCTTATCCGTCCTACAGGTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAc                             <  1:671536/150‑1 (MQ=255)
                                                                                                                         ggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGATTCGCATCCGACATTAGTGCGGCACTGGATTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACgt                           >  2:2550584/1‑149 (MQ=11)
                                                                                                                         ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACgg                           >  2:2550625/1‑150 (MQ=40)
                                                                                                                                   ccTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCg                 >  1:1748363/1‑150 (MQ=255)
                                                                                                                                   ccTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCAACg                 >  1:2164776/1‑150 (MQ=255)
                                                                                                                                    cTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGt                >  1:2274775/1‑150 (MQ=255)
                                                                                                                                    cTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGt                >  1:1999625/1‑150 (MQ=255)
                                                                                                                                     ttATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCggt               >  1:2039610/1‑148 (MQ=255)
                                                                                                                                     ttATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGtt               <  2:2283454/150‑1 (MQ=255)
                                                                                                                                     ttATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGtt               <  1:1832546/150‑1 (MQ=255)
                                                                                                                                     ttATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAGCAGATGCGCAAGGATCGCATCCGACATTCTTGCGGCACTGGAGGTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGGCTACGGTTGGGCGTCggg               >  1:822000/1‑148 (MQ=14)
                                                                                                                                     ttATCCGGCCTACAGTTGGGCATCGTTTCAGTCACTGTCGGTCGGATAAGATGCGCAAGAAACGCATACGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGTGCAACGtt               >  1:1628336/1‑150 (MQ=14)
                                                                                                                                     ttATCCGGCCTACAGTTGGGCATCGATTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGtt               <  1:535738/150‑1 (MQ=255)
                                                                                                                                     atATCCGGCCTACAGTTGGGCATAGTTAGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGtt               <  1:2427888/149‑1 (MQ=255)
                                                                                                                                       aTCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTg             >  1:1002970/1‑150 (MQ=255)
                                                                                                                                       aTCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTg             >  1:2401269/1‑150 (MQ=255)
                                                                                                                                       aTCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTg             >  1:281437/1‑150 (MQ=255)
                                                                                                                                        tCCGGCCTACAGTTGCGCATCCTCTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGTCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa            <  1:1181798/150‑1 (MQ=14)
                                                                                                                                        tCCGGCCTACAGGTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa            <  2:247479/150‑1 (MQ=255)
                                                                                                                                          cGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGt          >  2:115054/1‑150 (MQ=255)
                                                                                                                                           ggCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTc         >  2:1713576/1‑150 (MQ=255)
                                                                                                                                           ggCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTc         >  1:2028574/1‑150 (MQ=255)
                                                                                                                                                aCAGTTGGGCATCGTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGt    <  2:1258653/150‑1 (MQ=255)
                                                                                                                                                  aGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAg  <  1:578846/150‑1 (MQ=255)
                                                                                                                                                  aGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTGTGAGTCACTGTAg  >  2:1482190/1‑150 (MQ=255)
                                                                                                                                                  aGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATCATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAg  <  1:1124989/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAG  >  NC_000913/4153301‑4153596

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: