Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,551,460 T→C A72A (GCT→GCC hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,551,4600TC99.3% 935.2 / ‑4.7 280A72A (GCT→GCChemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (129/149);  total (131/149)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

Only 100 of 300 total aligned reads displayed.

CGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGT  >  NC_000913/2551370‑2551549
                                                                                          |                                                                                         
cggcggGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTc                                                                                         <  2:7094345/93‑1 (MQ=255)
  gcgggcggAGTCGGGTGTTGCGTAATGGTGGTGGTTTCGAAAACGCAGGGGTCAACTTTTCGCATGTCCACGGTGAGGCGATGTCTGCCTccc                                                                                       >  2:8455791/1‑92 (MQ=255)
   cgggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTCCCTCCGc                                                                                      >  2:4477433/1‑93 (MQ=255)
    gggcgTAGTCGGGTGTTGCGTAATGGTGATGTTTTCGAACAGGCAGGCGTCACCTTTTCGCTTGTCCACGGTGAGGCGATGCCTGCCTCCGcc                                                                                     <  2:5764205/93‑1 (MQ=255)
     ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCa                                                                                    >  1:5833395/1‑93 (MQ=255)
      gcgTAGTCGGGTGTTTCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCAc                                                                                   <  1:2628401/93‑1 (MQ=255)
      gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCAc                                                                                   <  1:239096/93‑1 (MQ=255)
      gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCAc                                                                                   >  2:4251952/1‑93 (MQ=255)
         tAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGc                                                                                >  2:3085846/1‑93 (MQ=255)
         gaGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGc                                                                                <  1:5560427/92‑1 (MQ=255)
          tgTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCt                                                                               <  2:7492525/92‑1 (MQ=255)
          aGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCt                                                                               >  1:1439795/1‑93 (MQ=255)
          aGTCGGGTGGTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCt                                                                               <  2:2831765/93‑1 (MQ=255)
           gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTc                                                                              >  1:2984483/1‑93 (MQ=255)
           gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTc                                                                              >  1:5684295/1‑93 (MQ=255)
           gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTc                                                                              >  2:6450773/1‑93 (MQ=255)
           gTCGGGGGTTGCGTAATGGAGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGACAGCCTCCGCCACCGCTc                                                                              >  2:4601249/1‑93 (MQ=255)
            tCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCa                                                                             >  2:4866182/1‑93 (MQ=255)
            tCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCa                                                                             <  1:610398/93‑1 (MQ=255)
             cGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCAt                                                                            <  1:4831474/93‑1 (MQ=255)
             cGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCAt                                                                            <  1:4662844/93‑1 (MQ=255)
              gggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATc                                                                           >  2:4142755/1‑93 (MQ=255)
              gggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATc                                                                           >  2:5710934/1‑93 (MQ=255)
                gtgtTGGGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGc                                                                         <  1:6755230/93‑1 (MQ=255)
                gtgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGc                                                                         <  2:4067981/93‑1 (MQ=255)
                gtgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGc                                                                         >  2:1628362/1‑93 (MQ=255)
                   ttGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTTCCTCCGCCACCGCTCATCGCCCg                                                                      >  2:6330352/1‑93 (MQ=255)
                   ttGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCg                                                                      <  2:2771354/93‑1 (MQ=255)
                   ttGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCg                                                                      <  2:2115938/93‑1 (MQ=255)
                      cGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGaa                                                                   <  2:3908650/93‑1 (MQ=255)
                        tAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACt                                                                 >  2:1227739/1‑93 (MQ=255)
                         aaTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACtt                                                                >  1:1372874/1‑93 (MQ=255)
                          aTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTg                                                               >  1:368998/1‑93 (MQ=255)
                             gtggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCg                                                            <  1:5299427/93‑1 (MQ=255)
                              tggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCgg                                                           >  2:1175734/1‑93 (MQ=255)
                               ggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCggg                                                          <  1:6556448/93‑1 (MQ=255)
                                 tgtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcg                                                        >  2:2520957/1‑93 (MQ=255)
                                  gtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcgc                                                       >  2:7280849/1‑93 (MQ=255)
                                   ttttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcgcg                                                      >  1:1182717/1‑92 (MQ=255)
                                        gAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTc                                                 <  2:5588907/93‑1 (MQ=255)
                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                               <  2:5911188/93‑1 (MQ=255)
                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                               <  2:2271699/93‑1 (MQ=255)
                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                               <  2:5055069/93‑1 (MQ=255)
                                            aggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCCCCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAgg                                             <  2:6891194/93‑1 (MQ=255)
                                             ggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGc                                            >  1:3625111/1‑93 (MQ=255)
                                             ggcaggcCTCAACTTTTCGCATGTCCCCGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGc                                            <  2:4316649/93‑1 (MQ=255)
                                              gcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCg                                           <  2:3673623/93‑1 (MQ=255)
                                                aggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGAt                                         <  1:2748256/93‑1 (MQ=255)
                                                aggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGAt                                         <  1:5910496/93‑1 (MQ=255)
                                                 ggcGTCAACTTTTCGCATGTCCACGGGGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATg                                        <  1:315554/93‑1 (MQ=255)
                                                  gcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                       <  2:3470375/93‑1 (MQ=255)
                                                  gcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                       <  1:5437074/93‑1 (MQ=255)
                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                      <  2:3969180/93‑1 (MQ=255)
                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                      <  2:8317705/93‑1 (MQ=255)
                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                      <  2:2504402/93‑1 (MQ=255)
                                                    gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgggg                                     >  2:6681214/1‑92 (MQ=255)
                                                    gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGc                                     <  1:1547130/93‑1 (MQ=255)
                                                    gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGc                                     >  2:6034256/1‑93 (MQ=255)
                                                     tCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCg                                    >  2:4528128/1‑93 (MQ=255)
                                                        aCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGttt                                 <  1:7067073/93‑1 (MQ=255)
                                                         cTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTc                                <  2:1748765/93‑1 (MQ=255)
                                                          gattCGCATGTCCACGGTGTGGCGATGCCCGCTTCCGCCACCGCTCAACGCCCGGAACTTGCCGGGCGCAGTTTCTAGGCGATGGGCGTTTCa                               >  1:7834856/3‑93 (MQ=255)
                                                           tttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                              <  1:4773956/93‑1 (MQ=255)
                                                           tttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                              <  2:8103786/93‑1 (MQ=255)
                                                           tttCGCAGGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                              <  1:7197518/93‑1 (MQ=255)
                                                             tCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTg                            >  2:6935118/1‑93 (MQ=255)
                                                             tCGCATGTCCACGGTGACGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTg                            <  2:3272427/93‑1 (MQ=255)
                                                              cGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTgg                           <  1:2945692/93‑1 (MQ=255)
                                                               gCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                          <  1:966478/93‑1 (MQ=255)
                                                                cATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTa                         >  1:1958061/1‑93 (MQ=255)
                                                                cATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTa                         <  2:7804315/93‑1 (MQ=255)
                                                                cATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTa                         <  1:7421607/93‑1 (MQ=255)
                                                                cATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTa                         >  1:5285492/1‑93 (MQ=255)
                                                                 aTGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAg                        >  2:6771346/1‑93 (MQ=255)
                                                                 aTGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAg                        >  2:8655806/1‑93 (MQ=255)
                                                                 aTGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAg                        >  2:4129509/1‑93 (MQ=255)
                                                                   gTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTg                      >  1:2109208/1‑93 (MQ=255)
                                                                      cACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCAt                   >  1:4092564/1‑93 (MQ=255)
                                                                       aCGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATc                  <  1:3468435/93‑1 (MQ=255)
                                                                       aCGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATc                  >  2:7614900/1‑93 (MQ=255)
                                                                        cGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATcc                 <  1:191574/93‑1 (MQ=255)
                                                                         ggTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCg                >  2:3779536/1‑93 (MQ=255)
                                                                         ggTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCg                <  1:1692797/93‑1 (MQ=255)
                                                                          gTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGc               >  1:6301784/1‑93 (MQ=255)
                                                                          gTGAGGCGATGCCTGCCTCCGCCACCACTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGc               >  1:3792046/1‑93 (MQ=255)
                                                                            gAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCAt             >  1:3287503/1‑93 (MQ=255)
                                                                              ggCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATaa           <  1:6336021/93‑1 (MQ=255)
                                                                                cGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAcc         <  2:5460920/93‑1 (MQ=255)
                                                                                cGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAcc         >  2:7930692/1‑93 (MQ=255)
                                                                                 gATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc        >  2:5315370/1‑93 (MQ=255)
                                                                                  aTGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCg       >  1:8131198/1‑93 (MQ=255)
                                                                                   tgcctgccTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGt      <  1:8692623/93‑1 (MQ=255)
                                                                                    gcctgcctCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTa     >  1:8324488/1‑93 (MQ=255)
                                                                                    gcctgcctCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTa     <  2:8032016/93‑1 (MQ=255)
                                                                                    gcctgcctCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTa     <  1:1012327/93‑1 (MQ=255)
                                                                                      cTGCCTCCGCCCGCGCTCATCGCCCGGAACTTGCCGTGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATg   <  2:302489/93‑1 (MQ=255)
                                                                                      cTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATg   >  1:509982/1‑93 (MQ=255)
                                                                                       tGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGt  >  1:8530363/1‑93 (MQ=255)
                                                                                       tGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGt  >  1:6783354/1‑93 (MQ=255)
                                                                                       tGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGt  <  1:7572768/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGT  >  NC_000913/2551370‑2551549

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