Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,024,827 G→T 69.6% D103Y (GAT→TAT)  pepQ → proline dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,024,8270GT69.6% 302.9 / 305.9 290D103Y (GAT→TAT) pepQproline dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (53/35);  new base T (106/96);  total (159/131)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.60e-01

Only 100 of 296 total aligned reads displayed.

CCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGTTATC  >  CP009273/4024680‑4024976
                                                                                                                                                   |                                                                                                                                                     
ccGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAatat                                                                                                                                                     <  2:104472/150‑1 (MQ=255)
ccGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAatat                                                                                                                                                     <  1:379059/150‑1 (MQ=255)
  gCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtgt                                                                                                                                                   >  1:4207133/1‑150 (MQ=255)
  gCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtgt                                                                                                                                                   >  1:1958227/1‑150 (MQ=255)
     aTTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                                <  1:2563029/150‑1 (MQ=255)
       tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAg                                                                                                                                              <  1:1727256/150‑1 (MQ=255)
       tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAg                                                                                                                                              <  1:1994264/150‑1 (MQ=255)
         aaaGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTg                                                                                                                                            <  1:3492330/150‑1 (MQ=255)
           aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                          <  2:961706/150‑1 (MQ=255)
           aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                          >  2:1755012/1‑150 (MQ=255)
            gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATc                                                                                                                                         <  2:1537006/150‑1 (MQ=255)
            gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATc                                                                                                                                         <  2:2681706/150‑1 (MQ=255)
            gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATc                                                                                                                                         >  1:691576/1‑150 (MQ=255)
            gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATc                                                                                                                                         <  1:4171363/150‑1 (MQ=255)
            gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATc                                                                                                                                         <  1:1357476/150‑1 (MQ=255)
             cGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATcg                                                                                                                                        <  2:2446456/150‑1 (MQ=255)
              gTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATcgc                                                                                                                                       <  2:783885/150‑1 (MQ=255)
                 ggTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCt                                                                                                                                    <  1:2385833/150‑1 (MQ=255)
                  gTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTg                                                                                                                                   >  1:1001627/1‑150 (MQ=255)
                  gTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTg                                                                                                                                   >  1:41347/1‑150 (MQ=255)
                  gTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTg                                                                                                                                   >  1:995078/1‑150 (MQ=255)
                  gTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTg                                                                                                                                   >  2:3196352/1‑150 (MQ=255)
                    gCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGACGTGATCGCGCTGcc                                                                                                                                 <  1:2220454/150‑1 (MQ=255)
                    gCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGcc                                                                                                                                 <  1:4030198/150‑1 (MQ=255)
                       ggTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGaa                                                                                                                              >  2:2178148/1‑150 (MQ=255)
                        gTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGaaa                                                                                                                             <  1:955724/150‑1 (MQ=255)
                             tCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGa                                                                                                                        >  2:43340/1‑150 (MQ=255)
                             tCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGa                                                                                                                        >  2:43419/1‑150 (MQ=255)
                                    ccAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCAtt                                                                                                                 <  1:4208012/150‑1 (MQ=255)
                                    ccAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCAtt                                                                                                                 <  1:1934942/150‑1 (MQ=255)
                                     cAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTg                                                                                                                <  1:2121238/150‑1 (MQ=255)
                                      aaaCTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAAAGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTgg                                                                                                               >  1:2937678/1‑150 (MQ=255)
                                       aaCTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGt                                                                                                              >  2:2272330/1‑150 (MQ=255)
                                          tgctggttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgc                                                                                                           <  1:1592342/150‑1 (MQ=255)
                                          tgctggttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgc                                                                                                           >  1:1758246/1‑150 (MQ=255)
                                             tggttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                                                                                        <  2:2615144/150‑1 (MQ=255)
                                              ggttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgc                                                                                                       >  2:4151641/1‑150 (MQ=255)
                                                 tgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgc                                                                                                    >  2:4542392/1‑150 (MQ=255)
                                                    tggtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctg                                                                                                 >  1:306827/1‑150 (MQ=255)
                                                     ggtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctgc                                                                                                <  2:4216959/150‑1 (MQ=255)
                                                     ggtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctgc                                                                                                >  2:4309533/1‑150 (MQ=255)
                                                      gtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctgct                                                                                               >  1:4032433/1‑150 (MQ=255)
                                                        ggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCctgctgcctgctgc                                                                                             >  2:2976122/1‑150 (MQ=255)
                                                        ggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCctgctgcctgctgc                                                                                             >  2:34333/1‑150 (MQ=255)
                                                         gATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcg                                                                                            >  2:3346187/1‑150 (MQ=255)
                                                         gATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcg                                                                                            >  2:361439/1‑150 (MQ=255)
                                                           tGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgcg                                                                                          >  1:908570/1‑150 (MQ=255)
                                                           tGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgcg                                                                                          <  2:1476425/150‑1 (MQ=255)
                                                            ggCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGAATACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGAGATCGCGCAGCCGAAAGCCGATGGCATTGGAAGCCAGCAGACTGCTgcgcgc                                                                                         >  2:2645835/1‑150 (MQ=255)
                                                                tGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCa                                                                                     <  2:1859149/150‑1 (MQ=255)
                                                                tGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCa                                                                                     <  2:1824799/150‑1 (MQ=255)
                                                                 gAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCaa                                                                                    >  1:2126156/1‑150 (MQ=255)
                                                                  aaCAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAat                                                                                   >  1:1426427/1‑150 (MQ=255)
                                                                     aaGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATc                                                                                >  2:282857/1‑150 (MQ=255)
                                                                     aaGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATc                                                                                >  1:877870/1‑150 (MQ=255)
                                                                          gAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTta                                                                           >  1:4283112/1‑150 (MQ=255)
                                                                          gAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTta                                                                           >  2:3131005/1‑150 (MQ=255)
                                                                              cTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATc                                                                       >  2:4429554/1‑150 (MQ=255)
                                                                                 tgGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                                    >  2:1636912/1‑150 (MQ=255)
                                                                                 tgGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                                    >  2:3896270/1‑150 (MQ=255)
                                                                                 tgGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                                    <  1:2273413/150‑1 (MQ=255)
                                                                                 tgGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                                    >  2:444431/1‑150 (MQ=255)
                                                                                  ggTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTc                                                                   >  1:2180104/1‑150 (MQ=255)
                                                                                  ggTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTc                                                                   <  2:4131431/150‑1 (MQ=255)
                                                                                    ttCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCg                                                                 <  1:1135090/150‑1 (MQ=255)
                                                                                        aTCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGc                                                             >  2:4053283/1‑150 (MQ=255)
                                                                                          cTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCg                                                           <  1:3980518/150‑1 (MQ=255)
                                                                                           tGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCgg                                                          >  2:1897530/1‑150 (MQ=255)
                                                                                             ccGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGAATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCTGTTATATCGGTCCGGTGCCGGac                                                        >  2:3496978/1‑149 (MQ=255)
                                                                                              cGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAAc                                                       <  2:2016402/150‑1 (MQ=255)
                                                                                                  tGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTgc                                                   <  2:3796409/150‑1 (MQ=255)
                                                                                                  tGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTgc                                                   <  1:3196352/150‑1 (MQ=255)
                                                                                                   gATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTgcg                                                  <  2:2713742/150‑1 (MQ=255)
                                                                                                          ggCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAAc                                           >  2:75481/1‑150 (MQ=255)
                                                                                                               aaCGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGt                                      >  1:3989615/1‑150 (MQ=255)
                                                                                                                 cGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATAGCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAt                                    >  2:3005627/1‑150 (MQ=255)
                                                                                                                  gTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAtt                                   <  2:1667240/150‑1 (MQ=255)
                                                                                                                  gTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAtt                                   >  2:3573539/1‑150 (MQ=255)
                                                                                                                   tCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTg                                  >  1:2808836/1‑150 (MQ=255)
                                                                                                                   tCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTg                                  <  1:1707661/150‑1 (MQ=255)
                                                                                                                     gAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAg                                <  1:2437792/150‑1 (MQ=255)
                                                                                                                     gAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAg                                >  2:1717276/1‑150 (MQ=255)
                                                                                                                      aaCCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAgg                               <  2:895673/150‑1 (MQ=255)
                                                                                                                        ccGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGcc                             >  1:1787593/1‑150 (MQ=255)
                                                                                                                         cGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCa                            >  1:3637811/1‑150 (MQ=255)
                                                                                                                           cTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCACCTGGGTATTGAGGCCAGc                          >  2:2082057/1‑150 (MQ=255)
                                                                                                                           cTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGc                          >  2:1737513/1‑150 (MQ=255)
                                                                                                                              ccGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAAc                       <  1:581685/150‑1 (MQ=255)
                                                                                                                               cGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACa                      >  2:1286692/1‑150 (MQ=255)
                                                                                                                                gACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACAt                     <  2:3768006/150‑1 (MQ=255)
                                                                                                                                    tctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAAc                 <  1:1877863/150‑1 (MQ=255)
                                                                                                                                     ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAAcc                >  2:2838715/1‑150 (MQ=255)
                                                                                                                                       ttCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTAATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCg              >  2:2914657/1‑150 (MQ=255)
                                                                                                                                        tCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGa             <  2:1179289/150‑1 (MQ=255)
                                                                                                                                          tGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGaaa           <  1:4542392/150‑1 (MQ=255)
                                                                                                                                          tGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGaaa           >  2:2563612/1‑150 (MQ=255)
                                                                                                                                              cTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgggg       <  2:3252413/150‑1 (MQ=255)
                                                                                                                                               tGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGt      <  1:3614864/150‑1 (MQ=255)
                                                                                                                                                gAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGtt     <  1:2043926/150‑1 (MQ=255)
                                                                                                                                                   gATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGTTATc  <  2:601030/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGGTTATC  >  CP009273/4024680‑4024976

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: