Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,715,659 C→A 22.7% A783E (GCG→GAG)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,715,6590CA22.7% 626.0 / 159.5 250A783E (GCG→GAG) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (103/90);  new base A (27/30);  total (130/120)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 253 total aligned reads displayed.

AGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATT  >  CP009273/2715510‑2715807
                                                                                                                                                     |                                                                                                                                                    
agaagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGc                                                                                                                                                      <  2:3977331/150‑1 (MQ=255)
agaagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGc                                                                                                                                                      >  1:4095884/1‑150 (MQ=255)
  aagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCgg                                                                                                                                                    >  2:4460996/1‑150 (MQ=255)
     aTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTac                                                                                                                                                 <  1:368531/150‑1 (MQ=255)
             caACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTAtt                                                                                                                                         <  2:460168/150‑1 (MQ=255)
             caACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTAtt                                                                                                                                         >  2:4813952/1‑150 (MQ=255)
             caACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAtt                                                                                                                                         <  1:1490975/150‑1 (MQ=255)
                      cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACg                                                                                                                                >  2:294229/1‑150 (MQ=255)
                       aaTTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACgg                                                                                                                               >  2:4296516/1‑150 (MQ=255)
                       aaTTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACgg                                                                                                                               >  2:1872966/1‑150 (MQ=255)
                              ttCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGATTTGACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGAt                                                                                                                        <  1:3025382/150‑1 (MQ=255)
                                 gcgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCt                                                                                                                     >  2:1475857/1‑150 (MQ=255)
                                  cgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTc                                                                                                                    <  1:4317952/150‑1 (MQ=255)
                                  cgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTc                                                                                                                    <  1:1864103/150‑1 (MQ=255)
                                  cgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTc                                                                                                                    >  2:4237633/1‑150 (MQ=255)
                                   gAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTATGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCg                                                                                                                   >  1:1168936/1‑150 (MQ=255)
                                   gAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCg                                                                                                                   >  1:3809305/1‑150 (MQ=255)
                                    aGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTCTTGATCAAACGGAAGAGATCCTCgg                                                                                                                  >  1:2541205/1‑150 (MQ=255)
                                     gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGc                                                                                                                 <  1:4363662/150‑1 (MQ=255)
                                     gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGc                                                                                                                 <  1:3544821/150‑1 (MQ=255)
                                     gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGc                                                                                                                 <  1:3463605/150‑1 (MQ=255)
                                     gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGc                                                                                                                 >  1:1005486/1‑150 (MQ=255)
                                     gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGc                                                                                                                 <  2:74188/150‑1 (MQ=255)
                                      tCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCg                                                                                                                >  1:1221576/1‑150 (MQ=255)
                                         ccAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa                                                                                                             >  1:1251643/1‑150 (MQ=255)
                                           aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                                                           >  2:4750093/1‑150 (MQ=255)
                                             agaagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcgcg                                                                                                         >  1:1780809/1‑150 (MQ=255)
                                                agaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGc                                                                                                      >  1:1250458/1‑150 (MQ=255)
                                                agaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGc                                                                                                      >  2:2437024/1‑150 (MQ=255)
                                                   tCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAt                                                                                                   <  2:3322772/150‑1 (MQ=255)
                                                      ttattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTc                                                                                                >  1:3293165/1‑150 (MQ=255)
                                                       tattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTcc                                                                                               <  1:1358838/150‑1 (MQ=255)
                                                          taCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGAt                                                                                            >  2:2416203/1‑150 (MQ=255)
                                                              gCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTg                                                                                        <  1:2128008/150‑1 (MQ=255)
                                                                tACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGAt                                                                                      <  1:306539/150‑1 (MQ=255)
                                                                 aCTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATa                                                                                     <  1:3881553/150‑1 (MQ=255)
                                                                  cTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATaa                                                                                    >  2:2797172/1‑150 (MQ=255)
                                                                  cTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATaa                                                                                    >  2:3272969/1‑150 (MQ=255)
                                                                  cTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATaa                                                                                    >  2:3213655/1‑150 (MQ=255)
                                                                      aGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATc                                                                                <  2:1978755/150‑1 (MQ=255)
                                                                        cGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGa                                                                              >  2:1793896/1‑150 (MQ=255)
                                                                         gagaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGAt                                                                             >  2:2551848/1‑150 (MQ=255)
                                                                         gagaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGAt                                                                             <  1:3687533/150‑1 (MQ=255)
                                                                         gagaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGAt                                                                             >  2:3453579/1‑150 (MQ=255)
                                                                            aTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGcc                                                                          >  1:2806174/1‑150 (MQ=255)
                                                                            aTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGcc                                                                          >  2:3147554/1‑150 (MQ=255)
                                                                            aTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGcc                                                                          <  1:3674101/150‑1 (MQ=255)
                                                                                  gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAttt                                                                    <  1:4704158/150‑1 (MQ=255)
                                                                                  gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAttt                                                                    <  2:1369585/150‑1 (MQ=255)
                                                                                       tACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGt                                                               >  2:1858080/1‑150 (MQ=255)
                                                                                       tACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGt                                                               <  1:1667869/150‑1 (MQ=255)
                                                                                       tACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGt                                                               >  2:1634106/1‑150 (MQ=255)
                                                                                        aCCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTa                                                              <  2:3903245/150‑1 (MQ=255)
                                                                                          ccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACt                                                            >  1:490980/1‑150 (MQ=255)
                                                                                          ccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGTCGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACt                                                            <  1:1921095/150‑1 (MQ=255)
                                                                                             tGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGt                                                         >  1:1216692/1‑150 (MQ=255)
                                                                                              gAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGtt                                                        <  1:2321963/150‑1 (MQ=255)
                                                                                              gAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGtt                                                        <  1:3760312/150‑1 (MQ=255)
                                                                                                aGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCg                                                      <  1:3049141/150‑1 (MQ=255)
                                                                                                aGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCg                                                      >  2:2377699/1‑150 (MQ=255)
                                                                                                aGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCg                                                      >  2:3784898/1‑150 (MQ=255)
                                                                                                aGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCg                                                      >  2:2722119/1‑150 (MQ=255)
                                                                                                  aTTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCt                                                    >  1:3532349/1‑150 (MQ=255)
                                                                                                  aTTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCt                                                    <  2:3535351/150‑1 (MQ=255)
                                                                                                   tttAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTc                                                   <  2:3663698/150‑1 (MQ=255)
                                                                                                    ttAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCg                                                  >  2:1136328/1‑150 (MQ=255)
                                                                                                    ttAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCg                                                  <  2:1902356/150‑1 (MQ=255)
                                                                                                     tAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCgg                                                 <  2:1704069/150‑1 (MQ=255)
                                                                                                             aTGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa                                         >  1:25057/1‑150 (MQ=255)
                                                                                                             aTGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa                                         >  1:957624/1‑150 (MQ=255)
                                                                                                              tGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAg                                        >  1:1025465/1‑150 (MQ=255)
                                                                                                              tGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAg                                        >  1:1070713/1‑150 (MQ=255)
                                                                                                              tGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAg                                        >  2:1446611/1‑150 (MQ=255)
                                                                                                              tGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAgg                                       >  1:1360714/1‑150 (MQ=255)
                                                                                                                 cacaGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGt                                     >  1:1443469/1‑150 (MQ=255)
                                                                                                                   caGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTa                                   <  1:936415/150‑1 (MQ=255)
                                                                                                                   caGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTa                                   >  2:3492904/1‑150 (MQ=255)
                                                                                                                       tCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCt                               <  1:1751444/150‑1 (MQ=255)
                                                                                                                         gACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTa                             <  2:842654/150‑1 (MQ=255)
                                                                                                                            tACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGt                          <  1:2245471/150‑1 (MQ=255)
                                                                                                                              cGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCg                        >  2:4233317/1‑150 (MQ=255)
                                                                                                                              cGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCg                        <  1:3738648/150‑1 (MQ=255)
                                                                                                                               gATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGa                       >  2:4468437/1‑150 (MQ=255)
                                                                                                                                aTCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGAc                      <  1:3991487/150‑1 (MQ=255)
                                                                                                                                      aaaTGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTTGGCTTAGGTCGACGCTTaa                <  2:1589422/150‑1 (MQ=255)
                                                                                                                                      aaaTGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTaa                <  2:4047719/150‑1 (MQ=255)
                                                                                                                                      aaaTGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACACCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTaa                >  1:3909837/1‑150 (MQ=255)
                                                                                                                                       aaTGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAAt               <  2:3152733/150‑1 (MQ=255)
                                                                                                                                               ttGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGt       <  1:4460996/150‑1 (MQ=255)
                                                                                                                                               ttGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGACCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGt       >  2:2337418/1‑150 (MQ=255)
                                                                                                                                               ttGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGt       <  1:2437024/150‑1 (MQ=255)
                                                                                                                                                tGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGtt      <  1:4009203/150‑1 (MQ=255)
                                                                                                                                                tGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGtt      <  2:1170446/150‑1 (MQ=255)
                                                                                                                                                 gTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGt       >  2:4369186/1‑148 (MQ=255)
                                                                                                                                                  tAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGa    <  2:257953/150‑1 (MQ=255)
                                                                                                                                                  tAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGa    >  1:265616/1‑150 (MQ=255)
                                                                                                                                                  tAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGa    >  2:2110625/1‑150 (MQ=255)
                                                                                                                                                   aGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt   >  1:2981291/1‑150 (MQ=255)
                                                                                                                                                   aGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt   <  2:141948/150‑1 (MQ=255)
                                                                                                                                                   aGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCAGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGAt   <  2:2919570/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATT  >  CP009273/2715510‑2715807

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: