Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,454,724 G→T 92.0% G472C (GGT→TGT)  paaH → 3‑hydroxyadipyl‑CoA dehydrogenase, NAD+‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,454,7240GT92.0% 1018.8 / 75.1 353G472C (GGT→TGT) paaH3‑hydroxyadipyl‑CoA dehydrogenase, NAD+‑dependent
Reads supporting (aligned to +/- strand):  ref base G (22/7);  new base T (156/167);  total (179/174)
Fisher's exact test for biased strand distribution p-value = 5.89e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 353 total aligned reads displayed.

GTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGGTCA  >  CP009273/1454577‑1454870
                                                                                                                                                   |                                                                                                                                                  
 tGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGtt                                                                                                                                                 <  1:4376833/150‑1 (MQ=255)
 tGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGtt                                                                                                                                                 <  2:1323797/150‑1 (MQ=255)
  gAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTa                                                                                                                                                >  2:711108/1‑150 (MQ=255)
    aTTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATg                                                                                                                                              >  2:373842/1‑150 (MQ=255)
      tATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAg                                                                                                                                            >  1:2498266/1‑150 (MQ=255)
        tCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTc                                                                                                                                          <  2:3686218/150‑1 (MQ=255)
        tCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTc                                                                                                                                          <  2:1621687/150‑1 (MQ=255)
        tCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCGTCTGGAGAGCTGTTATGAGTc                                                                                                                                          >  1:1071100/1‑150 (MQ=255)
         ccATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCa                                                                                                                                         <  2:447296/150‑1 (MQ=255)
         ccATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCa                                                                                                                                         >  2:3597872/1‑150 (MQ=255)
            tatGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATaa                                                                                                                                      <  2:4080627/150‑1 (MQ=255)
               ggCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGc                                                                                                                                   >  1:4649728/1‑150 (MQ=255)
                 cccACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCt                                                                                                                                 >  2:1480983/1‑150 (MQ=255)
                    aCTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGATCTGTTATGAGTCATAAGGCCTGGc                                                                                                                              >  1:1818654/1‑150 (MQ=255)
                    aCTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGc                                                                                                                              <  1:4181565/150‑1 (MQ=255)
                      ttGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTATGGCTCCTTGATAATATACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGATGCGCCATGGGGCGCTTATGGAGTGATGTTATGAGTCATAAGGCATGGCaa                                                                                                                            <  1:1361248/150‑1 (MQ=255)
                       tGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaa                                                                                                                           >  2:4725968/1‑150 (MQ=255)
                       tGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaa                                                                                                                           <  2:521353/150‑1 (MQ=255)
                         ccTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAAt                                                                                                                         <  1:1239357/150‑1 (MQ=255)
                           tGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGc                                                                                                                       >  1:3457614/1‑150 (MQ=255)
                           tGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGc                                                                                                                       >  2:3429043/1‑150 (MQ=255)
                                  cgcAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCa                                                                                                                <  2:3292411/150‑1 (MQ=255)
                                    cAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAAt                                                                                                              <  2:3142166/150‑1 (MQ=255)
                                        ttGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTAt                                                                                                          <  1:2672299/150‑1 (MQ=255)
                                        ttGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTAt                                                                                                          >  1:3295823/1‑150 (MQ=255)
                                           tctggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATTGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgag                                                                                                       <  1:2181055/149‑1 (MQ=255)
                                            ctggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgaga                                                                                                      >  2:1568368/1‑150 (MQ=255)
                                              ggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAAc                                                                                                    <  1:4534829/150‑1 (MQ=255)
                                                 aGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGAt                                                                                                 <  1:1661914/150‑1 (MQ=255)
                                                    gAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGcc                                                                                              >  1:1551549/1‑150 (MQ=255)
                                                     aaTATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCt                                                                                             <  2:4614987/150‑1 (MQ=255)
                                                      atatTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCACTGTATGAGAACGATGCCTg                                                                                            >  2:1288316/1‑150 (MQ=255)
                                                       tatTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTgc                                                                                           >  2:1205222/1‑150 (MQ=255)
                                                       tatTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTgc                                                                                           >  1:4082148/1‑150 (MQ=255)
                                                         ttAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTgcgc                                                                                         >  1:3921000/1‑150 (MQ=255)
                                                          tAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGcc                                                                                        <  1:130987/150‑1 (MQ=255)
                                                            aGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCaa                                                                                      <  1:2399540/150‑1 (MQ=255)
                                                             ggCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCaaa                                                                                     >  2:1863840/1‑150 (MQ=255)
                                                             ggCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCaaa                                                                                     >  2:420111/1‑150 (MQ=255)
                                                              gCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAg                                                                                    >  2:943511/1‑150 (MQ=255)
                                                              gCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAg                                                                                    >  1:1900509/1‑150 (MQ=255)
                                                               cTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAATGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgc                                                                                   >  1:646566/1‑150 (MQ=255)
                                                                 ccTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgcgc                                                                                 <  1:37841/150‑1 (MQ=255)
                                                                    tGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCGTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCCAGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTg                                                                              >  1:543700/1‑150 (MQ=255)
                                                                      aaaaTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGc                                                                            <  2:2339705/150‑1 (MQ=255)
                                                                        aaTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCAt                                                                          >  1:2649560/1‑150 (MQ=255)
                                                                        aaTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCAt                                                                          >  1:940552/1‑150 (MQ=255)
                                                                         aTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATc                                                                         <  2:195424/150‑1 (MQ=255)
                                                                         aTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCGTGGCCTc                                                                         >  1:3229111/1‑150 (MQ=255)
                                                                           cTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGa                                                                       >  1:102888/1‑150 (MQ=255)
                                                                             aCAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACa                                                                     >  2:2294994/1‑150 (MQ=255)
                                                                                   tCAATATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGATAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTAtct                                                               <  1:4489588/150‑1 (MQ=255)
                                                                                      cTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCaa                                                            <  2:4530487/150‑1 (MQ=255)
                                                                                          ggCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGa                                                        <  1:816702/150‑1 (MQ=255)
                                                                                           gCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGAt                                                       <  2:1135395/150‑1 (MQ=255)
                                                                                            cGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATg                                                      <  2:1413533/150‑1 (MQ=255)
                                                                                              aagaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGaa                                                    >  2:2700400/1‑150 (MQ=255)
                                                                                              aagaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGaa                                                    <  1:971616/150‑1 (MQ=255)
                                                                                              aagaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGaa                                                    >  1:3349963/1‑150 (MQ=255)
                                                                                               agaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAg                                                   >  2:865892/1‑150 (MQ=255)
                                                                                               agaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAg                                                   <  1:2777525/150‑1 (MQ=255)
                                                                                               agaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAg                                                   >  2:490318/1‑150 (MQ=255)
                                                                                                 aaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                                 >  1:1963937/1‑150 (MQ=255)
                                                                                                  aCGCTATCGCCCATGTTCAGTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCt                                                <  2:3838801/150‑1 (MQ=255)
                                                                                                    gCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCttt                                              >  1:2551696/1‑150 (MQ=255)
                                                                                                      tATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGc                                            <  1:3319683/150‑1 (MQ=255)
                                                                                                       aTCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCt                                           >  1:1925487/1‑150 (MQ=255)
                                                                                                       aTCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCt                                           >  1:3529244/1‑150 (MQ=255)
                                                                                                         cGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt                                         <  2:2686012/150‑1 (MQ=255)
                                                                                                          gCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTa                                        >  2:2017908/1‑150 (MQ=255)
                                                                                                          gCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTa                                        >  2:2057364/1‑150 (MQ=255)
                                                                                                           cccATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                                       <  1:1751654/150‑1 (MQ=255)
                                                                                                           cccATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAg                                       >  2:2535751/1‑150 (MQ=255)
                                                                                                             cATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTg                                     <  1:375114/150‑1 (MQ=255)
                                                                                                              aTGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGa                                    <  1:772275/150‑1 (MQ=255)
                                                                                                               tGTTCATTGCTGCGCCAACGAGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGAc                                   >  2:3035905/1‑150 (MQ=255)
                                                                                                                gTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGAcc                                  >  1:4124790/1‑150 (MQ=255)
                                                                                                                 ttCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                                 <  1:909168/150‑1 (MQ=255)
                                                                                                                 ttCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                                 >  2:1690352/1‑150 (MQ=255)
                                                                                                                 ttCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                                 >  2:4707875/1‑150 (MQ=255)
                                                                                                                   cATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                               >  1:3069397/1‑150 (MQ=255)
                                                                                                                   cATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                               >  2:2471767/1‑150 (MQ=255)
                                                                                                                   cATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                               >  2:4074045/1‑150 (MQ=255)
                                                                                                                          gcgcCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCa                        <  1:881392/150‑1 (MQ=255)
                                                                                                                            gcCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACt                      <  2:189046/150‑1 (MQ=255)
                                                                                                                              cAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGc                    >  1:3894844/1‑150 (MQ=255)
                                                                                                                               aaCGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGca                   <  2:4504920/150‑1 (MQ=255)
                                                                                                                               aaCGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGca                   <  2:1068621/150‑1 (MQ=255)
                                                                                                                                  gggCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaa                <  1:2102052/150‑1 (MQ=255)
                                                                                                                                     cgcTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATg             >  2:2886127/1‑150 (MQ=255)
                                                                                                                                     cgcTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATg             <  2:1780508/150‑1 (MQ=255)
                                                                                                                                      gcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGc            >  2:4423425/1‑150 (MQ=255)
                                                                                                                                      gcTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGc            >  2:3811519/1‑150 (MQ=255)
                                                                                                                                       cTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCt           >  2:1015318/1‑150 (MQ=255)
                                                                                                                                          cTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTaa        >  1:1989269/1‑150 (MQ=255)
                                                                                                                                          cTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTaa        >  2:3546476/1‑150 (MQ=255)
                                                                                                                                            ggAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACg      >  2:275315/1‑150 (MQ=255)
                                                                                                                                            ggAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACg      >  1:3829581/1‑150 (MQ=255)
                                                                                                                                              agagCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGGt    >  2:4352502/1‑150 (MQ=255)
                                                                                                                                               gagCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGGtc   <  1:1574808/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGGTCA  >  CP009273/1454577‑1454870

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