Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,556,088 G→A R204R (CGC→CGT b2441 ← ethanolamine ammonia‑lyase, heavy chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,556,0880GA99.0% 722.8 / ‑5.8 208R204R (CGC→CGTb2441ethanolamine ammonia‑lyase, heavy chain
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (114/91);  total (114/94)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 217 total aligned reads displayed.

CGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGATTTGCGC  >  NC_000913/2555998‑2556177
                                                                                          |                                                                                         
  cTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGc                                                                                       >  2:1378654/1‑93 (MQ=255)
  cTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGc                                                                                       >  1:5283566/1‑93 (MQ=255)
   ttCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATTACGCCATAGATGGTATCCAACACACGGCt                                                                                      >  2:2552413/1‑93 (MQ=255)
   ttCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCt                                                                                      >  2:5183795/1‑93 (MQ=255)
   ttCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCt                                                                                      <  2:4251678/93‑1 (MQ=255)
   ttCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCt                                                                                      >  2:5289920/1‑93 (MQ=255)
   ttCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCt                                                                                      <  2:5557870/93‑1 (MQ=255)
   ttCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCt                                                                                      >  1:5078226/1‑93 (MQ=255)
    tCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCtt                                                                                     >  1:1678684/1‑93 (MQ=255)
      gATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTaa                                                                                   <  2:3949562/93‑1 (MQ=255)
      gATCTGGGTGGTGACGTGCGCCAATACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTaa                                                                                   <  2:5393671/93‑1 (MQ=255)
      gATCTGGGTGATGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTaa                                                                                   <  2:1255232/93‑1 (MQ=255)
       aTCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGCCGATCACGCCATAGATGGTATCCAACACACGGCTTAAg                                                                                  >  1:6423603/1‑93 (MQ=255)
        tCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCAAACACACGGCTTAAGt                                                                                 >  2:7203727/1‑93 (MQ=255)
         cTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGtt                                                                                >  1:2353780/1‑93 (MQ=255)
          tGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGttt                                                                               <  2:3182429/93‑1 (MQ=255)
             gtggtgACGTGCACCAGTACACAGCCCTGAGTTGGAATGTTGAATTTGTCGATCACGCCATAAATGGTATCCAACACACGGCTTAAGTTTTcc                                                                            <  1:6021452/93‑1 (MQ=255)
              tggtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCa                                                                           >  2:3765454/1‑93 (MQ=255)
              tggtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCa                                                                           <  2:5705848/93‑1 (MQ=255)
              gggtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCa                                                                           <  2:4827202/92‑1 (MQ=255)
               ggtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAc                                                                          <  1:6937777/93‑1 (MQ=255)
               ggtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAc                                                                          <  2:6929202/93‑1 (MQ=255)
               ggtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAc                                                                          <  1:2055359/93‑1 (MQ=255)
                gtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAcg                                                                         >  2:1562281/1‑93 (MQ=255)
                gtgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAcg                                                                         >  2:3265814/1‑93 (MQ=255)
                 tgACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAcgt                                                                        <  2:2908670/93‑1 (MQ=255)
                  gACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGCTCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAcgtc                                                                       >  1:5671195/1‑93 (MQ=255)
                  gACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAcgtc                                                                       >  2:4528281/1‑93 (MQ=255)
                  gACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAcgtc                                                                       <  2:1174015/93‑1 (MQ=255)
                    cGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCAcgtcgt                                                                     >  1:2103641/1‑93 (MQ=255)
                      tGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCgtca                                                                   <  1:6328003/93‑1 (MQ=255)
                      tGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCgtca                                                                   <  2:7861891/93‑1 (MQ=255)
                        cgcCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCTATCACGCCATAGATGGTATCCACCACACGGCTTAAGTTTTCCACGTCgtcagt                                                                 <  2:7886984/93‑1 (MQ=255)
                         gcCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCgtcagtc                                                                <  1:7457064/93‑1 (MQ=255)
                          ccAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCgtcagtca                                                               >  2:2816352/1‑93 (MQ=255)
                          ccAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCgtcagtca                                                               <  2:6075347/93‑1 (MQ=255)
                          ccAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCgtcagtca                                                               >  2:7700153/1‑93 (MQ=255)
                           cAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCAc                                                              <  1:1054098/93‑1 (MQ=255)
                            aGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCAcc                                                             >  1:6533414/1‑93 (MQ=255)
                             gTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCg                                                            >  2:6188493/1‑93 (MQ=255)
                             gTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCg                                                            <  1:3158855/93‑1 (MQ=255)
                              tACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCgg                                                           >  1:1779757/1‑93 (MQ=255)
                               acacAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACCCGGCTTAAGTTTTCCACGTCGTCAGTCACCggg                                                          <  2:6588733/93‑1 (MQ=255)
                               acacAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCggg                                                          >  2:2293280/1‑93 (MQ=255)
                                  caGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTa                                                       <  1:6095739/93‑1 (MQ=255)
                                   aGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTaa                                                      <  2:1853476/93‑1 (MQ=255)
                                   aGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTaa                                                      >  1:7597668/1‑93 (MQ=255)
                                    gCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAAc                                                     <  1:3042303/93‑1 (MQ=255)
                                    gCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAAc                                                     >  2:5878754/1‑93 (MQ=255)
                                      ccTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGc                                                   >  1:2306908/1‑93 (MQ=255)
                                      ccTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGc                                                   >  1:2501674/1‑93 (MQ=255)
                                        tGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCg                                                 >  1:1071847/1‑93 (MQ=255)
                                          aGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGAt                                               <  2:4569293/93‑1 (MQ=255)
                                          aGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGAt                                               >  1:471420/1‑93 (MQ=255)
                                          aGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGAt                                               >  1:2231036/1‑93 (MQ=255)
                                           gTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATc                                              >  2:5997239/1‑93 (MQ=255)
                                            ttGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATcc                                             >  1:7115983/1‑92 (MQ=255)
                                            ttGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCa                                             >  1:4770895/1‑93 (MQ=255)
                                             tGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCAc                                            >  2:4943046/1‑93 (MQ=255)
                                             tGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCAc                                            <  1:5340575/93‑1 (MQ=255)
                                              gggATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCAcc                                           >  2:1031648/1‑93 (MQ=255)
                                                  tGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCAt                                       >  2:2059409/1‑93 (MQ=255)
                                                  tGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCAt                                       <  1:7035869/93‑1 (MQ=255)
                                                  tGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCAt                                       >  1:4936409/1‑93 (MQ=255)
                                                  tGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCAt                                       >  2:3788614/1‑93 (MQ=255)
                                                  tGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCAt                                       <  1:2698716/93‑1 (MQ=255)
                                                   gTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATc                                      >  2:3376661/1‑93 (MQ=255)
                                                    ttGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCg                                     >  2:1188543/1‑93 (MQ=255)
                                                    ttGAATTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCg                                     <  1:7534501/93‑1 (MQ=255)
                                                       aaTTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGccc                                  >  1:4088671/1‑93 (MQ=255)
                                                        aTTTGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGTATCGCCCa                                 >  2:2078575/1‑93 (MQ=255)
                                                          ttGTCGATCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCAcc                               >  1:944687/1‑93 (MQ=255)
                                                                aTCACGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAg                         >  2:1635752/1‑93 (MQ=255)
                                                                   ccGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGaa                      <  1:3476782/92‑1 (MQ=255)
                                                                   aCGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGaa                      >  2:845973/1‑93 (MQ=255)
                                                                   aCGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGaa                      <  2:6964772/93‑1 (MQ=255)
                                                                    ggCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGaaa                     <  1:551887/92‑1 (MQ=255)
                                                                    cGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGaaa                     >  2:263914/1‑93 (MQ=255)
                                                                    cGCCATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGaaa                     <  2:6141312/93‑1 (MQ=255)
                                                                      ccATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGc                   <  1:6859202/93‑1 (MQ=255)
                                                                      ccATAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGc                   >  2:3023997/1‑93 (MQ=255)
                                                                        aTAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGccc                 <  1:5586027/93‑1 (MQ=255)
                                                                         tAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCt                >  1:486786/1‑93 (MQ=255)
                                                                         tAGATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCt                >  2:7322568/1‑93 (MQ=255)
                                                                           gATGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTc              >  2:6249367/1‑93 (MQ=255)
                                                                            aTGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCg             <  1:4238525/93‑1 (MQ=255)
                                                                             tGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGt            >  1:3551529/1‑93 (MQ=255)
                                                                             tGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGt            <  1:7203797/93‑1 (MQ=255)
                                                                             tGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGt            >  1:7335595/1‑93 (MQ=255)
                                                                             tGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGt            >  1:1697144/1‑93 (MQ=255)
                                                                             tGGTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGt            <  1:3225783/93‑1 (MQ=255)
                                                                               gTATCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAg          >  1:194985/1‑93 (MQ=255)
                                                                                tATCCAACACACGGCTTAAGTTTTCCACGTCCTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGa         >  2:4679772/1‑93 (MQ=255)
                                                                                  tCCAACACACGGCTTAAGTTTTTCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATTGCCCACCCCGAAGGAAACCCCTTCGTAGCtt       >  1:4951289/1‑93 (MQ=255)
                                                                                  tCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGAtt       >  2:6348261/1‑93 (MQ=255)
                                                                                  tCCAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGAtt       >  2:7273683/1‑93 (MQ=255)
                                                                                   ccAACACACGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGAttt      <  1:6669178/93‑1 (MQ=255)
                                                                                      acacacGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGATTTgcg   >  1:4314321/1‑93 (MQ=255)
                                                                                       cacaCGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGATTTgcgc  >  2:7691139/1‑93 (MQ=255)
                                                                                       cacaCGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGATTTgcgc  >  2:2706757/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGCTTCGATCTGGGTGGTGACGTGCGCCAGTACACAGCCCTGAGTTGGGATGTTGAATTTGTCGATCACGCCATAGATGGTATCCAACACGCGGCTTAAGTTTTCCACGTCGTCAGTCACCGGGTTAACGCCGATCACCGCATCGCCCACCCCGAAGGAAAGCCCTTCGTAGATTTGCGC  >  NC_000913/2555998‑2556177

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: