Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,715,643 G→T 100% E778* (GAA→TAA)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,715,6430GT97.2% 804.0 / 13.0 248E778* (GAA→TAA) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (132/109);  total (135/113)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 248 total aligned reads displayed.

CCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGC  >  CP009273/2715495‑2715789
                                                                                                                                                    |                                                                                                                                                  
ccGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGa                                                                                                                                                   >  2:732661/1‑150 (MQ=255)
  gCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTaaa                                                                                                                                                 <  1:126496/150‑1 (MQ=255)
  gCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTaaa                                                                                                                                                 <  2:534219/150‑1 (MQ=255)
     cGATTTTGCCAGAAGATTAGCCACAACTTCAGCAATTCATTTAGCGAGTCACCAAAGAAGATCTTTATTAACGCTACTTTAGAGAGATCAACGAATTTACCCATGAATATTTAGACAACATGACACAGATATACTACTATCGGTAAATgg                                                                                                                                              >  2:3005763/1‑150 (MQ=255)
       aTTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTCCGATCGGTAAATGGCg                                                                                                                                            >  2:783520/1‑150 (MQ=255)
        ttttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGt                                                                                                                                           <  1:443267/150‑1 (MQ=255)
          ttGCCAGAAGATGAGCCACCACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGttt                                                                                                                                         <  2:701245/150‑1 (MQ=255)
          ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGttt                                                                                                                                         <  2:2901721/150‑1 (MQ=255)
          ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGttt                                                                                                                                         <  2:1436212/150‑1 (MQ=255)
            gCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGt                                                                                                                                       <  1:96247/150‑1 (MQ=255)
             ccAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTa                                                                                                                                      <  2:3623656/150‑1 (MQ=255)
                 aagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCgg                                                                                                                                  >  1:4070387/1‑150 (MQ=255)
                  agaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGt                                                                                                                                 <  2:2095582/150‑1 (MQ=255)
                      gAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTacga                                                                                                                             >  2:1541516/1‑150 (MQ=255)
                       aGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTacgac                                                                                                                            <  2:3292520/150‑1 (MQ=255)
                           acaACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTAt                                                                                                                        <  2:1807999/150‑1 (MQ=255)
                            caACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTAtt                                                                                                                       >  1:447854/1‑150 (MQ=255)
                             aaCTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTg                                                                                                                      >  2:614582/1‑150 (MQ=255)
                              aCTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGa                                                                                                                     >  2:2154787/1‑150 (MQ=255)
                              aCTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGa                                                                                                                     <  1:714482/150‑1 (MQ=255)
                               cTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGAt                                                                                                                    >  1:116336/1‑150 (MQ=255)
                                  cAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCaa                                                                                                                 >  1:3176566/1‑150 (MQ=255)
                                  cAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCaa                                                                                                                 >  1:3480511/1‑150 (MQ=255)
                                   aGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATTAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCaaa                                                                                                                <  2:683029/150‑1 (MQ=255)
                                     cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACg                                                                                                              >  2:3434225/1‑150 (MQ=255)
                                     cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACg                                                                                                              >  2:3434214/1‑150 (MQ=255)
                                     cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACg                                                                                                              >  1:1853862/1‑150 (MQ=255)
                                       aTTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGa                                                                                                            <  1:1167278/150‑1 (MQ=255)
                                       aTTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGa                                                                                                            >  2:2959997/1‑150 (MQ=255)
                                          cATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAaga                                                                                                         >  1:3732569/1‑150 (MQ=255)
                                             ttCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGAt                                                                                                      <  2:1440394/150‑1 (MQ=255)
                                               cgcgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATcc                                                                                                    <  1:599952/150‑1 (MQ=255)
                                                 cgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTc                                                                                                  <  2:3095602/150‑1 (MQ=255)
                                                  gAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCg                                                                                                 >  1:3944078/1‑150 (MQ=255)
                                                  gAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCg                                                                                                 >  1:3350990/1‑150 (MQ=255)
                                                    gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCACACGGAAGAGATCCTCGGc                                                                                               >  1:4044822/1‑150 (MQ=255)
                                                     tCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCg                                                                                              >  1:3368549/1‑150 (MQ=255)
                                                       aCCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTc                                                                                            >  1:1920217/1‑150 (MQ=255)
                                                       aCCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTc                                                                                            >  1:2911007/1‑150 (MQ=255)
                                                          aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTCCGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                                         >  2:4453617/1‑150 (MQ=255)
                                                          aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                                         <  1:2255817/150‑1 (MQ=255)
                                                          aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                                         >  2:4124078/1‑150 (MQ=255)
                                                          aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                                         >  2:2664835/1‑150 (MQ=255)
                                                             gaagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTCCGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGt                                                                                      >  2:1117704/1‑150 (MQ=255)
                                                             gaagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGt                                                                                      >  1:4413077/1‑150 (MQ=255)
                                                              aagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                                     <  2:3682621/150‑1 (MQ=255)
                                                              aagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                                     <  2:3251425/150‑1 (MQ=255)
                                                                gaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCg                                                                                   >  2:826661/1‑150 (MQ=255)
                                                                 aTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGa                                                                                  <  1:547947/150‑1 (MQ=255)
                                                                   cTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAtt                                                                                <  2:3791276/150‑1 (MQ=255)
                                                                        ttaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGa                                                                           >  1:1646684/1‑150 (MQ=255)
                                                                        ttaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACACGTGCGATTTCCGa                                                                           >  1:3593078/1‑150 (MQ=255)
                                                                         taCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGAt                                                                          >  1:2400913/1‑150 (MQ=255)
                                                                         taCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGAt                                                                          >  2:414025/1‑150 (MQ=255)
                                                                          aCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATc                                                                         >  2:1740893/1‑150 (MQ=255)
                                                                           ccGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTGAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATcc                                                                        <  1:2725920/150‑1 (MQ=255)
                                                                              cTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGa                                                                     >  1:1799429/1‑150 (MQ=255)
                                                                               tACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGAt                                                                    >  1:3444792/1‑150 (MQ=255)
                                                                                aCTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATa                                                                   >  2:4568755/1‑150 (MQ=255)
                                                                                aCTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATa                                                                   >  1:3484097/1‑150 (MQ=255)
                                                                                  tttAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAAc                                                                 <  1:3036111/150‑1 (MQ=255)
                                                                                    tAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACAt                                                               >  2:1911770/1‑150 (MQ=255)
                                                                                     aGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATc                                                              >  2:2375679/1‑150 (MQ=255)
                                                                                        gagaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTCCGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGAt                                                           >  2:1219347/1‑150 (MQ=255)
                                                                                        gagaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGAt                                                           >  1:353600/1‑150 (MQ=255)
                                                                                        gagaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGAt                                                           >  1:3928556/1‑150 (MQ=255)
                                                                                             cAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGa                                                      <  2:4264797/150‑1 (MQ=255)
                                                                                                 gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAttt                                                  <  2:4000388/150‑1 (MQ=255)
                                                                                                  aaTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                                 >  1:1544210/1‑150 (MQ=255)
                                                                                                   aTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGc                                                >  2:3022580/1‑150 (MQ=255)
                                                                                                     ttACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTg                                              >  1:3213545/1‑150 (MQ=255)
                                                                                                         ccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACt                                          >  2:3948115/1‑150 (MQ=255)
                                                                                                         ccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACt                                          >  2:1623512/1‑150 (MQ=255)
                                                                                                         ccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACt                                          <  1:3754188/150‑1 (MQ=255)
                                                                                                          cATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTg                                         <  2:388776/150‑1 (MQ=255)
                                                                                                                      gCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAt                             <  1:1013464/150‑1 (MQ=255)
                                                                                                                        cAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtct                           <  2:580942/150‑1 (MQ=255)
                                                                                                                        cAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtct                           >  1:2684931/1‑150 (MQ=255)
                                                                                                                         aaCATGACACAGATCGACAACGATCGGTAAATGGCGTTTGTAGCGGTACGACGGATTGATCAAACGGAAGAGATCCTCGGCGACACGCGTGCGATTTCAGAACCTGATAACAACGATGCCGAATTAGCTGTACTCGTACGCTCGGAtctc                          >  2:2066268/1‑150 (MQ=255)
                                                                                                                          aCATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCa                         <  2:2311000/150‑1 (MQ=255)
                                                                                                                           cATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaa                        <  2:2838284/150‑1 (MQ=255)
                                                                                                                           cATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaa                        <  2:2838103/150‑1 (MQ=255)
                                                                                                                           cATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaa                        >  2:1664228/1‑150 (MQ=255)
                                                                                                                               acacAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAggg                    >  1:533623/1‑150 (MQ=255)
                                                                                                                                 acaGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGtt                  >  1:3280525/1‑150 (MQ=255)
                                                                                                                                   aGATCGACTCCGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                >  2:2834827/1‑150 (MQ=255)
                                                                                                                                   aGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                >  1:1702669/1‑150 (MQ=255)
                                                                                                                                     aTCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGc              >  2:271116/1‑150 (MQ=255)
                                                                                                                                     aTCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGc              <  2:2501769/150‑1 (MQ=255)
                                                                                                                                      tCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCt             >  1:1739713/1‑150 (MQ=255)
                                                                                                                                      tCGACTACGATCGGTAAATGGCGTTTGGAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCt             <  2:3145246/150‑1 (MQ=255)
                                                                                                                                        gACTACGATCGGTAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGACACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATATGCTGTACTGGTTCGCTCGGATCTCAAAGTGTTAGGCTTa           >  2:2607810/1‑150 (MQ=255)
                                                                                                                                              gATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGa     >  1:2223359/1‑150 (MQ=255)
                                                                                                                                              gATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGa     >  1:3248377/1‑150 (MQ=255)
                                                                                                                                              gATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGa     >  1:2413055/1‑150 (MQ=255)
                                                                                                                                               aTCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGAc    >  2:4533482/1‑150 (MQ=255)
                                                                                                                                               aTCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGAc    <  1:965685/150‑1 (MQ=255)
                                                                                                                                               aTCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGAc    <  1:745087/150‑1 (MQ=255)
                                                                                                                                                tCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACg   >  2:2915874/1‑150 (MQ=255)
                                                                                                                                                 cGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGc  >  2:2032733/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGC  >  CP009273/2715495‑2715789

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: