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% 728.7 / ‑4.3 240G472C (GGT→TGT) paaH3‑hydroxyadipyl‑CoA dehydrogenase, NAD+‑dependent
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (115/123);  total (117/123)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 240 total aligned reads displayed.

AATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGG  >  CP009273/1454580‑1454867
                                                                                                                                                |                                                                                                                                               
 aTTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTa                                                                                                                                             >  1:1024786/1‑148 (MQ=255)
   tATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATg                                                                                                                                           >  1:3119427/1‑148 (MQ=255)
      ccATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGt                                                                                                                                        >  1:963051/1‑148 (MQ=255)
      ccATATCGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGt                                                                                                                                        <  1:12088/148‑1 (MQ=255)
       cATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTc                                                                                                                                       <  2:878568/148‑1 (MQ=255)
         tatGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCAt                                                                                                                                     >  2:1966992/1‑148 (MQ=255)
             gCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAgg                                                                                                                                 >  2:2142783/1‑148 (MQ=255)
              cccACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGc                                                                                                                                <  1:1932020/148‑1 (MQ=255)
              cccACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGc                                                                                                                                <  2:3351104/148‑1 (MQ=255)
               ccACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGcc                                                                                                                               >  1:2713172/1‑148 (MQ=255)
               ccACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGcc                                                                                                                               >  1:2712305/1‑148 (MQ=255)
               ccACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGcc                                                                                                                               <  1:689581/148‑1 (MQ=255)
                cACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCt                                                                                                                              <  2:1122679/148‑1 (MQ=255)
                 aCTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTg                                                                                                                             >  2:2247028/1‑148 (MQ=255)
                  cTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGATTCATAATGCCTgg                                                                                                                            >  1:3226796/1‑148 (MQ=255)
                    tGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCa                                                                                                                          >  1:1337813/1‑148 (MQ=255)
                     gCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaa                                                                                                                         <  2:449929/148‑1 (MQ=255)
                      ccTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaa                                                                                                                        >  2:2848499/1‑148 (MQ=255)
                      ccTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaa                                                                                                                        <  2:2032228/148‑1 (MQ=255)
                      ccTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaa                                                                                                                        <  1:243735/148‑1 (MQ=255)
                      ccTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaa                                                                                                                        >  2:209336/1‑148 (MQ=255)
                       cTGGGGAGCGCAACTTGTCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaaa                                                                                                                       <  2:260274/148‑1 (MQ=255)
                       cTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCaaaa                                                                                                                       <  2:2346657/148‑1 (MQ=255)
                         ggggAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATg                                                                                                                     <  1:2466565/148‑1 (MQ=255)
                           ggAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGcc                                                                                                                   >  1:2126723/1‑148 (MQ=255)
                           ggAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGcc                                                                                                                   >  1:332295/1‑148 (MQ=255)
                            gAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGccc                                                                                                                  <  2:2162910/148‑1 (MQ=255)
                            gAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGccc                                                                                                                  >  2:1970990/1‑148 (MQ=255)
                            gAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGccc                                                                                                                  >  2:2885468/1‑148 (MQ=255)
                                 cAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCa                                                                                                             <  2:3399386/148‑1 (MQ=255)
                                 cAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCa                                                                                                             <  2:3485026/148‑1 (MQ=255)
                                 cAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCa                                                                                                             <  1:1329067/148‑1 (MQ=255)
                                    cTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATg                                                                                                          >  2:3031838/1‑148 (MQ=255)
                                     ttGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGt                                                                                                         <  2:3179985/148‑1 (MQ=255)
                                      tggctggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTa                                                                                                        <  1:177195/148‑1 (MQ=255)
                                        gctggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATg                                                                                                      <  2:978604/148‑1 (MQ=255)
                                          tggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgag                                                                                                    >  2:1535557/1‑148 (MQ=255)
                                           ggcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgaga                                                                                                   <  2:2490046/148‑1 (MQ=255)
                                            gcAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGaa                                                                                                  <  2:2684093/148‑1 (MQ=255)
                                               gCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGa                                                                                               <  2:905869/148‑1 (MQ=255)
                                               gCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGa                                                                                               <  2:1653602/148‑1 (MQ=255)
                                                cGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGAt                                                                                              >  2:1387433/1‑148 (MQ=255)
                                                cGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGAt                                                                                              >  2:130076/1‑148 (MQ=255)
                                                 gAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATg                                                                                             <  1:217020/148‑1 (MQ=255)
                                                      ttAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTgc                                                                                        >  1:1478975/1‑148 (MQ=255)
                                                       tAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGAGCCAACGGGCGCTTCTGGATAGCTGTTATGAGTCATAACGCCTGGCAAAATGCCCATGCAATGTATTAGAACGAAGCCTgcg                                                                                       >  1:1991045/1‑148 (MQ=255)
                                                          ggCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCa                                                                                    >  2:10693/1‑148 (MQ=255)
                                                          ggCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCa                                                                                    >  2:10951/1‑148 (MQ=255)
                                                           gCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCaa                                                                                   >  2:1160828/1‑148 (MQ=255)
                                                               cTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgcg                                                                               <  1:1966175/148‑1 (MQ=255)
                                                                ttGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgcgc                                                                              <  2:2016476/148‑1 (MQ=255)
                                                                ttGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgcgc                                                                              >  1:2199731/1‑148 (MQ=255)
                                                                 tGAAAATCTACAGCTTCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCt                                                                             >  1:3086823/1‑148 (MQ=255)
                                                                 tGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATTCAATGTATGAGAACGATGCCTGCGCCAAAGCGCt                                                                             <  2:686066/148‑1 (MQ=255)
                                                                 tGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCt                                                                             >  1:1701343/1‑148 (MQ=255)
                                                                      aTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCa                                                                        >  1:2219217/1‑148 (MQ=255)
                                                                      aTCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCGTGGCa                                                                        >  1:408461/1‑148 (MQ=255)
                                                                        cTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATc                                                                      >  2:152108/1‑148 (MQ=255)
                                                                               aTCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTa                                                               <  1:3320266/148‑1 (MQ=255)
                                                                                tCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTAt                                                              >  1:369880/1‑148 (MQ=255)
                                                                                     aTGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCaa                                                         <  2:1024786/148‑1 (MQ=255)
                                                                                     aTGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCaa                                                         <  1:1274194/148‑1 (MQ=255)
                                                                                        gCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATgg                                                      <  1:1159183/148‑1 (MQ=255)
                                                                                           aagaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCAAGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAATATCTCAATGGATg                                                   >  1:3291558/1‑148 (MQ=255)
                                                                                             gaaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGACCGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGaa                                                 <  2:2664962/148‑1 (MQ=255)
                                                                                              aaCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAg                                                >  2:2527079/1‑148 (MQ=255)
                                                                                                cGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGc                                              >  1:2639290/1‑148 (MQ=255)
                                                                                                   tATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATTCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATATCAATGGATGAAGGCttt                                           >  1:18221/1‑148 (MQ=255)
                                                                                                   tATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCttt                                           <  1:2992491/148‑1 (MQ=255)
                                                                                                     tCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGc                                         >  2:409580/1‑148 (MQ=255)
                                                                                                      cGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCt                                        >  1:857264/1‑148 (MQ=255)
                                                                                                         ccATGTTCATTGCTGCGCCAACGGGCGCTTCGGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCAAGCAATGTATTAGAACGATGCCTGCGCCAAATCGCTTGGCATCGACATTATATCAATGGATGAAGGCTTTGCTGTa                                     >  1:1616732/1‑148 (MQ=255)
                                                                                                           aTGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGt                                   <  2:2582941/148‑1 (MQ=255)
                                                                                                           aTGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGt                                   <  2:3137207/148‑1 (MQ=255)
                                                                                                               tCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGAcc                               <  1:2629497/148‑1 (MQ=255)
                                                                                                               tCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGAcc                               <  2:2944182/148‑1 (MQ=255)
                                                                                                                cATTGCTGCGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCa                              <  2:3133025/148‑1 (MQ=255)
                                                                                                                   tgctgcGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGa                           <  1:464000/148‑1 (MQ=255)
                                                                                                                   tgctgcGCCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGa                           <  2:1744416/148‑1 (MQ=255)
                                                                                                                       gcgcCAACGGGCGCTTGTGGAGAGCTGTTATGAGTAATAATGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCAGCGCCAAAGCGCTTGGCATCGACCATATCACAAATGATGAAGGCTTTGCTGTAGTGACCATGACCGt                       >  1:3309130/1‑148 (MQ=255)
                                                                                                                       gcgcCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGt                       <  1:2243998/148‑1 (MQ=255)
                                                                                                                       gcgcCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAAGGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGt                       <  1:3259193/148‑1 (MQ=255)
                                                                                                                       gcgcCAACGGGCGCTTCTGGAGAGATGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGACCGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGt                       <  1:3387966/148‑1 (MQ=255)
                                                                                                                         gcCAACGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCa                     <  1:2247028/148‑1 (MQ=255)
                                                                                                                            aaCGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTg                  <  1:2268560/148‑1 (MQ=255)
                                                                                                                              cGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTCGTGACCATGACCGTCACTGca                >  2:2180644/1‑148 (MQ=255)
                                                                                                                              cGGGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGca                >  2:2259096/1‑148 (MQ=255)
                                                                                                                               gggCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGcac               >  2:273798/1‑148 (MQ=255)
                                                                                                                                ggCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGcaca              >  2:2430360/1‑148 (MQ=255)
                                                                                                                                 gcgcTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAAGGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaa             <  2:8186/148‑1 (MQ=255)
                                                                                                                                  cgcTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACaaa            <  2:3533832/148‑1 (MQ=255)
                                                                                                                                   gcTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAAt           <  2:1227825/148‑1 (MQ=255)
                                                                                                                                     ttCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGc         <  2:1478975/148‑1 (MQ=255)
                                                                                                                                      tCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCt        >  2:2360260/1‑148 (MQ=255)
                                                                                                                                      tCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCt        >  1:2223040/1‑148 (MQ=255)
                                                                                                                                       cTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCtt       <  2:2785550/148‑1 (MQ=255)
                                                                                                                                        tGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTa      <  2:963051/148‑1 (MQ=255)
                                                                                                                                         ggAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTaa     >  1:754882/1‑148 (MQ=255)
                                                                                                                                         ggAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTaa     >  1:2571679/1‑148 (MQ=255)
                                                                                                                                          gagaGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAAc    >  1:2203217/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGACCGTCACTGCACAAATGCTTAACGG  >  CP009273/1454580‑1454867

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