Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,904,553 T→A 39.5% P111P (CCA→CCT pstS ← periplasmic phosphate binding protein, high‑affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,904,5530TA39.5% 214.6 / 313.0 248P111P (CCA→CCTpstSperiplasmic phosphate binding protein, high‑affinity
Reads supporting (aligned to +/- strand):  ref base T (89/61);  new base A (63/35);  total (152/96)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

GCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAAC  >  CP009273/3904406‑3904701
                                                                                                                                                   |                                                                                                                                                    
gcggcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTg                                                                                                                                                     <  1:727634/149‑1 (MQ=255)
     gTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                                >  1:4382109/1‑149 (MQ=255)
      tactacAGCACTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGt                                                                                                                                               >  2:548486/1‑149 (MQ=255)
          acAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAc                                                                                                                                           <  2:4225903/149‑1 (MQ=255)
          acAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAAc                                                                                                                                           <  1:1661771/149‑1 (MQ=255)
           cAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAcc                                                                                                                                          <  1:496061/149‑1 (MQ=255)
             gCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGc                                                                                                                                        >  1:2728750/1‑149 (MQ=255)
             gCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTAAGCCCAGGAATGTTAACCGc                                                                                                                                        >  2:2598078/1‑149 (MQ=255)
              cAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGcc                                                                                                                                       <  2:426054/149‑1 (MQ=255)
               aaTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCa                                                                                                                                      <  2:763685/149‑1 (MQ=255)
               aaTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCa                                                                                                                                      >  1:1874746/1‑149 (MQ=255)
                aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAg                                                                                                                                     >  2:2419223/1‑149 (MQ=255)
                 tGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGc                                                                                                                                    >  2:1892773/1‑149 (MQ=255)
                   ttttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGcac                                                                                                                                  >  1:1955108/1‑149 (MQ=255)
                   ttttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGcac                                                                                                                                  <  1:4333025/149‑1 (MQ=255)
                           aGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgccgc                                                                                                                          >  2:3326471/1‑149 (MQ=255)
                            ggCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACgccgcc                                                                                                                         >  1:89292/1‑149 (MQ=255)
                             gCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCa                                                                                                                        >  1:1471562/1‑149 (MQ=255)
                                gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATc                                                                                                                     <  2:2221808/149‑1 (MQ=255)
                                gTTTCAGACCCGGATTCAGTTTGGCGATGGCCTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATc                                                                                                                     <  2:4619409/149‑1 (MQ=255)
                                  ttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCAc                                                                                                                   <  1:2449484/149‑1 (MQ=255)
                                    cAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACgg                                                                                                                 <  1:2140277/149‑1 (MQ=255)
                                       aCCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCg                                                                                                              >  1:1355103/1‑149 (MQ=255)
                                       aCCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCg                                                                                                              >  1:4145433/1‑149 (MQ=255)
                                         ccGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCggg                                                                                                            <  2:164414/149‑1 (MQ=255)
                                           ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGTGTTTTACCATCCAGCACCAGTTCGCCAGTCTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGaa                                                                                                          <  2:3665375/149‑1 (MQ=255)
                                           ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGaa                                                                                                          >  1:2242810/1‑149 (MQ=255)
                                            gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                                                                                         >  1:2116021/1‑149 (MQ=255)
                                            gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                                                                                         >  2:2503453/1‑149 (MQ=255)
                                             aTTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACt                                                                                                        >  2:1830336/1‑149 (MQ=255)
                                               tCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTgg                                                                                                      >  2:643789/1‑149 (MQ=255)
                                                  gTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                                                                                   >  2:4734552/1‑149 (MQ=255)
                                                   tttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACa                                                                                                  <  1:1590587/149‑1 (MQ=255)
                                                   tttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACa                                                                                                  <  1:3888802/149‑1 (MQ=255)
                                                      ggCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGATTGTTAACCGCCAGCTCCACGCCGCCAATCACGGTCGGGAACTGGAACAGAc                                                                                               <  1:1186162/149‑1 (MQ=255)
                                                      ggCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAc                                                                                               <  2:4133302/149‑1 (MQ=255)
                                                           tGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTcc                                                                                          <  1:2726627/149‑1 (MQ=255)
                                                             gCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTg                                                                                        >  1:2365961/1‑149 (MQ=255)
                                                               ttCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                      <  2:1499902/149‑1 (MQ=255)
                                                                tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGACCAGACCTTCCTGAGc                                                                                     >  1:2855217/1‑149 (MQ=255)
                                                                tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                                                                                     <  1:3984048/149‑1 (MQ=255)
                                                                tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                                                                                     <  1:2302234/149‑1 (MQ=255)
                                                                 catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGcc                                                                                    >  2:2530430/1‑149 (MQ=255)
                                                                 catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGcc                                                                                    >  2:2530418/1‑149 (MQ=255)
                                                                  atcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                                                                                   >  2:3536570/1‑149 (MQ=255)
                                                                  atcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                                                                                   <  2:3985265/149‑1 (MQ=255)
                                                                   tcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAg                                                                                  <  1:3350507/149‑1 (MQ=255)
                                                                     atcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGtt                                                                                >  2:622533/1‑149 (MQ=255)
                                                                       cccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGACCTGGAACAGACCTTCCTGAGCCAGtttt                                                                              >  2:4556135/1‑149 (MQ=255)
                                                                       cccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGtttt                                                                              <  1:821073/149‑1 (MQ=255)
                                                                          aCTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCg                                                                           >  2:1190663/1‑149 (MQ=255)
                                                                               ttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa                                                                      >  2:2769014/1‑149 (MQ=255)
                                                                                 gATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACa                                                                    >  1:1907477/1‑149 (MQ=255)
                                                                                  aTTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg                                                                   >  1:1835073/1‑149 (MQ=255)
                                                                                  aTTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg                                                                   >  1:2568626/1‑149 (MQ=255)
                                                                                          cAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGt                                                           <  1:3493140/149‑1 (MQ=255)
                                                                                            ggTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                         >  1:3820668/1‑149 (MQ=255)
                                                                                            ggTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                         >  2:1293862/1‑149 (MQ=255)
                                                                                            ggTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                         >  2:1286605/1‑149 (MQ=255)
                                                                                             gTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCag                                                        >  2:1448502/1‑149 (MQ=255)
                                                                                                gATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg                                                     >  1:2655665/1‑149 (MQ=255)
                                                                                                gATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg                                                     >  1:3410854/1‑149 (MQ=255)
                                                                                                  tGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa                                                   >  2:1554060/1‑149 (MQ=255)
                                                                                                  tGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa                                                   <  1:2922589/149‑1 (MQ=255)
                                                                                                   gTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCATGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCAc                                                  >  2:3161497/1‑149 (MQ=255)
                                                                                                        cGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCaaaa                                             >  1:525844/1‑149 (MQ=255)
                                                                                                         gAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAAt                                            >  2:1989968/1‑149 (MQ=255)
                                                                                                         gAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAAt                                            >  2:348202/1‑149 (MQ=255)
                                                                                                         gAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAAt                                            >  1:147623/1‑149 (MQ=255)
                                                                                                             gTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAAc                                        <  2:3372150/149‑1 (MQ=255)
                                                                                                              ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACg                                       <  2:4482793/149‑1 (MQ=255)
                                                                                                              ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACg                                       >  1:124644/1‑149 (MQ=255)
                                                                                                                 tACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTa                                    <  1:3542222/149‑1 (MQ=255)
                                                                                                                   ccATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTAtt                                  <  1:1224254/149‑1 (MQ=255)
                                                                                                                   ccATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTAtt                                  <  1:3091519/149‑1 (MQ=255)
                                                                                                                   ccATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTAtt                                  <  1:1691196/149‑1 (MQ=255)
                                                                                                                   ccATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTAtt                                  >  1:4389961/1‑149 (MQ=255)
                                                                                                                     aTCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAg                                <  2:1516785/149‑1 (MQ=255)
                                                                                                                     aTCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAg                                >  2:2913244/1‑149 (MQ=255)
                                                                                                                      tCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGc                               >  1:3233843/1‑149 (MQ=255)
                                                                                                                        cAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGa                             <  2:4324143/149‑1 (MQ=255)
                                                                                                                          gCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATa                           <  2:2335668/149‑1 (MQ=255)
                                                                                                                          gCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATa                           >  2:1385202/1‑149 (MQ=255)
                                                                                                                            aCCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAAt                         >  2:2562561/1‑149 (MQ=255)
                                                                                                                             ccAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATc                        <  1:1034732/149‑1 (MQ=255)
                                                                                                                             ccAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATc                        >  1:466374/1‑149 (MQ=255)
                                                                                                                                gTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGt                     <  2:2923178/149‑1 (MQ=255)
                                                                                                                                  tCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGttt                   <  1:4245337/149‑1 (MQ=255)
                                                                                                                                   cGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGtttt                  >  2:3730372/1‑149 (MQ=255)
                                                                                                                                   cGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGtttt                  >  1:4547636/1‑149 (MQ=255)
                                                                                                                                       aGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGc              >  2:2014904/1‑149 (MQ=255)
                                                                                                                                        gACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc             <  2:901324/149‑1 (MQ=255)
                                                                                                                                        gACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc             <  1:1045583/149‑1 (MQ=255)
                                                                                                                                        gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc             <  1:643330/149‑1 (MQ=255)
                                                                                                                                         aCTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACGGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCa            >  2:444845/1‑149 (MQ=255)
                                                                                                                                          cTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCAc           >  2:4636159/1‑149 (MQ=255)
                                                                                                                                          cTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCAc           >  1:4750496/1‑149 (MQ=255)
                                                                                                                                           ttCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCAcc          >  2:2777859/1‑149 (MQ=255)
                                                                                                                                             cAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCgg        >  1:4463895/1‑149 (MQ=255)
                                                                                                                                               gCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGgaa      >  2:401373/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAGAAC  >  CP009273/3904406‑3904701

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: