Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,906,096 T→A 6.7% K246N (AAA→AAT htpX ← putative endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,906,0960TA 6.7% 894.8 / 41.1 248K246N (AAA→AAThtpXputative endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (129/102);  new base A (11/6);  total (140/108)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 257 total aligned reads displayed.

CATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACC  >  CP009273/1905948‑1906243
                                                                                                                                                    |                                                                                                                                                   
  tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGtttt                                                                                                                                                   <  1:236092/149‑1 (MQ=255)
  tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGtttt                                                                                                                                                   <  2:360098/149‑1 (MQ=255)
  tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGtttt                                                                                                                                                   <  2:2030871/149‑1 (MQ=255)
    ttACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa                                                                                                                                                 <  2:3265137/149‑1 (MQ=255)
    ttACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa                                                                                                                                                 <  1:5035313/149‑1 (MQ=255)
     tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAg                                                                                                                                                <  1:3494160/149‑1 (MQ=255)
        ttCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGgcg                                                                                                                                             >  1:1532062/1‑149 (MQ=255)
            ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGc                                                                                                                                         >  1:1259407/1‑149 (MQ=255)
                ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAgcg                                                                                                                                     >  1:2773616/1‑149 (MQ=255)
                ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAgcg                                                                                                                                     >  1:810602/1‑149 (MQ=255)
                    cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGc                                                                                                                                 >  2:5192297/1‑149 (MQ=255)
                      cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCaa                                                                                                                               <  1:943465/149‑1 (MQ=255)
                       gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAAt                                                                                                                              <  1:4678757/149‑1 (MQ=255)
                        tACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATc                                                                                                                             >  2:4338405/1‑149 (MQ=255)
                           gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAtt                                                                                                                          >  1:275739/1‑149 (MQ=255)
                             agagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAtttt                                                                                                                        >  2:3712049/1‑149 (MQ=255)
                              gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtc                                                                                                                       >  2:902727/1‑149 (MQ=255)
                                gCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtctc                                                                                                                     <  1:1323488/149‑1 (MQ=255)
                                 cTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcg                                                                                                                    >  1:2430855/1‑149 (MQ=255)
                                  ttCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcgc                                                                                                                   <  1:3617413/149‑1 (MQ=255)
                                  ttCAATTCGTGTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACGTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcgc                                                                                                                   <  1:1830996/149‑1 (MQ=255)
                                    cAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGa                                                                                                                 <  1:1258252/149‑1 (MQ=255)
                                       ttCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGacca                                                                                                              <  2:1821550/149‑1 (MQ=255)
                                       ttCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGacca                                                                                                              >  1:1369778/1‑149 (MQ=255)
                                         cGTTTATCCAGCGGCGGGTGGGTCATGAACAACACACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCAAGCTGGTTGCTTCATGCGGTTCATAGCTGTTTTTCAGGCGCTGCAGCGCGGCAATCATATTCTCGCGaccaac                                                                                                            >  1:523205/1‑149 (MQ=255)
                                               tCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                      <  1:2532292/149‑1 (MQ=255)
                                               tCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                      >  2:2806741/1‑149 (MQ=255)
                                                ccAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGtttt                                                                                                     <  2:2696713/149‑1 (MQ=255)
                                                 cAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTg                                                                                                    >  2:4324484/1‑149 (MQ=255)
                                                         ggtgggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCg                                                                                            >  2:4715834/1‑149 (MQ=255)
                                                         ggtgggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCg                                                                                            <  1:1134727/149‑1 (MQ=255)
                                                            gggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCa                                                                                         <  1:1477250/149‑1 (MQ=255)
                                                              gtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATc                                                                                       <  1:2255337/149‑1 (MQ=255)
                                                                  tGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCa                                                                                   >  2:1207170/1‑149 (MQ=255)
                                                                    aacaacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATg                                                                                 >  1:1514404/1‑149 (MQ=255)
                                                                     acaacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATgg                                                                                <  2:4271690/149‑1 (MQ=255)
                                                                       aacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGaa                                                                              <  1:3082720/149‑1 (MQ=255)
                                                                        acTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAAt                                                                             <  1:4550921/149‑1 (MQ=255)
                                                                           cACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCa                                                                          >  2:5145449/1‑149 (MQ=255)
                                                                            aCTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCAc                                                                         <  2:2344909/149‑1 (MQ=255)
                                                                            aCTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCAc                                                                         <  2:5127605/149‑1 (MQ=255)
                                                                              tGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGa                                                                       >  2:1761561/1‑149 (MQ=255)
                                                                               gAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGAt                                                                      >  1:2187577/1‑149 (MQ=255)
                                                                                 gCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGATTCTTGCGGTTCATAGCTGGTTTTCAGGCTCTGCCGCGCGGCAATCATATTCTCGCGACCAACCAGTATTGCCGAACCGGCATAAGCATGGAATTCACGATGa                                                                    >  2:1695703/1‑149 (MQ=255)
                                                                                  cGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGAc                                                                   <  1:637401/149‑1 (MQ=255)
                                                                                   gATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGAcg                                                                  >  2:130105/1‑149 (MQ=255)
                                                                                        cGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCgaga                                                             >  2:2138265/1‑149 (MQ=255)
                                                                                         gACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGaa                                                            <  2:1057842/149‑1 (MQ=255)
                                                                                             tACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACcac                                                        <  1:4199551/149‑1 (MQ=255)
                                                                                               ccGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACAt                                                      <  2:1338054/149‑1 (MQ=255)
                                                                                                cGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATg                                                     >  1:2065751/1‑149 (MQ=255)
                                                                                                  ttAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGt                                                   >  1:510013/1‑149 (MQ=255)
                                                                                                   tAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTg                                                  <  1:1576096/149‑1 (MQ=255)
                                                                                                    aaTGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGa                                                 >  2:5001829/1‑149 (MQ=255)
                                                                                                    aaTGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGa                                                 <  1:709762/149‑1 (MQ=255)
                                                                                                    aaTGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGa                                                 <  2:633361/149‑1 (MQ=255)
                                                                                                     aTGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATTGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGAt                                                <  1:3056790/149‑1 (MQ=255)
                                                                                                        cAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATact                                             >  1:543357/1‑147 (MQ=255)
                                                                                                        cAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAAt                                             >  2:4949768/1‑149 (MQ=255)
                                                                                                         agagagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATg                                            >  1:4246448/1‑149 (MQ=255)
                                                                                                          gagagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGc                                           >  1:4908217/1‑149 (MQ=255)
                                                                                                           agagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCt                                          <  1:4589771/149‑1 (MQ=255)
                                                                                                           agagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCt                                          >  2:2879799/1‑149 (MQ=255)
                                                                                                            gagCCATCTTGCTGGTTGCTTCTTGCGGTTCATCGCTGGTTTTCAGGCGCTGCTGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCCTGGAATTCACGATGTCGCGAGAACCACATGGTGATAATGCTc                                         <  1:2623490/149‑1 (MQ=255)
                                                                                                            gagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTc                                         >  2:3217971/1‑149 (MQ=255)
                                                                                                             agCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCg                                        <  2:4805888/149‑1 (MQ=255)
                                                                                                              gCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGc                                       >  2:2940314/1‑149 (MQ=255)
                                                                                                              gCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGc                                       >  2:2663602/1‑149 (MQ=255)
                                                                                                                  tcatGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGa                                   >  1:4357370/1‑149 (MQ=255)
                                                                                                                  tcatGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGa                                   >  1:3241113/1‑149 (MQ=255)
                                                                                                                  tcatGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGa                                   >  1:3276860/1‑149 (MQ=255)
                                                                                                                    atGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                 >  1:1069836/1‑149 (MQ=255)
                                                                                                                    atGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                 >  1:1617502/1‑149 (MQ=255)
                                                                                                                    atGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                 >  2:118625/1‑149 (MQ=255)
                                                                                                                    atGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                 <  2:4115156/149‑1 (MQ=255)
                                                                                                                    atGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                 >  2:2169625/1‑149 (MQ=255)
                                                                                                                     tGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATa                                <  2:3902508/149‑1 (MQ=255)
                                                                                                                          gTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCaaa                           >  1:3841253/1‑149 (MQ=255)
                                                                                                                          gTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCaaa                           >  2:1495541/1‑149 (MQ=255)
                                                                                                                          gTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTATTGCCGAACCGGCATCAGCATGGAATTCACGATGAATCGAGAACCACATGGTGATAATGCTCGCCAGAATACCaaa                           >  1:1244987/1‑149 (MQ=255)
                                                                                                                           ttGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAac                          <  1:3681646/149‑1 (MQ=255)
                                                                                                                                tcttGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAg                     <  1:4861311/149‑1 (MQ=255)
                                                                                                                                  ttGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGtt                   <  1:4411283/149‑1 (MQ=255)
                                                                                                                                    gCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTcc                 <  2:668306/149‑1 (MQ=255)
                                                                                                                                       gTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGa              >  1:484903/1‑149 (MQ=255)
                                                                                                                                       gTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGa              >  2:3011485/1‑149 (MQ=255)
                                                                                                                                       gTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGa              >  1:1031494/1‑149 (MQ=255)
                                                                                                                                       gTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGATTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGa              >  2:2763159/1‑149 (MQ=255)
                                                                                                                                          cATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAAcc           >  2:2461264/1‑149 (MQ=255)
                                                                                                                                            tAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGt         <  1:3872524/149‑1 (MQ=255)
                                                                                                                                            tAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGt         <  2:2417121/149‑1 (MQ=255)
                                                                                                                                              gCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCCGAACCGTTg       >  2:4221195/1‑149 (MQ=255)
                                                                                                                                              gCTGGTATTCAGGCGCTGCAGCGCGGCAATCATATTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGACTACCAAACACCAGTTCCAGAACCGTTg       >  2:3469283/1‑149 (MQ=255)
                                                                                                                                               cTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGc      >  2:685337/1‑149 (MQ=255)
                                                                                                                                                 ggTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCaa    >  1:1909017/1‑149 (MQ=255)
                                                                                                                                                 ggTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCaa    >  2:4429345/1‑149 (MQ=255)
                                                                                                                                                 ggTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCaa    >  2:3807279/1‑149 (MQ=255)
                                                                                                                                                 ggTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCaa    <  1:1767562/149‑1 (MQ=255)
                                                                                                                                                 ggTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCaa    >  2:3068168/1‑149 (MQ=255)
                                                                                                                                                 ggTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCaa    <  2:753025/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACC  >  CP009273/1905948‑1906243

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: