Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,550,177 T→A 8.5% Q11L (CAG→CTG)  rtcA ← RNA 3'‑terminal phosphate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,550,1770TA 8.5% 1198.1 / 73.2 362Q11L (CAG→CTG) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (155/176);  new base A (18/13);  total (173/189)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 374 total aligned reads displayed.

TTTCCGTCGCCGCCTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTTCGTCGATCACT  >  CP009273/3550028‑3550326
                                                                                                                                                     |                                                                                                                                                     
 ttCCGTCGCCGCCTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCtg                                                                                                                                                      <  2:3528387/150‑1 (MQ=255)
  tCCGTCGCCGCCTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCtgt                                                                                                                                                     <  1:2034259/150‑1 (MQ=255)
    cGTCGCCGCCTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTgc                                                                                                                                                   <  2:327121/150‑1 (MQ=255)
    cGTCGCCGCCTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCAGTgc                                                                                                                                                   <  1:430604/150‑1 (MQ=255)
        gccgccTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCa                                                                                                                                               >  1:3484954/1‑150 (MQ=255)
          cgccTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATc                                                                                                                                             >  2:760259/1‑150 (MQ=255)
                  ccGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCCGCGCaa                                                                                                                                     >  2:2444209/1‑150 (MQ=255)
                  ccGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCaa                                                                                                                                     >  1:145242/1‑150 (MQ=255)
                    gcgGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAatc                                                                                                                                   <  1:3673327/150‑1 (MQ=255)
                      ggTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAatcat                                                                                                                                 >  1:3641550/1‑150 (MQ=255)
                       gTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAatcatc                                                                                                                                >  2:536829/1‑150 (MQ=255)
                         cAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCt                                                                                                                              >  2:4303500/1‑150 (MQ=255)
                             tGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTtc                                                                                                                          >  1:90537/1‑150 (MQ=255)
                              gCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGGGCGCCATCCAGCGCAATCATCCTTTtca                                                                                                                         <  2:2581838/150‑1 (MQ=255)
                               cTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTtcat                                                                                                                        >  2:1106951/1‑150 (MQ=255)
                               cTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATCGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTtcat                                                                                                                        <  2:2120968/150‑1 (MQ=255)
                                tGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTtcatc                                                                                                                       >  1:4443104/1‑150 (MQ=255)
                                   cgcAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCAtt                                                                                                                    <  2:4793798/150‑1 (MQ=255)
                                        cAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCCGCGCACTCATCCTGTTCATCATTTATcc                                                                                                               >  2:4744118/1‑150 (MQ=255)
                                             ccGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTa                                                                                                          <  2:1669958/150‑1 (MQ=255)
                                              cGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTAc                                                                                                         <  1:3272112/150‑1 (MQ=255)
                                                 tttCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGca                                                                                                      >  1:2820193/1‑150 (MQ=255)
                                                   tCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCCATCATCCTTTTCATCATTTATCCTTTTACGcaca                                                                                                    >  2:4373131/1‑150 (MQ=255)
                                                    cgcccgccCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCacac                                                                                                   <  2:4858119/150‑1 (MQ=255)
                                                    cgcccgccCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCacac                                                                                                   <  1:4581884/150‑1 (MQ=255)
                                                      cccgcccgGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCAcacca                                                                                                 <  2:2208913/150‑1 (MQ=255)
                                                      cccgcccgGCACGAAGGCTGGTGATGGTAAATGGCGGGCCGGGTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCAGTGCGCCATCCCGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCAcacca                                                                                                 <  2:3891515/150‑1 (MQ=255)
                                                       ccgcccgGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTACGCCCTGTGCGCCAGCCAGCGCAATCCTCCTTTTCAGCATTTATCCTTTTACGCAcaccac                                                                                                >  2:180517/1‑150 (MQ=255)
                                                        cgcccgGCACGTATGCTGGTGTTTCTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCTAACTATCTGCCCGCCACCTTCGCCCTGTGCGCCATCAAGCGCAATCAACCTTTTCATCATTTATCCTTTTACGCACaccacc                                                                                               >  2:580133/1‑150 (MQ=255)
                                                         gcccgGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCt                                                                                              >  2:3104079/1‑150 (MQ=255)
                                                         gcccgGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCt                                                                                              >  1:30522/1‑150 (MQ=255)
                                                          cccgGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTg                                                                                             <  1:1701833/150‑1 (MQ=255)
                                                          cccgGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTg                                                                                             >  2:1188074/1‑150 (MQ=255)
                                                            cgGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCCCTCATCCTTTTCATCATTTATCCTTTTACGCACACCACCGGAc                                                                                           >  2:5029275/1‑150 (MQ=255)
                                                               cACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCa                                                                                        >  2:1474988/1‑150 (MQ=255)
                                                               cACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCa                                                                                        <  2:83888/150‑1 (MQ=255)
                                                               cACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCa                                                                                        <  1:1465123/150‑1 (MQ=255)
                                                               cACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCa                                                                                        <  1:1051021/150‑1 (MQ=255)
                                                               cACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCa                                                                                        >  2:1762837/1‑150 (MQ=255)
                                                                 cGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAgg                                                                                      >  1:5027242/1‑150 (MQ=255)
                                                                  gAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAggg                                                                                     >  2:1362086/1‑150 (MQ=255)
                                                                  gAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAggg                                                                                     >  2:3447710/1‑150 (MQ=255)
                                                                   aaTGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGt                                                                                    >  2:20138/1‑150 (MQ=255)
                                                                    aTGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGtc                                                                                   >  1:4553215/1‑149 (MQ=255)
                                                                      gCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGtata                                                                                 >  1:4015220/1‑150 (MQ=255)
                                                                         ggTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGAt                                                                              >  1:1209749/1‑150 (MQ=255)
                                                                          gTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATa                                                                             >  1:4066318/1‑150 (MQ=255)
                                                                          gTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATa                                                                             >  1:2222170/1‑150 (MQ=255)
                                                                                gTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTcc                                                                       >  1:4981039/1‑150 (MQ=255)
                                                                                 tAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTcca                                                                      <  2:922417/150‑1 (MQ=255)
                                                                                 tAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTcca                                                                      <  1:1150787/150‑1 (MQ=255)
                                                                                 tAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTcca                                                                      >  2:4233431/1‑150 (MQ=255)
                                                                                   aaTGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTccacc                                                                    >  1:4847002/1‑150 (MQ=255)
                                                                                   aaTGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTccacc                                                                    >  2:2825608/1‑150 (MQ=255)
                                                                                   aaTGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTccacc                                                                    >  1:4182741/1‑150 (MQ=255)
                                                                                    aTGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTccacca                                                                   <  2:4484296/150‑1 (MQ=255)
                                                                                    aTGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTccacca                                                                   <  1:5148390/150‑1 (MQ=255)
                                                                                    aTGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTccacca                                                                   <  2:1930230/150‑1 (MQ=255)
                                                                                     tggctggcCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAg                                                                  >  1:588780/1‑150 (MQ=255)
                                                                                     tggctggcCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGACATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAg                                                                  <  1:1765993/150‑1 (MQ=255)
                                                                                         tggcCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCg                                                              >  2:366284/1‑150 (MQ=255)
                                                                                         tggcCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCg                                                              >  2:4661047/1‑150 (MQ=255)
                                                                                         tggcCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCg                                                              <  2:2318800/150‑1 (MQ=255)
                                                                                           gcCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCt                                                            <  2:2912720/150‑1 (MQ=255)
                                                                                              ggTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTtg                                                         >  1:4704575/1‑150 (MQ=255)
                                                                                                 tATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgc                                                      <  1:3600621/150‑1 (MQ=255)
                                                                                                  atcatcGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgcc                                                     <  2:1448221/150‑1 (MQ=255)
                                                                                                  atcatcGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCGTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgcc                                                     >  2:1595999/1‑150 (MQ=255)
                                                                                                   tcatcGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgccg                                                    <  1:998311/150‑1 (MQ=255)
                                                                                                   tcatcGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgccg                                                    >  1:628327/1‑150 (MQ=255)
                                                                                                    catcGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgccgc                                                   <  1:1646667/150‑1 (MQ=255)
                                                                                                     atcGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgccgcc                                                  >  1:1319912/1‑150 (MQ=255)
                                                                                                     atcGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTgccgcc                                                  >  2:4607373/1‑150 (MQ=255)
                                                                                                         aCAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCa                                              >  1:119519/1‑150 (MQ=255)
                                                                                                           aGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCAcc                                            <  1:4887052/150‑1 (MQ=255)
                                                                                                               tCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCAt                                        <  1:3069194/150‑1 (MQ=255)
                                                                                                                 aGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCa                                      >  2:3703971/1‑150 (MQ=255)
                                                                                                                 aGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCa                                      <  1:2730235/150‑1 (MQ=255)
                                                                                                                   cgcCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAat                                    >  1:2541618/1‑150 (MQ=255)
                                                                                                                    gcCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAata                                   <  2:4400244/150‑1 (MQ=255)
                                                                                                                     ccGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAatat                                  >  2:2406355/1‑150 (MQ=255)
                                                                                                                      cGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATc                                 >  1:2717256/1‑150 (MQ=255)
                                                                                                                       gAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCt                                <  1:3167166/150‑1 (MQ=255)
                                                                                                                        aGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCtt                               <  1:3245392/150‑1 (MQ=255)
                                                                                                                           gcAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTtat                            <  1:536829/150‑1 (MQ=255)
                                                                                                                           gcAGGATCTGCCCGCCACCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTtat                            >  1:1885824/1‑150 (MQ=255)
                                                                                                                             aGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATAc                          >  2:2976226/1‑150 (MQ=255)
                                                                                                                                  cTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCAt                     >  2:4918648/1‑150 (MQ=255)
                                                                                                                                  cTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCAt                     <  1:3812607/150‑1 (MQ=255)
                                                                                                                                   tGCCCGCCACCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATc                    <  1:3130803/150‑1 (MQ=255)
                                                                                                                                       cGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGa                <  2:2668585/150‑1 (MQ=255)
                                                                                                                                       cGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGa                >  1:2293937/1‑150 (MQ=255)
                                                                                                                                        gCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCgggct               >  2:3968021/1‑148 (MQ=255)
                                                                                                                                           aCCTTCGCCCAGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTtc            <  2:232535/150‑1 (MQ=255)
                                                                                                                                             cTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTtcgt          <  1:742845/150‑1 (MQ=255)
                                                                                                                                              ttCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTtcgtc         <  1:209965/150‑1 (MQ=255)
                                                                                                                                               tCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTtcgtcg        >  2:5081891/1‑150 (MQ=255)
                                                                                                                                                cGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCAACGGGATTTCGTCGa       >  2:3311395/1‑150 (MQ=255)
                                                                                                                                                   ccTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTTCGTCGATCa    <  2:418570/150‑1 (MQ=255)
                                                                                                                                                   ccTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTTCGTCGATCa    >  1:769411/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCCGTCGCCGCCTTTACCGCGGTCAGATGCTGGCGCAACAGCCCCGGTTTCGCCCGCCCGGCACGAATGCTGGTGATGGTAAATGGCTGGCCGGTTATCATCGACAGGCTCAGCGCCGAGCGCAGGATCTGCCCGCCACCTTCGCCCTGTGCGCCATCCAGCGCAATCATCCTTTTCATCATTTATCCTTTTACGCACACCACCTGACGCAGGGTATAGATAACTTCCACCAGATCGCTTTGTGCCGCCATCACCGCATCAATATCTTTATACGCCATCGGGATTTCGTCGATCACT  >  CP009273/3550028‑3550326

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: