Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,454,724 G→T 100% 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,7240GT99.5% 635.3 / ‑4.3 200G472C (GGT→TGT) paaH3‑hydroxyadipyl‑CoA dehydrogenase, NAD+‑dependent
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (90/109);  total (91/109)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 200 total aligned reads displayed.

AATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGGT  >  CP009273/1454580‑1454868
                                                                                                                                                |                                                                                                                                                
aaTTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGtt                                                                                                                                               <  2:3132647/148‑1 (MQ=255)
 aTTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTa                                                                                                                                              >  2:1823396/1‑148 (MQ=255)
   tATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATg                                                                                                                                            <  2:940782/148‑1 (MQ=255)
    aTCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGa                                                                                                                                           >  1:763681/1‑148 (MQ=255)
      ccATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGt                                                                                                                                         >  2:957500/1‑148 (MQ=255)
      ccATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGt                                                                                                                                         <  1:2995765/148‑1 (MQ=255)
      ccATATGGCCCACTTGCCTGGGGAGCGCAACTAGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGt                                                                                                                                         >  1:2023176/1‑148 (MQ=255)
       cATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTc                                                                                                                                        <  2:311076/148‑1 (MQ=255)
            ggCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAg                                                                                                                                   >  1:665640/1‑148 (MQ=255)
             gCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAgg                                                                                                                                  <  1:2862465/148‑1 (MQ=255)
                 aCTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTg                                                                                                                              >  2:2499803/1‑148 (MQ=255)
                  cTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGCGAGCTGTTATGAGTCATAAGGCCTgg                                                                                                                             >  1:98184/1‑148 (MQ=255)
                       cTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaaa                                                                                                                        <  2:243986/148‑1 (MQ=255)
                        tGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAAt                                                                                                                       >  1:2796337/1‑148 (MQ=255)
                         ggggAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATg                                                                                                                      <  1:164195/148‑1 (MQ=255)
                         ggggAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATg                                                                                                                      <  2:2888046/148‑1 (MQ=255)
                         ggggAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATg                                                                                                                      <  1:1125213/148‑1 (MQ=255)
                           ggAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGcc                                                                                                                    <  2:3218143/148‑1 (MQ=255)
                            gAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTAATAAGGCCTGGCAAAATGccc                                                                                                                   >  1:2636264/1‑148 (MQ=255)
                                 cAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCa                                                                                                              <  1:678860/148‑1 (MQ=255)
                                 cAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCa                                                                                                              >  1:2737083/1‑148 (MQ=255)
                                   aCTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAAt                                                                                                            >  1:2298175/1‑148 (MQ=255)
                                    cTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATg                                                                                                           <  1:516787/148‑1 (MQ=255)
                                      tggctggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTa                                                                                                         <  1:330869/148‑1 (MQ=255)
                                      tggctggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTa                                                                                                         <  1:631447/148‑1 (MQ=255)
                                       ggctggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTAt                                                                                                        >  2:1632816/1‑148 (MQ=255)
                                           ggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgaga                                                                                                    <  2:2412217/148‑1 (MQ=255)
                                           ggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgaga                                                                                                    <  1:1183739/148‑1 (MQ=255)
                                              aGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACg                                                                                                 <  1:2975218/148‑1 (MQ=255)
                                              aGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACg                                                                                                 <  1:2978415/148‑1 (MQ=255)
                                               gCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGa                                                                                                <  2:2418809/148‑1 (MQ=255)
                                                cGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGAt                                                                                               >  2:1557393/1‑148 (MQ=255)
                                                  aaTATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGc                                                                                             <  2:1017511/148‑1 (MQ=255)
                                                   atatTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGcc                                                                                            <  2:2695022/148‑1 (MQ=255)
                                                     atTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGATAACGATGCCTg                                                                                          >  2:2060457/1‑148 (MQ=255)
                                                       tAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTgcg                                                                                        >  2:3257193/1‑148 (MQ=255)
                                                        aaGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTgcgc                                                                                       <  2:1110193/148‑1 (MQ=255)
                                                          ggCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCa                                                                                     <  2:1075745/148‑1 (MQ=255)
                                                            cTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCaaa                                                                                   >  2:2155077/1‑148 (MQ=255)
                                                             tCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCaacg                                                                                  >  2:2222907/1‑146 (MQ=255)
                                                              ccTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgc                                                                                 >  1:1031623/1‑148 (MQ=255)
                                                              ccTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgc                                                                                 >  2:3376417/1‑148 (MQ=255)
                                                              ccTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgc                                                                                 >  2:2327722/1‑148 (MQ=255)
                                                              ccTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgc                                                                                 >  1:2153636/1‑148 (MQ=255)
                                                                ttGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgcgc                                                                               <  2:3354534/148‑1 (MQ=255)
                                                                 tGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCt                                                                              <  1:1079054/148‑1 (MQ=255)
                                                                 tGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCt                                                                              <  2:1226861/148‑1 (MQ=255)
                                                                    aaaTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTgg                                                                           >  1:2006695/1‑148 (MQ=255)
                                                                    aaaTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTgg                                                                           >  2:1291639/1‑148 (MQ=255)
                                                                       tCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCAt                                                                        <  1:1168186/148‑1 (MQ=255)
                                                                        cTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATc                                                                       >  1:1881834/1‑148 (MQ=255)
                                                                        cTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATc                                                                       >  2:1379468/1‑148 (MQ=255)
                                                                        cTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCGGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCCTGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATc                                                                       >  1:3384805/1‑148 (MQ=255)
                                                                            aGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACa                                                                   <  1:1789436/148‑1 (MQ=255)
                                                                              cATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAtt                                                                 >  2:78647/1‑148 (MQ=255)
                                                                              cATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAtt                                                                 <  1:2392129/148‑1 (MQ=255)
                                                                              cATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAtt                                                                 <  2:1774158/148‑1 (MQ=255)
                                                                              cATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAtt                                                                 <  1:219738/148‑1 (MQ=255)
                                                                              cATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAtt                                                                 <  2:1939047/148‑1 (MQ=255)
                                                                               aTCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTa                                                                <  2:3069135/148‑1 (MQ=255)
                                                                               aTCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTa                                                                <  2:2780009/148‑1 (MQ=255)
                                                                               aTCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTa                                                                <  2:2636264/148‑1 (MQ=255)
                                                                                   cTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTAtctc                                                            >  2:1949602/1‑148 (MQ=255)
                                                                                    tATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                           >  1:468208/1‑148 (MQ=255)
                                                                                    tATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                           >  2:797029/1‑148 (MQ=255)
                                                                                    tATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCa                                                           >  1:3063913/1‑148 (MQ=255)
                                                                                         cGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGa                                                      <  2:350035/148‑1 (MQ=255)
                                                                                         cGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGa                                                      <  1:1675889/148‑1 (MQ=255)
                                                                                          gaagaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGAt                                                     <  1:897823/148‑1 (MQ=255)
                                                                                           aagaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATg                                                    >  2:2994942/1‑148 (MQ=255)
                                                                                           aagaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAGTATCTCAATGGATg                                                    >  2:1980173/1‑148 (MQ=255)
                                                                                            agaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGa                                                   <  1:950121/148‑1 (MQ=255)
                                                                                            agaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGa                                                   <  1:3071067/148‑1 (MQ=255)
                                                                                            agaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGa                                                   <  2:3021676/148‑1 (MQ=255)
                                                                                               aCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAgg                                                <  2:2796337/148‑1 (MQ=255)
                                                                                                cGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAAAGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATATCAATGGATGAAGGc                                               >  1:2124155/1‑148 (MQ=255)
                                                                                                cGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAAGGGATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                               >  2:100526/1‑148 (MQ=255)
                                                                                                  cTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCtt                                             <  1:1557393/148‑1 (MQ=255)
                                                                                                    aTCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTg                                           <  1:2867141/148‑1 (MQ=255)
                                                                                                        cccATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt                                       >  1:1013799/1‑148 (MQ=255)
                                                                                                        cccATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt                                       <  2:3036856/148‑1 (MQ=255)
                                                                                                         ccATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTa                                      <  1:1496176/148‑1 (MQ=255)
                                                                                                         ccATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTa                                      >  1:2440385/1‑148 (MQ=255)
                                                                                                         ccATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTAATAAGGCCTGGCAAAATTCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTa                                      >  1:653843/1‑148 (MQ=255)
                                                                                                            tGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTg                                   <  2:260572/148‑1 (MQ=255)
                                                                                                            tGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTg                                   >  2:2838490/1‑148 (MQ=255)
                                                                                                             gTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGa                                  >  2:3469395/1‑148 (MQ=255)
                                                                                                             gTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGa                                  >  2:547752/1‑148 (MQ=255)
                                                                                                                cATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                               <  2:2389761/148‑1 (MQ=255)
                                                                                                                  ttGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGa                                  <  1:1820168/143‑1 (MQ=255)
                                                                                                                  ttGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATg                             >  2:1440204/1‑148 (MQ=255)
                                                                                                                            aaCGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTg                   <  1:2499803/148‑1 (MQ=255)
                                                                                                                              cGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGca                 <  2:2923893/148‑1 (MQ=255)
                                                                                                                              cGGGCGCGTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGca                 <  2:729792/148‑1 (MQ=255)
                                                                                                                                  cgcTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa             <  1:3148908/148‑1 (MQ=255)
                                                                                                                                      tCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATTCt         >  2:1924120/1‑148 (MQ=255)
                                                                                                                                      tCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCt         <  1:3376417/148‑1 (MQ=255)
                                                                                                                                         ggAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTaa      >  1:85317/1‑148 (MQ=255)
                                                                                                                                           agagCTGTTATGAGTCATAAGGCCTGACAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACg    >  2:1589468/1‑148 (MQ=255)
                                                                                                                                             agCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGGt  >  1:1461034/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGGT  >  CP009273/1454580‑1454868

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: ‑