Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,531,282 A→G D276G (GAC→GGC)  pckA → phosphoenolpyruvate carboxykinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,531,2820AG98.6% 645.1 / ‑1.4 212D276G (GAC→GGC) pckAphosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (97/112);  total (98/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

Only 100 of 228 total aligned reads displayed.

TCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTATCCGTC  >  NC_000913/3531192‑3531371
                                                                                          |                                                                                         
tcttcGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCg                                                                                         >  1:4641611/1‑93 (MQ=255)
   tcGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGgcggcg                                                                                      >  1:157856/1‑93 (MQ=255)
     ggCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGGTGgcggcgcg                                                                                    >  1:1240344/1‑91 (MQ=255)
     ggCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCgtg                                                                                    <  1:3750201/93‑1 (MQ=255)
     ggCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCgtg                                                                                    >  1:223596/1‑93 (MQ=255)
      gCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCgtgt                                                                                   <  1:4158830/93‑1 (MQ=255)
      gCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCgtgt                                                                                   >  2:402690/1‑93 (MQ=255)
       ccTTTCCGGCACCGGTAAAACCCCCCTTTCCACCGACCGGAAACGTCGCCGGATTGGCGATGACGAACCCGGCTGGGACGATGGCGGCGTGtt                                                                                  <  2:1471670/93‑1 (MQ=255)
       ccTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGtt                                                                                  <  2:1219017/93‑1 (MQ=255)
       ccTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGtt                                                                                  >  2:2269259/1‑93 (MQ=255)
        cTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGttt                                                                                 <  1:3469614/93‑1 (MQ=255)
        cTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGttt                                                                                 >  1:3400104/1‑93 (MQ=255)
          ttCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTaa                                                                               <  1:671287/93‑1 (MQ=255)
          ttCCGGCACCGATAAAACCCCCCTTTCCGCCGACCCGAAACGGCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTaa                                                                               <  2:1743074/93‑1 (MQ=255)
           tCCGGCACCGGTAAAACCACCCTTGCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAAc                                                                              <  1:2876204/93‑1 (MQ=255)
              ggCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTc                                                                           >  2:2664339/1‑93 (MQ=255)
              ggCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTc                                                                           <  1:3600684/93‑1 (MQ=255)
              ggCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGGTTAACTTc                                                                           >  2:2988014/1‑93 (MQ=255)
               gCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCg                                                                          >  1:3636975/1‑93 (MQ=255)
                cACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTCCGa                                                                         <  2:768718/93‑1 (MQ=255)
                 cccGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGaa                                                                        <  1:899345/92‑1 (MQ=255)
                 aCCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGaa                                                                        <  1:2064513/93‑1 (MQ=255)
                 aCCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGaa                                                                        <  1:1173883/93‑1 (MQ=255)
                  ccGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAg                                                                       >  1:4058996/1‑93 (MQ=255)
                  ccGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAg                                                                       <  1:366084/93‑1 (MQ=255)
                   cGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAgg                                                                      <  2:1100287/93‑1 (MQ=255)
                   cGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAgg                                                                      >  2:2880079/1‑93 (MQ=255)
                    ggTAAAACCACCCTTTCCGCCGCCCCGAACCGTCGCCTGATTGGCGATGGCGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAggc                                                                     <  2:817331/93‑1 (MQ=255)
                    ggTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAggc                                                                     >  2:445313/1‑93 (MQ=255)
                     gTAAAGCCACCCTTTCCACCGACGCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAggcg                                                                    <  2:1591110/93‑1 (MQ=255)
                     gTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAggcg                                                                    >  2:3841997/1‑93 (MQ=255)
                       aaaaCCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAggcggc                                                                  <  1:183268/93‑1 (MQ=255)
                        aaaCCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAAGCGgct                                                                 >  1:4440666/1‑93 (MQ=255)
                          accaccCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGgctgc                                                               >  2:104007/1‑93 (MQ=255)
                           ccaccCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGgctgct                                                              >  1:4564653/1‑93 (MQ=255)
                              cccTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACg                                                           >  2:3749625/1‑93 (MQ=255)
                               ccTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGc                                                          >  1:3929605/1‑93 (MQ=255)
                                   tCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCaaaa                                                      >  1:2649968/1‑93 (MQ=255)
                                    cgaCCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCaaaaa                                                     <  1:2995433/91‑1 (MQ=255)
                                     cACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAAc                                                    <  2:59062/93‑1 (MQ=255)
                                     cACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAAc                                                    >  2:4072407/1‑93 (MQ=255)
                                     cACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAAc                                                    <  1:1274467/93‑1 (MQ=255)
                                     cACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAAc                                                    >  2:3962331/1‑93 (MQ=255)
                                      aCCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACt                                                   >  1:3107999/1‑93 (MQ=255)
                                       ccGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTa                                                  >  2:2570253/1‑93 (MQ=255)
                                       ccGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTa                                                  <  2:2824599/93‑1 (MQ=255)
                                       ccGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTa                                                  <  1:953571/93‑1 (MQ=255)
                                         gACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATc                                                <  2:345889/93‑1 (MQ=255)
                                         gACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATc                                                <  1:3151780/93‑1 (MQ=255)
                                          gcccGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCa                                               <  1:531176/92‑1 (MQ=255)
                                          aCCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCa                                               <  2:2320296/93‑1 (MQ=255)
                                           cccGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCaa                                              >  1:1715192/1‑93 (MQ=255)
                                            ccGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTAtacag                                             >  1:1839018/1‑89 (MQ=255)
                                            ccGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAg                                             >  2:199907/1‑93 (MQ=255)
                                             cGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGc                                            >  1:1531270/1‑93 (MQ=255)
                                              gAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCt                                           >  1:635122/1‑93 (MQ=255)
                                                aaCGTCGCCTGATTGGCGATGACGAACACGCCTGGGACGATGGCGGCGTGTTTAACGTCGAACGCGGCTGCCACGCAAAAACTATCAAGCTGt                                         >  1:1055532/1‑93 (MQ=255)
                                                 aCGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTc                                        <  2:1862336/93‑1 (MQ=255)
                                                   gTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGa                                      >  1:2089783/1‑93 (MQ=255)
                                                   gTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTGAAATTCGAAGGCGGCTGCTACGCCAAAACTATCAAGCTGCCGa                                      >  2:4656767/1‑93 (MQ=255)
                                                    tCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGaa                                     >  2:1529835/1‑93 (MQ=255)
                                                     cGCCTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGaaa                                    <  1:4592112/93‑1 (MQ=255)
                                                        cTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAaagaa                                 >  2:1142699/1‑93 (MQ=255)
                                                        cTGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAaagaa                                 >  1:1001316/1‑93 (MQ=255)
                                                         tGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAaagaag                                <  2:4283393/93‑1 (MQ=255)
                                                         tGATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAaagaag                                >  2:3418256/1‑93 (MQ=255)
                                                          gATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGc                               >  2:1314523/1‑93 (MQ=255)
                                                          gATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGc                               <  2:670684/93‑1 (MQ=255)
                                                          gATTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGc                               >  1:590782/1‑93 (MQ=255)
                                                           aTTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCg                              >  1:806069/1‑93 (MQ=255)
                                                           aTTGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCg                              <  2:1191990/93‑1 (MQ=255)
                                                           aTTGGCGATGACGAACACGGCTGGGACGATGGCGGCGGGTTTTACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCg                              >  1:4684825/1‑93 (MQ=255)
                                                            ttGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCgg                             >  2:1898399/1‑93 (MQ=255)
                                                            ttGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCgg                             >  2:3366907/1‑93 (MQ=255)
                                                             tGGCGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGa                            <  1:1804942/93‑1 (MQ=255)
                                                              tgCGATGACGAACACGGCTGGGACGATGTCAGCGTGGATAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGaa                           <  1:3803126/92‑1 (MQ=255)
                                                                cGATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAAcc                         <  2:2791563/93‑1 (MQ=255)
                                                                 gATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCt                        <  1:2319858/93‑1 (MQ=255)
                                                                 gATGACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCt                        >  2:3811538/1‑93 (MQ=255)
                                                                    gACGAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGaa                     >  2:3923794/1‑93 (MQ=255)
                                                                       gAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATc                  <  2:2449298/93‑1 (MQ=255)
                                                                       gAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATc                  >  2:1611228/1‑93 (MQ=255)
                                                                       gAACACGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATc                  <  1:2180744/93‑1 (MQ=255)
                                                                         acacGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTa                <  2:4416264/93‑1 (MQ=255)
                                                                         acacGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTa                <  1:2712957/93‑1 (MQ=255)
                                                                          cacGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTAc               <  2:3127279/93‑1 (MQ=255)
                                                                            cGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACaa             <  1:2540195/93‑1 (MQ=255)
                                                                            cGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACaa             >  1:2848727/1‑93 (MQ=255)
                                                                            cGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACaa             >  2:370482/1‑93 (MQ=255)
                                                                            cGGCTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACaa             <  1:1407677/93‑1 (MQ=255)
                                                                            cGGCTGGGACGATGGCGGCGTGTTTAAATTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGGACCTGAAATCTACaa             >  2:4640316/1‑93 (MQ=255)
                                                                               cTGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGc          <  2:2296256/93‑1 (MQ=255)
                                                                                tGGGACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCt         >  2:3089310/1‑93 (MQ=255)
                                                                                 gggACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTa        >  2:3260950/1‑93 (MQ=255)
                                                                                 gggACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTa        <  1:1667394/93‑1 (MQ=255)
                                                                                  ggACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTAt       <  1:1876250/93‑1 (MQ=255)
                                                                                   gACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTATc      <  1:1851240/93‑1 (MQ=255)
                                                                                   gACGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTATc      >  2:2520905/1‑93 (MQ=255)
                                                                                    aCGATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTATcc     <  1:712022/93‑1 (MQ=255)
                                                                                      gATGGCGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTATCcgt   >  2:4580799/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCGGAACCTGAAATCTACAACGCTATCCGTC  >  NC_000913/3531192‑3531371

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: