Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,046,958 T→A 6.9% E414D (GAA→GAT pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,046,9580TA 6.9% 1083.4 / 51.9 325E414D (GAA→GATpepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base T (147/154);  new base A (13/10);  total (161/164)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 332 total aligned reads displayed.

ACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACC  >  CP009273/3046809‑3047107
                                                                                                                                                     |                                                                                                                                                     
 cTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGttt                                                                                                                                                      <  2:1660174/150‑1 (MQ=255)
    aaCCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGAGTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGTTTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCgg                                                                                                                                                   >  1:1914327/1‑150 (MQ=255)
    aaCCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACAGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCgg                                                                                                                                                   <  1:285143/150‑1 (MQ=255)
      ccGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTa                                                                                                                                                 >  1:1068317/1‑150 (MQ=255)
        gggAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAAt                                                                                                                                               <  1:2888803/150‑1 (MQ=255)
         ggAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATc                                                                                                                                              >  2:3711968/1‑150 (MQ=255)
          gAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCCACGCTTCGAGTTCTTCCGGCTGTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGGTTCGGTAATca                                                                                                                                             >  2:518663/1‑150 (MQ=255)
            aaTAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTATCGGTAATcaca                                                                                                                                           >  1:694533/1‑150 (MQ=255)
             aTAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAAACACaa                                                                                                                                          >  2:105103/1‑150 (MQ=255)
                  cagcgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATgtcg                                                                                                                                     <  2:3606302/150‑1 (MQ=255)
                   agcgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATgtcgt                                                                                                                                    <  2:312604/150‑1 (MQ=255)
                   agcgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATgtcgt                                                                                                                                    >  2:49776/1‑150 (MQ=255)
                    gcgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATgtcgtc                                                                                                                                   <  2:1741597/150‑1 (MQ=255)
                     cgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCt                                                                                                                                  <  2:1536066/150‑1 (MQ=255)
                      gccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCtt                                                                                                                                 >  2:442250/1‑150 (MQ=255)
                      gccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCtt                                                                                                                                 <  1:1541989/150‑1 (MQ=255)
                       ccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTc                                                                                                                                >  2:3802865/1‑150 (MQ=255)
                        cagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCa                                                                                                                               <  2:4267597/150‑1 (MQ=255)
                            gTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTATTCAATAc                                                                                                                           <  2:4929899/150‑1 (MQ=255)
                                cgcccgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAAt                                                                                                                       <  2:205309/150‑1 (MQ=255)
                                  cccgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGGGGGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTGCAATACGAATGc                                                                                                                     >  1:2810062/1‑150 (MQ=255)
                                   ccgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGcc                                                                                                                    >  1:4063085/1‑150 (MQ=255)
                                   ccgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGcc                                                                                                                    <  1:1361072/150‑1 (MQ=255)
                                   ccgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTATCGGTAATCACAATGTCGTCTTCAATACGAATGcc                                                                                                                    <  2:2160345/150‑1 (MQ=255)
                                    cgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCg                                                                                                                   <  1:3886999/150‑1 (MQ=255)
                                    cgccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTATCGGTAATCACAATGTCGTCTTCAATACGAATGCCg                                                                                                                   <  2:4166869/150‑1 (MQ=255)
                                        aTGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATAc                                                                                                               <  2:1160781/150‑1 (MQ=255)
                                          gccgccACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACcg                                                                                                             <  1:3112070/150‑1 (MQ=255)
                                          gccgccACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACcg                                                                                                             <  1:2959026/150‑1 (MQ=255)
                                           ccgccACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACcgc                                                                                                            >  2:1876431/1‑150 (MQ=255)
                                               cACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGata                                                                                                        <  2:2719162/150‑1 (MQ=255)
                                               cACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTATCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGata                                                                                                        >  1:1646582/1‑150 (MQ=255)
                                                   gacgaTGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGt                                                                                                    >  1:1182371/1‑150 (MQ=255)
                                                    acgaTGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGtt                                                                                                   >  1:3674634/1‑150 (MQ=255)
                                                       atgatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTg                                                                                                <  2:3990346/150‑1 (MQ=255)
                                                       atgatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTg                                                                                                <  1:3389730/150‑1 (MQ=255)
                                                       atgatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTg                                                                                                <  2:2945503/150‑1 (MQ=255)
                                                       atgatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTg                                                                                                >  2:3900927/1‑150 (MQ=255)
                                                        tgatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTgg                                                                                               <  1:1787987/150‑1 (MQ=255)
                                                        tgatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTgg                                                                                               >  1:2310031/1‑150 (MQ=255)
                                                         gatTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGc                                                                                              <  1:579106/150‑1 (MQ=255)
                                                             aCGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACtt                                                                                          <  2:1261064/150‑1 (MQ=255)
                                                             aCGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATAGTGTTCTGGCACtt                                                                                          >  1:2395434/1‑150 (MQ=255)
                                                               gCTCATTGCTGTCTCGCAGCAACCAACCACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATAACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCt                                                                                        <  1:3985067/150‑1 (MQ=255)
                                                                cTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTg                                                                                       <  1:922287/150‑1 (MQ=255)
                                                                cTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTg                                                                                       >  2:2738440/1‑150 (MQ=255)
                                                                  cATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCa                                                                                     >  1:1563511/1‑150 (MQ=255)
                                                                  cATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCa                                                                                     <  2:5127892/150‑1 (MQ=255)
                                                                   aTTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCAt                                                                                    <  1:5051089/150‑1 (MQ=255)
                                                                    ttGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATc                                                                                   <  1:438/150‑1 (MQ=255)
                                                                    ttGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATc                                                                                   >  1:3102942/1‑150 (MQ=255)
                                                                     tGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATcc                                                                                  <  2:461887/150‑1 (MQ=255)
                                                                         ttCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGgcg                                                                              <  1:2250995/150‑1 (MQ=255)
                                                                          tctcGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGgcgc                                                                             >  2:2797386/1‑150 (MQ=255)
                                                                           ctcGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCa                                                                            <  1:2479443/150‑1 (MQ=255)
                                                                              gcagcaACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAata                                                                         <  1:621366/150‑1 (MQ=255)
                                                                                agcaACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAAtata                                                                       >  2:4338131/1‑150 (MQ=255)
                                                                                  caACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACa                                                                     >  2:3626336/1‑150 (MQ=255)
                                                                                  caACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACa                                                                     <  2:4473375/150‑1 (MQ=255)
                                                                                   aaCCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAg                                                                    >  1:4995943/1‑150 (MQ=255)
                                                                                    aCCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGc                                                                   >  1:2121303/1‑150 (MQ=255)
                                                                                    aCCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGc                                                                   <  2:4883167/150‑1 (MQ=255)
                                                                                     ccATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGcc                                                                  >  1:817872/1‑150 (MQ=255)
                                                                                      cATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGTGCAATATACAGccc                                                                 <  2:2534754/150‑1 (MQ=255)
                                                                                       aTCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCt                                                                <  2:1097995/150‑1 (MQ=255)
                                                                                       aTCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCt                                                                <  1:3746473/150‑1 (MQ=255)
                                                                                        tCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTg                                                               <  2:1575285/150‑1 (MQ=255)
                                                                                         cAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTgg                                                              <  2:1777801/150‑1 (MQ=255)
                                                                                          aaCGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGc                                                             <  1:3549552/150‑1 (MQ=255)
                                                                                           aCGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGct                                                            <  2:3063880/150‑1 (MQ=255)
                                                                                             gCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGctct                                                          >  1:2383735/1‑150 (MQ=255)
                                                                                              cTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTa                                                         >  1:4205758/1‑150 (MQ=255)
                                                                                              cTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTa                                                         >  1:764058/1‑150 (MQ=255)
                                                                                                  gATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGt                                                     >  2:211781/1‑150 (MQ=255)
                                                                                                   aTTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTc                                                    >  1:3306227/1‑150 (MQ=255)
                                                                                                    tttCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCa                                                   <  1:2676498/150‑1 (MQ=255)
                                                                                                      tcttcCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTATCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGt                                                 <  1:3231586/150‑1 (MQ=255)
                                                                                                       cttcCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTa                                                >  2:200479/1‑150 (MQ=255)
                                                                                                        ttcCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTAc                                               >  1:956372/1‑150 (MQ=255)
                                                                                                         tcCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTAcc                                              <  1:3220421/150‑1 (MQ=255)
                                                                                                          ccGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCa                                             <  2:856095/150‑1 (MQ=255)
                                                                                                           cGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCAt                                            >  2:332275/1‑150 (MQ=255)
                                                                                                           cGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCAt                                            >  1:4579425/1‑150 (MQ=255)
                                                                                                            ggCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATg                                           >  1:4197690/1‑150 (MQ=255)
                                                                                                              cTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGcc                                         >  1:4734739/1‑150 (MQ=255)
                                                                                                               tttttCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGccc                                        <  1:4873392/150‑1 (MQ=255)
                                                                                                                     accacGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTcc                                  >  2:55927/1‑150 (MQ=255)
                                                                                                                         cGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATACCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGaa                              >  1:2532783/1‑150 (MQ=255)
                                                                                                                            tGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATgc                           >  2:3271741/1‑150 (MQ=255)
                                                                                                                                 gTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCgcgag                      <  2:2993156/150‑1 (MQ=255)
                                                                                                                                  tGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCgcgagc                     <  1:4971499/150‑1 (MQ=255)
                                                                                                                                   gAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCgcgagcg                    >  1:1698782/1‑150 (MQ=255)
                                                                                                                                      gTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATc                 <  1:3900927/150‑1 (MQ=255)
                                                                                                                                       ttttCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATcc                <  2:710097/150‑1 (MQ=255)
                                                                                                                                        tttCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCt               <  2:822841/150‑1 (MQ=255)
                                                                                                                                         ttCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTg              <  1:2213457/150‑1 (MQ=255)
                                                                                                                                            gTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGAcc           <  1:836700/150‑1 (MQ=255)
                                                                                                                                                 cGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATaaa      >  2:1981158/1‑150 (MQ=255)
                                                                                                                                                 cGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATaaa      >  1:1017333/1‑150 (MQ=255)
                                                                                                                                                     ttCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACAcc  >  2:3399401/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACC  >  CP009273/3046809‑3047107

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: