Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,453,830 A→G intergenic (‑188/‑97) paaZ ← / → paaA oxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase/ring 1,2‑phenylacetyl‑CoA epoxidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,453,8300AG100.0% 225.2 / NA 66intergenic (‑188/‑97)paaZ/paaAoxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase/ring 1,2‑phenylacetyl‑CoA epoxidase subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/43);  total (23/43)

AAATCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTAT  >  NC_000913_3_pae_pgi/1453684‑1453976
                                                                                                                                                  |                                                                                                                                                  
aaaTCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGTGAAATCCGTCACAGTTCAAACATACAAAGttt                                                                                                                                                 <  2:65057/150‑1 (MQ=255)
aaaTCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGttt                                                                                                                                                 <  1:333987/150‑1 (MQ=255)
 aaTCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTtg                                                                                                                                                <  4:270960/150‑1 (MQ=255)
     gcgAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGAt                                                                                                                                            >  4:200500/1‑150 (MQ=255)
        aaTCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGAtttt                                                                                                                                         <  2:812069/150‑1 (MQ=255)
         aTCATAAAAAGGAAAGGATCTTTTAACGCACTGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAAAATTGTGCGAAATCCGTCACAGTTCAAACATACAAAGTATGTGATTTTa                                                                                                                                        <  2:84337/150‑1 (MQ=255)
         aTCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTa                                                                                                                                        <  1:351558/150‑1 (MQ=255)
         aTCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTa                                                                                                                                        >  2:502017/1‑150 (MQ=255)
             tAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTa                                                                                                                                    >  4:515105/1‑150 (MQ=255)
              aaaaaGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTaa                                                                                                                                   <  1:836926/150‑1 (MQ=255)
                   gCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTAtt                                                                                                                              <  2:58893/150‑1 (MQ=255)
                   gCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTAtt                                                                                                                              <  2:226997/150‑1 (MQ=255)
                    cAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATtg                                                                                                                             <  3:50512/150‑1 (MQ=255)
                     aaaGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATtgt                                                                                                                            <  4:239111/150‑1 (MQ=255)
                       aGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTgtgt                                                                                                                          <  2:184927/150‑1 (MQ=255)
                         gATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTaa                                                                                                                        >  3:335817/1‑150 (MQ=255)
                              tttAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTc                                                                                                                   >  3:387347/1‑150 (MQ=255)
                                  aCGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTa                                                                                                                                        <  1:666369/125‑1 (MQ=255)
                                  aCGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTa                                                                                                                                        >  2:666369/1‑125 (MQ=255)
                                  aCGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATaa                                                                                                               >  4:405614/1‑150 (MQ=255)
                                   cGAACTGTTAAATATGCGATCAGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATATTACTTCACTATTGTGTAACTTTCATaaa                                                                                                              >  4:486370/1‑150 (MQ=255)
                                   cGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATaaa                                                                                                              >  1:740036/1‑150 (MQ=255)
                                    gAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATaaaa                                                                                                             >  3:327611/1‑150 (MQ=255)
                                     aaaTGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAAc                                                                                                            >  1:555453/1‑150 (MQ=255)
                                      aaTGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACa                                                                                                           <  1:454444/150‑1 (MQ=255)
                                        tGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAt                                                                                                         <  3:506111/150‑1 (MQ=255)
                                         gTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAtg                                                                                                        >  2:17870/1‑150 (MQ=255)
                                         gTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAtg                                                                                                        <  3:161780/150‑1 (MQ=255)
                                          ttAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAtgt                                                                                                       >  2:652552/1‑150 (MQ=255)
                                            aaCTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGa                                                                                                     <  4:339121/150‑1 (MQ=255)
                                                       cTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGttttt                                                                                          <  1:488374/150‑1 (MQ=255)
                                                        tGTATAGCAACTGCGGAAAACATTAATGCACTGCTAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGtttttt                                                                                         >  4:521537/1‑149 (MQ=255)
                                                         gTATAGCAACTGCGGGAAACATTAAGTCACTGATAAATAATGATTTATAAAAAGAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGAGTTTACGTAACTATTGTGTAACTTTCATAAAACAATGTGAGTCGTGTTTttaa                                                                                        <  2:695881/150‑1 (MQ=255)
                                                          tataGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTttaat                                                                                       <  2:211311/150‑1 (MQ=255)
                                                          tataGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTttaat                                                                                       >  4:320666/1‑150 (MQ=255)
                                                          tataCCAACTGAGGAAAACATTAATGCAATGATAAATGATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTttaat                                                                                       <  3:194608/150‑1 (MQ=255)
                                                            taGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTttaatta                                                                                     <  2:758640/150‑1 (MQ=255)
                                                             aGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTttaattaa                                                                                    <  4:537730/150‑1 (MQ=255)
                                                               cAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTaattaatt                                                                                  <  2:582714/150‑1 (MQ=255)
                                                               cAACTGCGAAAAACATTAATGCACTTATAAATAATTATTTATAAAAATAGGGTACTAAATCCTTCACATTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTaattaatt                                                                                  <  1:338913/150‑1 (MQ=255)
                                                                   tGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACg                                                                              >  2:75167/1‑150 (MQ=255)
                                                                      ggAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGaaa                                                                           >  1:766749/1‑150 (MQ=255)
                                                                       gAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGaaaa                                                                          >  3:58749/1‑150 (MQ=255)
                                                                          aaCATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGATATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTg                                                                       <  1:232915/150‑1 (MQ=255)
                                                                             aTTAATGCACTGATAAATACTGATTTATAAACATAGGGTGCGTGTTTAGTCCCAGTACGAACATGCTAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGaa                                                                    <  1:308288/150‑1 (MQ=255)
                                                                             aTTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAt                                                                                                         <  4:530042/113‑1 (MQ=255)
                                                                             aTTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAt                                                                                                         >  3:530042/1‑113 (MQ=255)
                                                                                       tGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGt                                                          <  2:145214/150‑1 (MQ=255)
                                                                                          gagataatGATTTATAAAAATAGGGTGCGAATTCCGTCACAGTTCAAACATAAAAAGTTTGTGATTTTAGTTAACTATTGTGTAACTTTCATAAAACAATGTGAGTCGTGTTTTTAATTAATTCTCGAAAACTGGAATCGTAAAGTtgat                                                       <  2:627752/147‑1 (MQ=255)
                                                                                                  gATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGAc                                               >  1:726264/1‑150 (MQ=255)
                                                                                                        taAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCaaga                                         <  2:600752/150‑1 (MQ=255)
                                                                                                          aaaaTAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCaagaag                                       <  3:298067/150‑1 (MQ=255)
                                                                                                           aaaTAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAagaaga                                      >  2:254663/1‑150 (MQ=255)
                                                                                                            aaTAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAAGGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAgaagaa                                     <  2:72394/150‑1 (MQ=255)
                                                                                                                 ggTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCtt                                <  3:97205/150‑1 (MQ=255)
                                                                                                                 ggTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCtt                                <  2:487577/150‑1 (MQ=255)
                                                                                                                  gTGCGAAATCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCttt                               <  3:320666/150‑1 (MQ=255)
                                                                                                                          tCCGTCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACg                       <  1:167294/150‑1 (MQ=255)
                                                                                                                              tCACAGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACACAAGAAGAACGCTTTGAGCAACGGATa                   <  1:502017/150‑1 (MQ=255)
                                                                                                                                 caGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGcc                <  3:405614/150‑1 (MQ=255)
                                                                                                                                  aGTTCAAACATACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGccc               <  2:766749/150‑1 (MQ=255)
                                                                                                                                           aTACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACgg      <  4:327611/150‑1 (MQ=255)
                                                                                                                                           aTACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACgg      <  1:17870/150‑1 (MQ=255)
                                                                                                                                            tACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAg                  <  4:250720/137‑1 (MQ=255)
                                                                                                                                            tACAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAg                  >  3:250720/1‑137 (MQ=255)
                                                                                                                                             aCAAAGTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCt    >  1:388121/1‑150 (MQ=255)
                                                                                                                                               aaagTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTACTTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGGGACCCAAGAAGAACGCTTTGAGCAACGGAAAGAACAGGAAACGGCTAt  >  2:458388/1‑150 (MQ=255)
                                                                                                                                               aaagTTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTAt  >  4:280576/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAATCGCGAATCATAAAAAGCAAAGGATCTTTTAACGAAATGTTAACTATGCGATCTGTATAGCAACTGCGGAAAACATTAATGCACTGATAAATAATGATTTATAAAAATAGGGTGCGAAATCCGTCACAGTTCAAACATACAAAATTTGTGATTTTACTTAACTATTGTGTAACTTTCATAAAACAATGTGATTCGTGTTTTTAATTAATTCACGAAAACTGGAATCGTAAAGGTGATGACGTGACCCAAGAAGAACGCTTTGAGCAACGGATAGCCCAGGAAACGGCTAT  >  NC_000913_3_pae_pgi/1453684‑1453976

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