Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,097,506 G→C 5.8% S137T (AGT→ACT)  gpx → hydroperoxy fatty acid reductase Gpx1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,097,5060GC 5.8% 222.2 / 4.5 69S137T (AGT→ACT) gpxhydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base C (2/2);  total (34/35)
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.92e-01

GCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCA  >  NC_002947/2097359‑2097654
                                                                                                                                                   |                                                                                                                                                    
gCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTg                                                                                                                                                    <  1:2261900/150‑1 (MQ=255)
  gTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGc                                                                                                                                                  <  2:461293/150‑1 (MQ=255)
     aGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTCTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAg                                                                                                                                               <  2:1400424/150‑1 (MQ=255)
      gTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAgg                                                                                                                                              >  2:1929814/1‑150 (MQ=255)
            cGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCTGGCCACTGGCCGGGTcgacc                                                                                                                                        <  1:53103/150‑5 (MQ=255)
            cGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTATGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAgc                                                                                                                                        >  1:436109/1‑150 (MQ=255)
              cTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGATCAACATCTGTTCGTCGCACTCAAGCAACTCGCCCCGTGTATTTTTTTTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAgcgc                                                                                                                                      >  2:1465363/1‑150 (MQ=255)
                   tCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGc                                                                                                                                 >  1:820540/1‑150 (MQ=255)
                    cGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTTGCCAGAAGATCTATGGGTACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcc                                                                                                                                <  1:1486819/150‑1 (MQ=255)
                     gCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGccc                                                                                                                               >  1:1823785/1‑150 (MQ=255)
                     gCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGccc                                                                                                                               >  1:1941417/1‑150 (MQ=255)
                             gAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTg                                                                                                                                                    <  2:1433104/121‑1 (MQ=255)
                             gAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTg                                                                                                                                                    >  1:1433104/1‑121 (MQ=255)
                                gTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAg                                                                                                                    >  2:1922779/1‑150 (MQ=255)
                                     cGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGCCTACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCcaca                                                                                                               >  1:1970458/1‑150 (MQ=255)
                                         ccGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGcc                                                                                                           >  2:964732/1‑150 (MQ=255)
                                          cGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGccc                                                                                                          >  2:1922936/1‑150 (MQ=255)
                                           gggCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCt                                                                                                         <  2:1651910/150‑1 (MQ=255)
                                                 tCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccaccac                                                                                                                        >  2:752509/1‑129 (MQ=255)
                                                 tCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccaccac                                                                                                                        <  1:752509/129‑1 (MQ=255)
                                                 tCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGc                                                                                                   <  1:1289804/150‑1 (MQ=255)
                                                 tCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCACGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGc                                                                                                   >  2:69860/1‑150 (MQ=255)
                                                     ccTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCaa                                                                                                                     >  1:1562585/1‑128 (MQ=255)
                                                     ccTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCaa                                                                                                                     <  2:1562585/128‑1 (MQ=255)
                                                            tcgtcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGc                                                                                        <  1:1399287/150‑1 (MQ=255)
                                                               tcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGa                                                                                                <  1:1972113/139‑1 (MQ=255)
                                                               tcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGa                                                                                                >  2:1972113/1‑139 (MQ=255)
                                                               tcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGtt                                                                                     <  1:2371409/150‑1 (MQ=255)
                                                                cgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGttt                                                                                    <  2:1571610/150‑1 (MQ=255)
                                                                cgAACTCAAGCAACGCGCCCCGGGTATTTGGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGttt                                                                                    <  2:834185/150‑1 (MQ=255)
                                                                     tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGccc                                                                                                                               >  2:979338/1‑102 (MQ=255)
                                                                     tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGccc                                                                                                                               <  1:979338/102‑1 (MQ=255)
                                                                     tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCa                                                                               <  2:720654/150‑1 (MQ=255)
                                                                      cAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAg                                                                              <  2:1823785/150‑1 (MQ=255)
                                                                       aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGc                                                                             >  1:1846214/1‑150 (MQ=255)
                                                                           gaCGCGCCCCGGGGATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCGATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCGGCt                                                                         >  2:1001967/2‑150 (MQ=255)
                                                                             cgcgCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGa                                                                       <  1:1922936/150‑1 (MQ=255)
                                                                             cgcgCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGa                                                                       <  2:2161346/150‑1 (MQ=255)
                                                                              gcgcCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGaa                                                                      <  1:1568420/150‑1 (MQ=255)
                                                                               cgcCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAg                                                                     <  2:568177/150‑1 (MQ=255)
                                                                                  cccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAg                                                                                                     >  1:791202/1‑115 (MQ=255)
                                                                                  cccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAg                                                                                                     <  2:791202/115‑1 (MQ=255)
                                                                                   ccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGa                                                                                                >  1:666753/1‑119 (MQ=255)
                                                                                   ccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGa                                                                                                <  2:666753/119‑1 (MQ=255)
                                                                                   ccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCa                                                                 >  2:1377021/1‑150 (MQ=255)
                                                                                    cGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCAc                                                                <  2:245705/150‑1 (MQ=255)
                                                                                                gtggCCGGAAGATCAAGGGGAATTTCACCACGTTTCTGGTCGACCCGTGCAGTGGCGAGGTCAAGCGCTATGCCCCGACGACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACTCGATGGTCGAc                                                    <  2:229107/147‑1 (MQ=255)
                                                                                                     cAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAggg                                               >  1:174706/1‑150 (MQ=255)
                                                                                                      agaagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGc                                              >  2:2009544/1‑150 (MQ=255)
                                                                                                         agaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAg                                           <  1:1377021/150‑1 (MQ=255)
                                                                                                          gaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCGGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGCAGGGAACCCAATGGTTGACCAGGGCCAGt                                          <  2:1970458/150‑1 (MQ=255)
                                                                                                                gTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGctct                                    >  1:859927/1‑150 (MQ=255)
                                                                                                                      cTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcgg                              <  2:2362146/150‑1 (MQ=255)
                                                                                                                          aCCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGa                          <  2:1008381/150‑1 (MQ=255)
                                                                                                                          aCCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGa                          >  2:1531696/1‑150 (MQ=255)
                                                                                                                            cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTAt                                                                                                                                   >  2:1981693/1‑43 (MQ=255)
                                                                                                                            cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTAt                                                                                                                                   <  1:1981693/43‑1 (MQ=255)
                                                                                                                             aaGTTTCTGGGTTACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATg                       <  1:1922779/150‑1 (MQ=255)
                                                                                                                              aGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATgg                      <  2:1979753/150‑1 (MQ=255)
                                                                                                                                      gTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCa              >  2:2406774/1‑150 (MQ=255)
                                                                                                                                      gTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCa              >  2:2201311/1‑150 (MQ=255)
                                                                                                                                      gTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCa              >  2:2121085/1‑150 (MQ=255)
                                                                                                                                       ttgACCTGGCCACTGGCCGGGTCGAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATc                                                                                            <  2:527254/69‑1 (MQ=25)
                                                                                                                                       ttgACCTGGCCACTGGCCGGGTCGAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATc                                                                                            >  1:527254/3‑71 (MQ=25)
                                                                                                                                         gACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCAGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGt           >  2:1801585/1‑150 (MQ=255)
                                                                                                                                          aCCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTa          >  2:1918775/1‑150 (MQ=255)
                                                                                                                                            ccGNCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGt        >  2:1505923/1‑150 (MQ=255)
                                                                                                                                               gCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCCATGCCCTGTCTGGTCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGt     >  2:63840/1‑150 (MQ=255)
                                                                                                                                                  aGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCa  <  1:2201311/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCA  >  NC_002947/2097359‑2097654

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