Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 862,272 T→C T164A (ACT→GCT)  ybiY ← putative pyruvate formate‑lyase activating enzyme 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913862,2720TC99.6% 795.8 / ‑5.0 254T164A (ACT→GCT) ybiYputative pyruvate formate‑lyase activating enzyme 2
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (148/105);  total (148/106)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

Only 100 of 263 total aligned reads displayed.

CGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGC  >  NC_000913/862182‑862361
                                                                                          |                                                                                         
 gtgtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATg                                                                                        >  1:7279368/1‑93 (MQ=255)
 gtgtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATg                                                                                        >  1:3977572/1‑93 (MQ=255)
  tgtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGa                                                                                       >  2:1357080/1‑93 (MQ=255)
   gtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGaa                                                                                      >  1:390643/1‑93 (MQ=255)
   gtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGaa                                                                                      <  1:2354157/93‑1 (MQ=255)
    ttttAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAt                                                                                     >  2:7325803/1‑93 (MQ=255)
    ttttAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAt                                                                                     <  2:6352991/93‑1 (MQ=255)
     tttAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtg                                                                                    >  2:108713/1‑93 (MQ=255)
     tttAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtg                                                                                    >  2:3596522/1‑93 (MQ=255)
      ttAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgc                                                                                   >  2:1915302/1‑93 (MQ=255)
       tAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcc                                                                                  <  2:777527/93‑1 (MQ=255)
       tAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcc                                                                                  <  2:6191660/93‑1 (MQ=255)
        aaaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcct                                                                                 >  1:1051062/1‑93 (MQ=255)
        aaaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcct                                                                                 <  1:1318997/93‑1 (MQ=255)
        aaaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcct                                                                                 >  2:6370826/1‑93 (MQ=255)
         aaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctg                                                                                >  1:4111332/1‑93 (MQ=255)
         aaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctc                                                                                >  2:4815478/1‑92 (MQ=255)
          aTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgc                                                                               <  2:4869500/93‑1 (MQ=255)
          aTCAGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgc                                                                               >  2:2742671/1‑93 (MQ=255)
           tCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgcc                                                                              <  2:3104490/93‑1 (MQ=255)
           tCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgcc                                                                              >  1:6094313/1‑93 (MQ=255)
           acGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgcc                                                                              <  2:5435723/92‑1 (MQ=255)
           acGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgcc                                                                              <  1:1738657/92‑1 (MQ=255)
            cGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATgcctgcct                                                                             >  2:5641964/1‑93 (MQ=255)
             ggCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTc                                                                            >  1:4550944/1‑93 (MQ=255)
              gCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCg                                                                           <  2:5287692/93‑1 (MQ=255)
               cAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCGCATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGt                                                                          <  1:4954900/93‑1 (MQ=255)
                aaGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtg                                                                         >  2:6585847/1‑93 (MQ=255)
                 aGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtgg                                                                        >  2:5752437/1‑93 (MQ=255)
                 aGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtgg                                                                        >  1:6275055/1‑93 (MQ=255)
                  gAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtggg                                                                       >  2:627322/1‑92 (MQ=255)
                  gAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtggc                                                                       <  1:3848813/93‑1 (MQ=255)
                        gATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCtt                                                                 >  2:3454646/1‑93 (MQ=255)
                         aTCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTg                                                                >  1:3346038/1‑93 (MQ=255)
                            gATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGccg                                                             >  2:4411637/1‑91 (MQ=255)
                            gATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAg                                                             >  2:4654768/1‑93 (MQ=255)
                            gATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAg                                                             <  1:5386380/93‑1 (MQ=255)
                             atatAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGt                                                            <  2:1634881/93‑1 (MQ=255)
                             atatAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGt                                                            >  2:7905713/1‑93 (MQ=255)
                              tataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTa                                                           <  2:7368923/93‑1 (MQ=255)
                              tataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTa                                                           >  2:892304/1‑93 (MQ=255)
                              tataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTa                                                           <  1:2882846/93‑1 (MQ=255)
                              tataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTa                                                           >  1:4827875/1‑93 (MQ=255)
                               ataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAg                                                          >  1:3314814/1‑93 (MQ=255)
                               ataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAg                                                          <  2:3733746/93‑1 (MQ=255)
                                   ggCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgcc                                                      >  2:5340015/1‑93 (MQ=255)
                                   ggCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgcc                                                      <  1:6759939/93‑1 (MQ=255)
                                    gCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgcca                                                     >  1:423660/1‑93 (MQ=255)
                                    gCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgcca                                                     >  1:783450/1‑93 (MQ=255)
                                    gCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgcca                                                     >  1:1659803/1‑93 (MQ=255)
                                    gCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgcca                                                     <  2:7168277/93‑1 (MQ=255)
                                          aaGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgcca                                               <  1:820429/93‑1 (MQ=255)
                                           aGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgccat                                              <  2:3312702/93‑1 (MQ=255)
                                           aGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgccat                                              >  2:1780301/1‑93 (MQ=255)
                                           aGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgccat                                              <  1:7643136/93‑1 (MQ=255)
                                            ggggCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCAtt                                             <  2:6986892/93‑1 (MQ=255)
                                             gggCGATATATTGCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCAttt                                            <  2:3799053/93‑1 (MQ=255)
                                               gCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTcc                                          >  1:1865215/1‑93 (MQ=255)
                                                 gATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCgg                                        <  2:6791521/93‑1 (MQ=255)
                                                 gATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCgg                                        >  1:5659506/1‑93 (MQ=255)
                                                 gATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCgg                                        >  1:440469/1‑93 (MQ=255)
                                                   tatatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTGGCGCCATCGCCATTTCCGGCt                                      >  1:2605720/1‑93 (MQ=255)
                                                   tatatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCt                                      >  2:6838314/1‑93 (MQ=255)
                                                    atatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                     <  2:1408364/93‑1 (MQ=255)
                                                    atatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                     >  1:6788912/1‑93 (MQ=255)
                                                     tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGc                                    >  2:2269275/1‑93 (MQ=255)
                                                       tttCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCAt                                  >  2:3647510/1‑93 (MQ=255)
                                                          ccACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATaaa                               >  1:3477060/1‑93 (MQ=255)
                                                          ccACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATaaa                               <  1:6610403/93‑1 (MQ=255)
                                                          ccACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATaaa                               <  2:4718382/93‑1 (MQ=255)
                                                          ccACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATaaa                               <  2:4670831/93‑1 (MQ=255)
                                                           cACGGCACATGCAGACAGGTTTCTACCGCAGCATGATTGCCTGACTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATACAg                              <  2:7037344/93‑1 (MQ=255)
                                                           cACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAg                              <  1:2156739/93‑1 (MQ=255)
                                                           cACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAg                              <  2:1878691/93‑1 (MQ=255)
                                                              ggCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGc                           <  1:2111659/93‑1 (MQ=255)
                                                               gCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCCGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCt                          >  1:6215454/1‑93 (MQ=255)
                                                                 acaTGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTca                        >  2:4508971/1‑93 (MQ=255)
                                                                  caTGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcac                       <  2:830124/93‑1 (MQ=255)
                                                                    tGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                     >  2:1277691/1‑93 (MQ=255)
                                                                     gCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcaccac                    >  1:7162279/1‑93 (MQ=255)
                                                                       aGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCAccc                  <  1:5933472/93‑1 (MQ=255)
                                                                       aGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCAccc                  >  1:7606942/1‑93 (MQ=255)
                                                                        gACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCg                 >  2:2993146/1‑93 (MQ=255)
                                                                        gACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGGAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCg                 >  1:1483091/1‑93 (MQ=255)
                                                                         aCAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGa                >  1:3840041/1‑93 (MQ=255)
                                                                            ggTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAg             >  1:7497742/1‑93 (MQ=255)
                                                                            ggTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAg             >  2:2893482/1‑93 (MQ=255)
                                                                            ggTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCCCACCACCCGAAAg             >  1:6529511/1‑93 (MQ=255)
                                                                             gTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGc            >  1:4737171/1‑93 (MQ=255)
                                                                              tttCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCg           >  2:1320910/1‑93 (MQ=255)
                                                                                tCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGtt         >  1:6927098/1‑93 (MQ=255)
                                                                                tCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGtt         <  1:100637/93‑1 (MQ=255)
                                                                                  tACCGCAGCATGAATGCCTGCCTCGTGGCTGGCGTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTaa       >  1:5804458/1‑93 (MQ=255)
                                                                                   aCCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTaaa      <  2:6034241/93‑1 (MQ=255)
                                                                                   aCCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTaaa      <  2:4124869/93‑1 (MQ=255)
                                                                                    ccGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGGTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAc     <  1:1923524/93‑1 (MQ=255)
                                                                                    ccGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAc     <  2:4987696/93‑1 (MQ=255)
                                                                                    ccGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAc     >  2:2126165/1‑93 (MQ=255)
                                                                                      gcagcaTGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAccg   >  1:5144521/1‑93 (MQ=255)
                                                                                      gcagcaTGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAccg   >  1:7782622/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGC  >  NC_000913/862182‑862361

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