Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 495,020 T→A 100% E137V (GAA→GTA)  aes ← acetyl esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273495,0200TA97.7% 696.2 / 4.3 216E137V (GAA→GTA) aesacetyl esterase
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base A (108/103);  total (109/107)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
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 216 total aligned reads displayed.

GCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAAT  >  CP009273/494875‑495164
                                                                                                                                                 |                                                                                                                                                
gCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTAcc                                                                                                                                                <  1:457442/148‑1 (MQ=255)
 cAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACct                                                                                                                                               <  1:635664/148‑1 (MQ=255)
 cAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACct                                                                                                                                               <  1:3142283/148‑1 (MQ=255)
  aaTCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACctc                                                                                                                                              >  2:448012/1‑148 (MQ=255)
    tCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTa                                                                                                                                            >  1:920134/1‑148 (MQ=255)
    tCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTa                                                                                                                                            <  1:202951/148‑1 (MQ=255)
     cGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTAt                                                                                                                                           <  2:2512493/148‑1 (MQ=255)
       aTCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCg                                                                                                                                         >  1:3345165/1‑148 (MQ=255)
        tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        <  2:1616036/148‑1 (MQ=255)
        tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        >  1:2530289/1‑148 (MQ=255)
         cTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCt                                                                                                                                       >  2:30025/1‑148 (MQ=255)
         cTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCt                                                                                                                                       >  2:1850915/1‑148 (MQ=255)
           gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTg                                                                                                                                     <  2:2188403/148‑1 (MQ=255)
           gTTGATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTg                                                                                                                                     <  2:1016925/148‑1 (MQ=255)
            tttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGc                                                                                                                                    >  1:1250162/1‑148 (MQ=255)
             ttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCg                                                                                                                                   <  1:1179870/148‑1 (MQ=255)
                tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaa                                                                                                                                >  1:2935722/1‑148 (MQ=255)
                tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaa                                                                                                                                >  1:2935429/1‑148 (MQ=255)
                 cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCCATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaa                                                                                                                               >  2:190272/1‑148 (MQ=255)
                 cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaa                                                                                                                               >  2:190266/1‑148 (MQ=255)
                  aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaaa                                                                                                                              <  1:49443/148‑1 (MQ=255)
                  aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaaa                                                                                                                              >  1:2901030/1‑148 (MQ=255)
                    gCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAAcg                                                                                                                            <  1:3242818/148‑1 (MQ=255)
                    gCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAAcg                                                                                                                            <  1:1813870/148‑1 (MQ=255)
                     cAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAAcgc                                                                                                                           <  2:849969/148‑1 (MQ=255)
                      aaCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAAcgcg                                                                                                                          <  1:3170890/148‑1 (MQ=255)
                       aCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACgcgc                                                                                                                         >  1:349094/1‑148 (MQ=255)
                        ccACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGACATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCt                                                                                                                        >  1:3451416/1‑148 (MQ=255)
                           caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     >  1:1293855/1‑148 (MQ=255)
                           caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     <  1:3310059/148‑1 (MQ=255)
                            aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAg                                                                                                                    <  1:3204066/148‑1 (MQ=255)
                             aCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGAAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAgg                                                                                                                   <  2:3375982/148‑1 (MQ=255)
                                  cTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGaaa                                                                                                              >  1:2531839/1‑148 (MQ=255)
                                    ggCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAgg                                                                                                            >  1:2667579/1‑148 (MQ=255)
                                     gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAggg                                                                                                           >  1:2147987/1‑148 (MQ=255)
                                       gAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgtg                                                                                                         <  1:133946/148‑1 (MQ=255)
                                        aGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgtgt                                                                                                        <  2:321330/148‑1 (MQ=255)
                                         gcgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAACACGCGCTTCAGGTGAAAGGGTGTa                                                                                                       >  2:2578442/1‑148 (MQ=255)
                                         gcgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTa                                                                                                       <  1:2921998/148‑1 (MQ=255)
                                               gCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAAt                                                                                                 <  1:2454520/148‑1 (MQ=255)
                                                    gCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCaa                                                                                            >  1:1579802/1‑148 (MQ=255)
                                                       ccTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCa                                                                                         <  2:1645095/148‑1 (MQ=255)
                                                         tGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCAcc                                                                                       <  2:1966984/148‑1 (MQ=255)
                                                          gCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCg                                                                                      <  2:2946028/148‑1 (MQ=255)
                                                             gAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTAc                                                                                   <  2:2890020/148‑1 (MQ=255)
                                                             gAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTAc                                                                                   >  1:253282/1‑148 (MQ=255)
                                                              aaTCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACa                                                                                  <  1:3089814/148‑1 (MQ=255)
                                                                      gCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGt                                                                          <  2:282878/148‑1 (MQ=255)
                                                                          aGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCt                                                                      >  1:1900119/1‑148 (MQ=255)
                                                                          aGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCt                                                                      >  2:557944/1‑148 (MQ=255)
                                                                          aGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCt                                                                      >  2:358971/1‑148 (MQ=255)
                                                                             cAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGc                                                                   <  2:2608918/148‑1 (MQ=255)
                                                                             cAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGc                                                                   >  1:1978834/1‑148 (MQ=255)
                                                                                    ggACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAgg                                                            >  1:1407064/1‑148 (MQ=255)
                                                                                    ggACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAgg                                                            >  1:2481323/1‑148 (MQ=255)
                                                                                      aCATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcg                                                          >  1:3550022/1‑148 (MQ=255)
                                                                                         tatTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCat                                                       >  1:737817/1‑148 (MQ=255)
                                                                                         tatTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCat                                                       >  2:3179438/1‑148 (MQ=255)
                                                                                          atTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatg                                                      <  2:3317811/148‑1 (MQ=255)
                                                                                           ttgatttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatga                                                     <  2:1821898/148‑1 (MQ=255)
                                                                                           ttgatttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatga                                                     <  2:2095966/148‑1 (MQ=255)
                                                                                               tttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCg                                                 >  2:1418520/1‑148 (MQ=255)
                                                                                                ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  1:382504/1‑148 (MQ=255)
                                                                                                ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  2:741697/1‑148 (MQ=255)
                                                                                                 tgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGAt                                               >  1:95979/1‑148 (MQ=255)
                                                                                                  gatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATc                                              <  2:1352408/148‑1 (MQ=255)
                                                                                                      aTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTgg                                          <  2:2066124/148‑1 (MQ=255)
                                                                                                      aTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTgg                                          <  1:2979629/148‑1 (MQ=255)
                                                                                                       tcctccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTggg                                         <  1:1488567/148‑1 (MQ=255)
                                                                                                        cctccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGt                                        >  2:1169999/1‑148 (MQ=255)
                                                                                                         ctccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTa                                       <  2:3083879/148‑1 (MQ=255)
                                                                                                          tccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTAt                                      <  2:1497060/148‑1 (MQ=255)
                                                                                                           ccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATc                                     >  1:2550355/1‑148 (MQ=255)
                                                                                                              ccTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgag                                  >  1:233682/1‑148 (MQ=255)
                                                                                                              ccTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgag                                  >  1:1652192/1‑148 (MQ=255)
                                                                                                              ccTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgag                                  >  1:1717389/1‑148 (MQ=255)
                                                                                                                  ctgGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTg                              <  1:2942222/148‑1 (MQ=255)
                                                                                                                     gtggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCg                           >  1:1877858/1‑148 (MQ=255)
                                                                                                                     gtggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCg                           <  1:311415/148‑1 (MQ=255)
                                                                                                                        gAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCgaga                        >  2:2809301/1‑148 (MQ=255)
                                                                                                                         aaaTAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGaa                       <  2:3541925/148‑1 (MQ=255)
                                                                                                                         aaaTAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGaa                       <  2:3494917/148‑1 (MQ=255)
                                                                                                                           aTAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                     >  2:1181951/1‑148 (MQ=255)
                                                                                                                           aTAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                     >  2:2336390/1‑148 (MQ=255)
                                                                                                                            tAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaa                    >  2:544993/1‑148 (MQ=255)
                                                                                                                             aacaacAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaaa                   >  1:2039219/1‑148 (MQ=255)
                                                                                                                                 acaaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCg               <  1:448013/148‑1 (MQ=255)
                                                                                                                                 acaaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCg               >  1:896521/1‑148 (MQ=255)
                                                                                                                                  caaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGc              >  1:2982246/1‑148 (MQ=255)
                                                                                                                                   aaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGcc             <  2:831772/148‑1 (MQ=255)
                                                                                                                                     gcagcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTc           <  2:2716751/148‑1 (MQ=255)
                                                                                                                                      cagcCACATTTACCTCTATCGCTTGCGGAAAACGCGCTTCTGGTGAAAGGGTGTAATCAATACCAATCACCGTACTTTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTcc          <  2:920134/148‑1 (MQ=255)
                                                                                                                                       agcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGACTAAAACCGCCTCCa         >  2:885902/1‑148 (MQ=255)
                                                                                                                                       agcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCa         >  2:2666842/1‑148 (MQ=255)
                                                                                                                                         ccACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATg       >  2:66679/1‑148 (MQ=255)
                                                                                                                                         ccACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATg       >  1:78032/1‑148 (MQ=255)
                                                                                                                                            caATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaa    >  2:3188214/1‑148 (MQ=255)
                                                                                                                                             aaTTACCTCTATCGCTTGCTGAAAACGCTCTTCTGGTGAAAGGGTGTAATCAATACCAATCACCGTACTTTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATGGTGGGTATCGAGTTTGCCGACAATAAAACCGCCTCCATGCaaa   <  1:613188/148‑1 (MQ=255)
                                                                                                                                              aTTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaat  <  2:2481323/148‑1 (MQ=255)
                                                                                                                                              aTTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaat  >  1:2727556/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAAT  >  CP009273/494875‑495164

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