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,2720TC98.0% 306.3 / ‑2.5 99T164A (ACT→GCT) ybiYputative pyruvate formate‑lyase activating enzyme 2
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (43/54);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 105 total aligned reads displayed.

AACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGC  >  NC_000913/862180‑862361
                                                                                            |                                                                                         
aaCGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGt                                                                                           <  2:1580248/93‑1 (MQ=255)
  cGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCAt                                                                                         >  1:798269/1‑93 (MQ=255)
  cGTGTTTGAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCAt                                                                                         <  2:2884751/93‑1 (MQ=255)
  agtgtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCAt                                                                                         <  1:1531881/92‑1 (MQ=255)
   gtgtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATg                                                                                        >  2:3089034/1‑93 (MQ=255)
    tgtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGa                                                                                       <  2:2467533/93‑1 (MQ=255)
     gtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGaa                                                                                      <  2:4275704/93‑1 (MQ=255)
     gtTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGaa                                                                                      <  2:3752324/93‑1 (MQ=255)
      ttttAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAt                                                                                     <  1:3758194/93‑1 (MQ=255)
         tAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcc                                                                                  <  1:2294616/93‑1 (MQ=255)
         tAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcc                                                                                  <  2:4113770/93‑1 (MQ=255)
         tAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcc                                                                                  <  1:2656657/93‑1 (MQ=255)
          aaaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcct                                                                                 <  2:3330446/93‑1 (MQ=255)
          aaaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcct                                                                                 <  2:217359/93‑1 (MQ=255)
          aaaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcct                                                                                 >  1:2955261/1‑93 (MQ=255)
           aaTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctg                                                                                >  1:2820592/1‑93 (MQ=255)
            aTCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgc                                                                               >  2:290559/1‑93 (MQ=255)
             tCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgcc                                                                              <  1:872948/93‑1 (MQ=255)
             tCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAAtgcctgcc                                                                              <  1:4569641/93‑1 (MQ=255)
              cGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGCATgcctgcct                                                                             >  2:58717/1‑93 (MQ=255)
              cGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATgcctgcct                                                                             <  2:75490/93‑1 (MQ=255)
               ggCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTc                                                                            >  1:4416744/1‑93 (MQ=255)
                gCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCg                                                                           >  1:4140701/1‑93 (MQ=255)
                 cAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGt                                                                          <  1:129927/93‑1 (MQ=255)
                 cAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGt                                                                          <  2:2215145/93‑1 (MQ=255)
                  aaGAAACAGATCGATATAGGGGGGAGAAGGGGCGATATATTTCCACGGCACACGCAGACAGGCTTCTACCGCAGCAGGAATGACTTCCTCGtg                                                                         >  2:667434/1‑93 (MQ=255)
                  aaGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtg                                                                         >  2:667435/1‑93 (MQ=255)
                  aaGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtg                                                                         >  2:274654/1‑93 (MQ=255)
                   aGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtgg                                                                        >  2:2241898/1‑93 (MQ=255)
                   aGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtgg                                                                        <  1:2336981/93‑1 (MQ=255)
                     aaaCAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtggct                                                                      >  1:2333890/1‑93 (MQ=255)
                     aaaCAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtggct                                                                      >  2:2719322/1‑93 (MQ=255)
                      aaCAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtggctg                                                                     >  2:2521725/1‑93 (MQ=255)
                      aaCAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGCGTTCCAGCGCCGCATGAATGCCTGCCTCGtggatg                                                                     >  2:2642801/1‑93 (MQ=255)
                       aCAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGtggctgg                                                                    >  1:2114166/1‑93 (MQ=255)
                          gATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCtt                                                                 >  2:3911875/1‑93 (MQ=255)
                           aTCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTg                                                                >  1:2285935/1‑93 (MQ=255)
                            tCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGc                                                               <  2:2992734/93‑1 (MQ=255)
                              gATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAg                                                             <  2:2106435/93‑1 (MQ=255)
                              gATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAg                                                             <  1:1941672/93‑1 (MQ=255)
                               atatAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGt                                                            >  1:4415703/1‑93 (MQ=255)
                                tataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTa                                                           <  2:4386658/93‑1 (MQ=255)
                                 ataGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAg                                                          <  1:4746377/93‑1 (MQ=255)
                                  taGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgc                                                         <  1:2452406/93‑1 (MQ=255)
                                  taGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgc                                                         >  2:1450156/1‑93 (MQ=255)
                                  taGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgc                                                         >  2:2955933/1‑93 (MQ=255)
                                    gggCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                       >  1:685636/1‑93 (MQ=255)
                                      gCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgcca                                                     <  2:3156995/93‑1 (MQ=255)
                                       cAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccat                                                    <  2:196268/93‑1 (MQ=255)
                                       cAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccat                                                    <  2:3974087/93‑1 (MQ=255)
                                       cAGAGAAGGGGCGATATATTGCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccat                                                    <  2:3974088/93‑1 (MQ=255)
                                        agagAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatc                                                   >  2:1356743/1‑93 (MQ=255)
                                        agagAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatc                                                   <  1:33694/93‑1 (MQ=255)
                                         gagaAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcg                                                  >  2:2893312/1‑93 (MQ=255)
                                         gagaAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcg                                                  >  1:4448241/1‑93 (MQ=255)
                                           gaAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgcc                                                >  1:2692761/1‑93 (MQ=255)
                                           gaAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgcc                                                >  2:573185/1‑93 (MQ=255)
                                            aaGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGcgccatcgcca                                               <  1:4525980/93‑1 (MQ=255)
                                              ggggCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCAtt                                             >  2:1962724/1‑93 (MQ=255)
                                              ggggCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCAtt                                             >  1:3864434/1‑93 (MQ=255)
                                                ggCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTc                                           >  1:1037012/1‑93 (MQ=255)
                                                 gCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTcc                                          >  1:644905/1‑93 (MQ=255)
                                                      atatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                     <  1:3488693/93‑1 (MQ=255)
                                                       tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGc                                    >  1:3812488/1‑93 (MQ=255)
                                                       tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGc                                    >  2:2428561/1‑93 (MQ=255)
                                                       tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGc                                    >  1:3925005/1‑93 (MQ=255)
                                                         tttCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCAt                                  >  1:3826540/1‑93 (MQ=255)
                                                         tttCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCAt                                  <  2:3959415/93‑1 (MQ=255)
                                                         tttCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCAt                                  >  2:1209666/1‑93 (MQ=255)
                                                         ttgccACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCAt                                  <  2:3959414/90‑1 (MQ=255)
                                                          ttCCACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATa                                 <  2:3014024/93‑1 (MQ=255)
                                                            ccACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATaaa                               <  2:2365624/93‑1 (MQ=255)
                                                             cACGGCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAg                              >  2:605285/1‑93 (MQ=255)
                                                                ggCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGc                           >  2:2738742/1‑93 (MQ=255)
                                                                 gCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCt                          >  2:1525671/1‑93 (MQ=255)
                                                                 gCACATGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCt                          <  1:426902/93‑1 (MQ=255)
                                                                    caTGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcac                       >  2:3154276/1‑93 (MQ=255)
                                                                     aTGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacc                      >  2:4607533/1‑93 (MQ=255)
                                                                     aTGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGCCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacc                      >  2:1023661/1‑93 (MQ=255)
                                                                      tGCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCCTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                     <  1:4346658/93‑1 (MQ=255)
                                                                       gCAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcaccac                    >  2:531471/1‑93 (MQ=255)
                                                                        cAGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCaccacc                   <  2:291399/93‑1 (MQ=255)
                                                                         aGACAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCAccc                  <  1:2335331/93‑1 (MQ=255)
                                                                            cAGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGaa               >  2:429611/1‑93 (MQ=255)
                                                                             aGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGaaa              <  2:2194048/93‑1 (MQ=255)
                                                                             aGGTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGaaa              <  2:1692204/93‑1 (MQ=255)
                                                                              ggTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAg             <  2:1108749/93‑1 (MQ=255)
                                                                              ggTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAg             <  1:3901383/93‑1 (MQ=255)
                                                                              ggTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAg             <  1:1931992/93‑1 (MQ=255)
                                                                               gTTTCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGc            <  1:2684844/93‑1 (MQ=255)
                                                                                tttCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCg           <  1:1891337/93‑1 (MQ=255)
                                                                                tttCTACCGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCg           >  1:1909152/1‑93 (MQ=255)
                                                                                    tACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTaa       <  2:4111480/93‑1 (MQ=255)
                                                                                      ccGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAc     >  2:1572026/1‑93 (MQ=255)
                                                                                       cGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAcc    <  1:3963181/93‑1 (MQ=255)
                                                                                       cGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAcc    <  1:2434166/93‑1 (MQ=255)
                                                                                       cGCAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAcc    <  1:913160/93‑1 (MQ=255)
                                                                                        gcagcaTGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAccg   >  2:2483290/1‑93 (MQ=255)
                                                                                        gcagcaTGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAcc  >  2:705377/1‑92 (MQ=255)
                                                                                         cAGCATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAAccgc  <  2:516283/93‑1 (MQ=255)
                                                                                            |                                                                                         
AACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGC  >  NC_000913/862180‑862361

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