Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,715,643 G→T 87.1% 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,6430GT87.1% 989.2 / 142.1 371E778* (GAA→TAA) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/18);  new base T (183/140);  total (213/158)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 373 total aligned reads displayed.

CGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCT  >  CP009273/2715496‑2715790
                                                                                                                                                   |                                                                                                                                                   
cGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGaa                                                                                                                                                   <  1:5066566/150‑1 (MQ=255)
 gCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTaaa                                                                                                                                                  <  2:4804857/150‑1 (MQ=255)
 gCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTaaa                                                                                                                                                  >  1:1845882/1‑150 (MQ=255)
  cccGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAAt                                                                                                                                                 <  1:5543352/150‑1 (MQ=255)
   ccGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATg                                                                                                                                                <  2:1654226/150‑1 (MQ=255)
     gATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGc                                                                                                                                              >  2:1859820/1‑150 (MQ=255)
      aTTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCg                                                                                                                                             <  1:689268/150‑1 (MQ=255)
      aTTTTGCCAGAAGATGAGCCACAACATCAGCAATACATTTAGCGAGTCACCAAAGAAGAACTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGAATTATCCAACATGACACAGATCGACTCAGATCGGTAAATGGCg                                                                                                                                             >  1:3314306/1‑150 (MQ=255)
          tGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTg                                                                                                                                         <  1:4479181/150‑1 (MQ=255)
           gCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGt                                                                                                                                        <  2:2734690/150‑1 (MQ=255)
             cAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAg                                                                                                                                      >  1:103928/1‑150 (MQ=255)
             cAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAg                                                                                                                                      <  1:3160598/150‑1 (MQ=255)
                      aGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTacgac                                                                                                                             <  2:1205599/150‑1 (MQ=255)
                       gCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTacgacg                                                                                                                            <  2:782430/150‑1 (MQ=255)
                        ccACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGt                                                                                                                           >  2:1068796/1‑150 (MQ=255)
                              cTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                                                                                                                     <  2:1199842/150‑1 (MQ=255)
                               ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATc                                                                                                                    >  2:1040865/1‑150 (MQ=255)
                               ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATc                                                                                                                    <  1:5161287/150‑1 (MQ=255)
                                 cAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCaa                                                                                                                  >  1:3523000/1‑150 (MQ=255)
                                  aGCAATTCATTTCGCGAGTCACCAAATAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCaaa                                                                                                                 <  2:5618063/150‑1 (MQ=255)
                                       ttCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGaa                                                                                                            <  1:529484/150‑1 (MQ=255)
                                        tCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAag                                                                                                           >  2:1731719/1‑150 (MQ=255)
                                        tCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAag                                                                                                           <  1:2104688/150‑1 (MQ=255)
                                          aTTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAagag                                                                                                         >  1:4288614/1‑150 (MQ=255)
                                          aTTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAagag                                                                                                         <  2:543741/150‑1 (MQ=255)
                                              cgcgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATcc                                                                                                     >  1:2658470/1‑150 (MQ=255)
                                                  aGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCgg                                                                                                 >  2:2255081/1‑150 (MQ=255)
                                                   gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGc                                                                                                >  2:4483525/1‑150 (MQ=255)
                                                       ccAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa                                                                                            >  1:2932458/1‑150 (MQ=255)
                                                       ccAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa                                                                                            >  2:1840582/1‑150 (MQ=255)
                                                        cAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAc                                                                                           <  1:4979111/150‑1 (MQ=255)
                                                         aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                                          >  1:1732636/1‑150 (MQ=255)
                                                             aagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                                      <  2:94812/150‑1 (MQ=255)
                                                             aagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                                      >  2:3930563/1‑150 (MQ=255)
                                                             aagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAATCGGAAGAGATCCTCGGCGTCACGCGTg                                                                                      >  2:2276880/1‑150 (MQ=255)
                                                             aagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                                      >  2:3120832/1‑150 (MQ=255)
                                                                aTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGa                                                                                   <  1:3636249/150‑1 (MQ=255)
                                                                 tCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAt                                                                                  <  1:224525/150‑1 (MQ=255)
                                                                 tCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAt                                                                                  <  1:1982466/150‑1 (MQ=255)
                                                                  cTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATTGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAtt                                                                                 >  2:3460395/1‑150 (MQ=255)
                                                                  cTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAtt                                                                                 <  1:1728358/150‑1 (MQ=255)
                                                                  cTTGATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAtt                                                                                 <  2:4115157/150‑1 (MQ=255)
                                                                    ttattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCTTCACGCGTGCGATTTc                                                                               >  2:3626854/1‑150 (MQ=255)
                                                                    ttattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTc                                                                               >  1:909481/1‑150 (MQ=255)
                                                                    ttattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTc                                                                               >  1:3733230/1‑150 (MQ=255)
                                                                     tattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTcc                                                                              >  1:4177666/1‑150 (MQ=255)
                                                                        taCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGAt                                                                           >  2:5525000/1‑150 (MQ=255)
                                                                         aCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATc                                                                          <  1:3394063/150‑1 (MQ=255)
                                                                           cGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCt                                                                        >  2:2196910/1‑150 (MQ=255)
                                                                             cTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGa                                                                      >  1:1108863/1‑150 (MQ=255)
                                                                               aCTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATa                                                                    <  1:2867116/150‑1 (MQ=255)
                                                                                  ttAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACa                                                                 <  1:4760258/150‑1 (MQ=255)
                                                                                    aGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATc                                                               >  2:5490646/1‑150 (MQ=255)
                                                                                    aGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACAACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTAATGATCAAACGGAAGAGATCCTCGGCGTCAAGCGTGCGATTTCCGTTCCTGAAAACATc                                                               >  1:5202267/1‑150 (MQ=255)
                                                                                       gagaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGAt                                                            >  2:3347733/1‑150 (MQ=255)
                                                                                        agaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATg                                                           >  1:366975/1‑150 (MQ=255)
                                                                                              aCGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAt                                                     >  1:4842585/1‑150 (MQ=255)
                                                                                               cGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAtt                                                    >  2:5540831/1‑150 (MQ=255)
                                                                                                gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAttt                                                   <  2:1579056/150‑1 (MQ=255)
                                                                                                gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAttt                                                   <  1:606757/150‑1 (MQ=255)
                                                                                                gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGAAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCTGCGTCACGCGTGCGATTTCCGTTCCTGATAACATCGATGCCGAAttt                                                   >  2:1395950/1‑150 (MQ=255)
                                                                                                 aaTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGACTAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                                  >  2:5223502/1‑150 (MQ=255)
                                                                                                     tACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGt                                              <  2:4365019/150‑1 (MQ=255)
                                                                                                      aCCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTa                                             >  1:1499375/1‑150 (MQ=255)
                                                                                                       cccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTAc                                            >  2:505723/1‑150 (MQ=255)
                                                                                                         cATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTg                                          <  1:3941043/150‑1 (MQ=255)
                                                                                                         cATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTg                                          <  1:1376358/150‑1 (MQ=255)
                                                                                                          aTGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTgg                                         >  2:845846/1‑150 (MQ=255)
                                                                                                            gAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGtt                                       >  2:536868/1‑150 (MQ=255)
                                                                                                            gAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGtt                                       <  1:2272864/150‑1 (MQ=255)
                                                                                                                aTTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCt                                   <  2:4597678/150‑1 (MQ=255)
                                                                                                                aTTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCt                                   <  1:1208260/150‑1 (MQ=255)
                                                                                                                   tAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCgg                                <  2:2831608/150‑1 (MQ=255)
                                                                                                                   tAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTAGTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATACCATCGATGCCGAATTTGCTGTACTGGTTCGCTCgg                                >  2:28726/1‑150 (MQ=255)
                                                                                                                    aGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGa                               >  1:3845366/1‑150 (MQ=255)
                                                                                                                     gCCAACATGCCACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAt                              <  2:469982/150‑1 (MQ=255)
                                                                                                                       cAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGGTCGCTCGGAtct                            >  2:3782789/1‑150 (MQ=255)
                                                                                                                        aaCATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtctc                           >  2:245588/1‑150 (MQ=255)
                                                                                                                        aaCATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtctc                           >  2:2412658/1‑150 (MQ=255)
                                                                                                                         aCATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCa                          >  2:2627912/1‑150 (MQ=255)
                                                                                                                          cATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaa                         >  1:2353175/1‑150 (MQ=255)
                                                                                                                          cATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaa                         <  1:774205/150‑1 (MQ=255)
                                                                                                                           aTGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa                        >  2:2915871/1‑150 (MQ=255)
                                                                                                                           aTGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGAACAAAAGGAAGAGATCATCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATAGATGCCGAAAATTCTGTACTGGTTCGCTCGGATCTCaaa                        >  1:3780445/1‑150 (MQ=255)
                                                                                                                              acacAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAggg                     >  2:486198/1‑150 (MQ=255)
                                                                                                                               cacaGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGt                    >  1:1573612/1‑150 (MQ=255)
                                                                                                                                  aGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                 >  2:3757322/1‑150 (MQ=255)
                                                                                                                                  aGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                 <  1:5490646/150‑1 (MQ=255)
                                                                                                                                  aGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                 >  2:1873514/1‑150 (MQ=255)
                                                                                                                                   gATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAgg                >  1:1487102/1‑150 (MQ=255)
                                                                                                                                     tCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCt              <  1:5146453/150‑1 (MQ=255)
                                                                                                                                        aCTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAg           >  1:5342489/1‑150 (MQ=255)
                                                                                                                                         cTACGATCGGTAAATGGCGTTTGTAGCGGTACGAAGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAgg          <  1:2999309/150‑1 (MQ=255)
                                                                                                                                          tACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGt         <  2:3221495/150‑1 (MQ=255)
                                                                                                                                          tACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGt         <  2:889668/150‑1 (MQ=255)
                                                                                                                                          tACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGt         <  2:2065419/150‑1 (MQ=255)
                                                                                                                                           aCGATCGGTTAATGGCGTTTGTAGCGGTACGACGTAGTGATCAAAAGGAAGAGATCCTCGGAGTCACGCGTGCGATTTCCGATCCTGATAACATCGATTCCGAATTTGCTGTACTGGGTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTc        <  2:3076662/150‑1 (MQ=255)
                                                                                                                                           aCGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTc        <  2:3996153/150‑1 (MQ=255)
                                                                                                                                           aCGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTc        <  2:3712171/150‑1 (MQ=255)
                                                                                                                                                 gggAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCt  >  2:29948/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCT  >  CP009273/2715496‑2715790

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: