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,0200TA94.5% 714.8 / 24.8 234E137V (GAA→GTA) aesacetyl esterase
Reads supporting (aligned to +/- strand):  ref base T (1/12);  new base A (99/122);  total (100/134)
Fisher's exact test for biased strand distribution p-value = 8.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 234 total aligned reads displayed.

CGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAAT  >  CP009273/494874‑495164
                                                                                                                                                  |                                                                                                                                                
 gCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTAcc                                                                                                                                                <  1:2741915/148‑1 (MQ=255)
   aaTCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACctc                                                                                                                                              <  2:2904018/148‑1 (MQ=255)
      cGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTAt                                                                                                                                           >  1:682588/1‑148 (MQ=255)
        aTCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGGAATAACAACAAGCAGCCACAATTACCTCTATCg                                                                                                                                         >  1:1389786/1‑148 (MQ=255)
        aTCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCg                                                                                                                                         >  1:3660451/1‑148 (MQ=255)
         tCTGTTTATCACGCAACCACGACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        >  2:641817/1‑148 (MQ=255)
         tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        >  1:3431372/1‑148 (MQ=255)
         tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        >  1:3051665/1‑148 (MQ=255)
         tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        <  1:1646699/148‑1 (MQ=255)
         tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        >  1:3361889/1‑148 (MQ=255)
         tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        >  1:2556023/1‑148 (MQ=255)
          cTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCt                                                                                                                                       <  1:954464/148‑1 (MQ=255)
          cTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCt                                                                                                                                       >  1:2509236/1‑148 (MQ=255)
          cTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCt                                                                                                                                       >  2:611912/1‑148 (MQ=255)
            gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTg                                                                                                                                     <  2:686448/148‑1 (MQ=255)
            gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTg                                                                                                                                     >  2:1694628/1‑148 (MQ=255)
             tttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATGTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGc                                                                                                                                    >  1:3562406/1‑148 (MQ=255)
                 tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaa                                                                                                                                >  2:398610/1‑148 (MQ=255)
                  cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaa                                                                                                                               <  2:2214345/148‑1 (MQ=255)
                  cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaa                                                                                                                               <  1:3225772/148‑1 (MQ=255)
                  cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaa                                                                                                                               <  1:851467/148‑1 (MQ=255)
                   aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCATCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaaa                                                                                                                              <  2:261398/148‑1 (MQ=255)
                   aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaaa                                                                                                                              <  1:769029/148‑1 (MQ=255)
                    cGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAAc                                                                                                                             >  2:1109574/1‑148 (MQ=255)
                      cAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAAcgc                                                                                                                           >  2:165751/1‑148 (MQ=255)
                        aCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACgcgc                                                                                                                         >  1:3491328/1‑148 (MQ=255)
                        aCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACgcgc                                                                                                                         >  1:1340111/1‑148 (MQ=255)
                          cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCtt                                                                                                                       <  1:2404781/148‑1 (MQ=255)
                          cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCtt                                                                                                                       <  1:3012758/148‑1 (MQ=255)
                          cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCtt                                                                                                                       <  2:3679593/148‑1 (MQ=255)
                          cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCtt                                                                                                                       <  2:583924/148‑1 (MQ=255)
                          cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCtt                                                                                                                       >  2:3602000/1‑148 (MQ=255)
                            caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACACGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     <  1:3137675/148‑1 (MQ=255)
                            caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     <  2:3297221/148‑1 (MQ=255)
                            caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     <  1:2296785/148‑1 (MQ=255)
                             aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAg                                                                                                                    >  2:296644/1‑148 (MQ=255)
                               cGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGt                                                                                                                  <  1:2602724/148‑1 (MQ=255)
                               cGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGt                                                                                                                  >  2:3402737/1‑148 (MQ=255)
                                gCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTg                                                                                                                 <  1:2690957/148‑1 (MQ=255)
                                  aCTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGaa                                                                                                               <  1:1523971/148‑1 (MQ=255)
                                   cTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGaaa                                                                                                              <  2:2073721/148‑1 (MQ=255)
                                      gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAggg                                                                                                           >  2:1287842/1‑148 (MQ=255)
                                          gcgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTa                                                                                                       <  2:2157016/148‑1 (MQ=255)
                                          gcgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTa                                                                                                       <  1:338692/148‑1 (MQ=255)
                                                gCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAAt                                                                                                 <  2:2814510/148‑1 (MQ=255)
                                                gCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAAt                                                                                                 <  1:934572/148‑1 (MQ=255)
                                                 cATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATa                                                                                                <  1:3097255/148‑1 (MQ=255)
                                                 cATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATa                                                                                                <  1:1719902/148‑1 (MQ=255)
                                                  aTGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATAc                                                                                               <  1:1273687/148‑1 (MQ=255)
                                                     gCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAACGGGTGTACTCCATACCaa                                                                                            >  2:3221634/1‑148 (MQ=255)
                                                      cACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAAt                                                                                           <  1:1485190/148‑1 (MQ=255)
                                                         cTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCAc                                                                                        <  2:2090140/148‑1 (MQ=255)
                                                         cTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCAc                                                                                        <  2:1325068/148‑1 (MQ=255)
                                                          tGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCAcc                                                                                       >  1:1850398/1‑148 (MQ=255)
                                                                       gCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGt                                                                          >  2:74858/1‑148 (MQ=255)
                                                                           aGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCt                                                                      >  1:374479/1‑148 (MQ=255)
                                                                           aGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCt                                                                      >  1:2299237/1‑148 (MQ=255)
                                                                            gCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCtt                                                                     >  2:2594369/1‑148 (MQ=255)
                                                                               aaTGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGcc                                                                  >  1:636708/1‑148 (MQ=255)
                                                                                           atTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatg                                                      >  1:503365/1‑148 (MQ=255)
                                                                                           atTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatg                                                      <  1:1175910/148‑1 (MQ=255)
                                                                                            ttgatttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatga                                                     <  1:1868707/148‑1 (MQ=255)
                                                                                            ttgatttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatga                                                     <  2:1030625/148‑1 (MQ=255)
                                                                                                tttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCg                                                 <  2:3100404/148‑1 (MQ=255)
                                                                                                 ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                <  2:3644761/148‑1 (MQ=255)
                                                                                                 ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                <  1:296661/148‑1 (MQ=255)
                                                                                                 ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                <  2:1752597/148‑1 (MQ=255)
                                                                                                 ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  1:3353282/1‑148 (MQ=255)
                                                                                                 ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  2:1677771/1‑148 (MQ=255)
                                                                                                    atAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCg                                             >  2:3194168/1‑148 (MQ=255)
                                                                                                     tAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGt                                            >  1:2786678/1‑148 (MQ=255)
                                                                                                     tAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGt                                            >  1:2786472/1‑148 (MQ=255)
                                                                                                      aaTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTg                                           <  1:1787995/148‑1 (MQ=255)
                                                                                                      aaTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTg                                           >  2:3175023/1‑148 (MQ=255)
                                                                                                       aTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTgg                                          <  2:3431372/148‑1 (MQ=255)
                                                                                                        tcctccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTggg                                         >  2:236534/1‑148 (MQ=255)
                                                                                                             cGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCg                                    <  1:2905412/148‑1 (MQ=255)
                                                                                                               ccTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgag                                  <  1:3039756/148‑1 (MQ=255)
                                                                                                               ccTGCTGGTGGAAATAACAACAAGCAGCCACAAGTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgag                                  <  2:1967235/148‑1 (MQ=255)
                                                                                                                ctgctgGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgaga                                 <  2:1340701/148‑1 (MQ=255)
                                                                                                                    tggtggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGc                             <  1:221239/148‑1 (MQ=255)
                                                                                                                        ggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAGGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCgag                         >  2:798473/1‑148 (MQ=255)
                                                                                                                          aaaTAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGaa                       >  1:1957274/1‑148 (MQ=255)
                                                                                                                               acaacaAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                <  2:2410410/118‑1 (MQ=255)
                                                                                                                               acaacaAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  1:2410410/1‑118 (MQ=255)
                                                                                                                               acaacaAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaaaa                  >  2:918484/1‑148 (MQ=255)
                                                                                                                               acaacaAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaaaa                  <  1:1765020/148‑1 (MQ=255)
                                                                                                                               acaacaAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaaaa                  <  1:641817/148‑1 (MQ=255)
                                                                                                                                  acaaGCAGCCACAATTACCTCAATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGGGCATGATGCGATCGTGGTTATCGAGATTGCCGAGAATAAAACCg               <  1:1125559/148‑1 (MQ=255)
                                                                                                                                   caaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGc              <  2:1053256/148‑1 (MQ=255)
                                                                                                                                   caaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGc              >  1:702575/1‑148 (MQ=255)
                                                                                                                                   caaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGc              >  1:10434/1‑148 (MQ=255)
                                                                                                                                       ccgcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAATGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTcc          <  1:3484276/146‑1 (MQ=255)
                                                                                                                                        agacacaAGTACCTCTATCGCTTGCGTAAAACGCTCTTCAGGTGAATGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCTGCAGGCGCATGATGCGATCGTGGTTATCGAGATTGCCTTGAATAAAACCGCCTCCa         <  2:3508251/145‑1 (MQ=255)
                                                                                                                                          ccACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATg       <  2:1859263/148‑1 (MQ=255)
                                                                                                                                          ccACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATg       <  2:1859242/148‑1 (MQ=255)
                                                                                                                                              aaTTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaa   >  2:3377382/1‑148 (MQ=255)
                                                                                                                                              aaTTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaa   <  2:1187502/148‑1 (MQ=255)
                                                                                                                                              aaTTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaa   <  1:611912/148‑1 (MQ=255)
                                                                                                                                               aTTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaat  >  1:2544609/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAAT  >  CP009273/494874‑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: ‑