Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,097,506 G→C 6.0% 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 6.0% 262.0 / 7.8 85S137T (AGT→ACT) gpxhydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base G (43/35);  new base C (3/2);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

AGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCAT  >  NC_002947/2097364‑2097655
                                                                                                                                              |                                                                                                                                                     
ctgatctcGCTTTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCCCCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAg                                                                                                                                                <  1:733756/143‑1 (MQ=255)
  ttttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGt                                                                                                                                              >  2:249320/1‑150 (MQ=255)
      ccGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGTCCAGGt                                                                                                                                              >  2:1795004/1‑146 (MQ=255)
      ccGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGTCCAGGt                                                                                                                                              <  1:1795004/146‑1 (MQ=255)
       cGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAgc                                                                                                                                         <  1:249320/150‑1 (MQ=255)
         cTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAgcgc                                                                                                                                       >  2:2168198/1‑150 (MQ=255)
            tttCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTAt                                                                                                                                    <  1:119124/150‑1 (MQ=255)
               cGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcc                                                                                                                                 >  2:885626/1‑150 (MQ=255)
                gCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGccc                                                                                                                                <  2:1876043/150‑1 (MQ=255)
                 cAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcccc                                                                                                                               >  2:126173/1‑150 (MQ=255)
                      tCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccacca                                                                                                                          >  1:2100647/1‑150 (MQ=255)
                        gAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTg                                                                                                                                                     <  1:681503/121‑1 (MQ=255)
                        gAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTg                                                                                                                                                     >  2:681503/1‑121 (MQ=255)
                        gAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccaccacc                                                                                                                        >  2:239476/1‑150 (MQ=255)
                              aaCGGCCCGGGCGCTCACCCTCTGTTGGTCGAACTCAAGCAACGCGTGTCGGGTGTGTTGGGTAGCCGGAAGATCAAGGGGGACTTCACCAAGTTTCTGGTCGACCCGGCCAGGGGCCAGGTCAAGCGCTATGCCCCCACCACTaaagca                                                                                                                  >  1:356924/1‑146 (MQ=255)
                                 ggCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAg                                                                                                               <  2:327069/150‑1 (MQ=255)
                                   cccGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcc                                                                                                                                 >  1:1021385/1‑130 (MQ=255)
                                   cccGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcc                                                                                                                                 <  2:1021385/130‑1 (MQ=255)
                                            tCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGc                                                                                                    >  2:2179223/1‑150 (MQ=255)
                                            tCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGc                                                                                                    >  1:1930264/1‑150 (MQ=255)
                                              aCCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCaa                                                                                                                                           >  1:844268/1‑109 (MQ=255)
                                              aCCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCaa                                                                                                                                           <  2:844268/109‑1 (MQ=255)
                                                      ttCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAg                                                                                          >  2:878721/1‑150 (MQ=255)
                                                         gtcgAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGt                                                                                       >  1:11458/1‑150 (MQ=255)
                                                            gAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCt                                                                                                                                      <  1:644491/100‑1 (MQ=255)
                                                            gAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCt                                                                                                                                      >  2:644491/1‑100 (MQ=255)
                                                            gAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGc                                                                                         <  2:203522/145‑1 (MQ=255)
                                                            gAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGc                                                                                         >  1:203522/1‑145 (MQ=255)
                                                             aaCTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCACACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGc                                                                                   >  2:1685366/1‑150 (MQ=255)
                                                                  aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGGCATCGAGCGTTTGCTCAGc                                                                              >  1:39785/1‑150 (MQ=255)
                                                                  aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCt                                                                                  >  1:371001/1‑146 (MQ=255)
                                                                  aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCt                                                                                  <  2:371001/146‑1 (MQ=255)
                                                                  aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGc                                                                              >  2:1549505/1‑150 (MQ=255)
                                                                    gCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCa                                                                                                                                                 <  2:189291/81‑1 (MQ=255)
                                                                    gCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCa                                                                                                                                                 >  1:189291/1‑81 (MQ=255)
                                                                        cgcgCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGa                                                                        >  2:455284/1‑150 (MQ=255)
                                                                         gcgcCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGGCATCGAGCGTTTGCTCAGCCGCTGaa                                                                       <  2:39785/150‑1 (MQ=255)
                                                                           gcCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                                     <  1:2015020/150‑1 (MQ=255)
                                                                             cccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTAACTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGc                                                                   >  1:175155/1‑150 (MQ=255)
                                                                              ccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCg                                                                                                   >  2:2265944/1‑117 (MQ=255)
                                                                              ccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCg                                                                                                   <  1:2265944/117‑1 (MQ=255)
                                                                               cGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCAc                                                                 <  1:126173/150‑1 (MQ=255)
                                                                                ggggATTTTGGCTAGCCAGAAGATCAAGTGGACCTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAATCCACAGGCCCTGGCGGCGGACACCGAGCGTTTGCTCAGCCGCTGAAGGGCAcc                                                                >  2:473310/1‑150 (MQ=255)
                                                                                     ttttGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATg                                                           <  2:393853/150‑1 (MQ=255)
                                                                                     ttttGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTAAAGCGATATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATg                                                           <  2:80156/150‑1 (MQ=255)
                                                                                           gTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGTAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCAc                                                                 >  2:315503/1‑138 (MQ=255)
                                                                                           gTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCAc                                                                 <  1:315503/138‑1 (MQ=255)
                                                                                               ccAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAgg                                                 <  1:2168147/150‑1 (MQ=255)
                                                                                                 agaagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGc                                               <  1:2179223/150‑1 (MQ=255)
                                                                                                  gaagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGc                                                                              <  2:1766816/118‑1 (MQ=255)
                                                                                                  gaagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGc                                                                              >  1:1766816/1‑118 (MQ=255)
                                                                                                   aagcTCAAGTGGAACTTCACCAAGTTTCTGGTCGCCCAGGCCATTGATCAGTTCAAGCGCTATGCCCCCACCCAAAAGCCACAGGCCCTGGAGGCGGACATCNAGCGTTNGCTCAGCCGCTGNAGGGCACCCAATGGTCGACCAGGGCCa                                             <  1:1685366/150‑1 (MQ=255)
                                                                                                    agaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAg                                            <  2:2100647/150‑1 (MQ=255)
                                                                                                         aaGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGct                                       >  2:455367/1‑150 (MQ=255)
                                                                                                           gTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCACACCACCAAGCCACAGGCCCTGGGGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGctct                                     >  1:1169298/1‑150 (MQ=255)
                                                                                                            tGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGc                                                                                                                                  >  2:1487605/1‑56 (MQ=255)
                                                                                                            tGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGc                                                                                                                                  <  1:1487605/56‑1 (MQ=255)
                                                                                                              gAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTcg                                  >  1:1412539/1‑150 (MQ=255)
                                                                                                                aCTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGAc                                                     >  1:1791202/1‑129 (MQ=255)
                                                                                                                aCTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGAc                                                     <  2:1791202/129‑1 (MQ=255)
                                                                                                                 cTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcgg                               >  2:2084608/1‑150 (MQ=255)
                                                                                                                 cTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcgg                               <  2:11458/150‑1 (MQ=255)
                                                                                                                  ttCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcggc                              >  2:1408074/1‑150 (MQ=255)
                                                                                                                     aCCAAGTTTCTGGTCGACCCGGCCAATGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGa                           >  2:1337689/1‑150 (MQ=255)
                                                                                                                       cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGa                                                      <  1:213383/121‑1 (MQ=255)
                                                                                                                       cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGa                                                      >  2:213383/1‑121 (MQ=255)
                                                                                                                       cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAAt                         >  1:207016/1‑150 (MQ=255)
                                                                                                                          gTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGc                      <  2:1412539/150‑1 (MQ=255)
                                                                                                                           tttCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTc                                                        <  2:1863984/115‑1 (MQ=255)
                                                                                                                           tttCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTc                                                        >  1:1863984/1‑115 (MQ=255)
                                                                                                                             tCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTg                   <  1:2168198/150‑1 (MQ=255)
                                                                                                                               tGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCg                                                       >  2:307856/1‑112 (MQ=255)
                                                                                                                               tGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCg                                                       <  1:307856/112‑1 (MQ=255)
                                                                                                                               tGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc                 <  2:439987/150‑1 (MQ=255)
                                                                                                                                 gTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCa               >  2:2064636/1‑150 (MQ=255)
                                                                                                                                 cttgACCTGGCCACTGGCCGGGTCGAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCa               >  2:1672056/4‑150 (MQ=255)
                                                                                                                                  ttgACCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg              <  1:2047076/148‑1 (MQ=255)
                                                                                                                                  ttgACCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg              >  2:1406991/3‑150 (MQ=255)
                                                                                                                                    gACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATg                                                           >  1:885130/1‑103 (MQ=255)
                                                                                                                                    gACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATg                                                           <  2:885130/103‑1 (MQ=255)
                                                                                                                                    gACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGt            >  2:1026581/1‑150 (MQ=255)
                                                                                                                                     aCCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTa           <  1:455284/150‑1 (MQ=255)
                                                                                                                                       ctggCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAAt                         >  1:1562106/3‑134 (MQ=255)
                                                                                                                                       ctggCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAAt                         <  2:1562106/132‑1 (MQ=255)
                                                                                                                                         ggCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACACCGAGCGTTTGCTCAGCCGCTGACGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGGCg       >  2:269432/1‑150 (MQ=255)
                                                                                                                                              gTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCGGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTc    >  1:666633/1‑148 (MQ=255)
                                                                                                                                              gTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCGGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTc    <  2:666633/148‑1 (MQ=255)
                                                                                                                                              gTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGACGGGTAGCGCAAGACAATGCTCCGCTAGGAGTGCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGGTCAt  >  1:273994/1‑150 (MQ=255)
                                                                                                                                              gTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCAt  <  1:1549505/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
AGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCCTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTCAT  >  NC_002947/2097364‑2097655

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: