Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,904,487 G→T 23.4% I133I (ATC→ATA 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,4870GT23.4% 538.7 / 136.8 220I133I (ATC→ATApstSperiplasmic phosphate binding protein, high‑affinity
Reads supporting (aligned to +/- strand):  ref base G (86/82);  new base T (16/36);  total (102/118)
Fisher's exact test for biased strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TTTTCCACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAG  >  CP009273/3904339‑3904636
                                                                                                                                                    |                                                                                                                                                     
   tCCACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTAttt                                                                                                                                                   <  1:2620744/150‑1 (MQ=255)
     cACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTg                                                                                                                                                 <  1:923840/150‑1 (MQ=255)
     cACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTg                                                                                                                                                 >  1:3798422/1‑150 (MQ=255)
     cACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTg                                                                                                                                                 >  2:1543931/1‑150 (MQ=255)
            cGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGt                                                                                                                                          >  2:3125351/1‑150 (MQ=255)
            cGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGt                                                                                                                                          <  1:1386609/150‑1 (MQ=255)
            cGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGAGTTTGCCCAGGt                                                                                                                                          >  2:3561650/1‑150 (MQ=255)
             gTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTa                                                                                                                                         <  1:436957/150‑1 (MQ=255)
             gTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTa                                                                                                                                         <  1:1856357/150‑1 (MQ=255)
                cACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGAt                                                                                                                                      <  1:3933965/150‑1 (MQ=255)
                cACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGAt                                                                                                                                      <  1:415728/150‑1 (MQ=255)
                cACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGAt                                                                                                                                      <  1:4318322/150‑1 (MQ=255)
                    ttCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCg                                                                                                                                  <  1:880868/150‑1 (MQ=255)
                     tCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGc                                                                                                                                 <  2:1714974/150‑1 (MQ=255)
                      cGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGcc                                                                                                                                <  2:2349248/150‑1 (MQ=255)
                      cGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGcc                                                                                                                                <  1:4123505/150‑1 (MQ=255)
                          aGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAgg                                                                                                                            <  1:4811240/150‑1 (MQ=255)
                           ggTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAggg                                                                                                                           >  2:4164261/1‑150 (MQ=255)
                                       aGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCa                                                                                                               <  1:4804995/150‑1 (MQ=255)
                                       aGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCa                                                                                                               <  1:1238596/150‑1 (MQ=255)
                                        gACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAg                                                                                                              >  2:562583/1‑150 (MQ=255)
                                          cGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCa                                                                                                            >  2:4012721/1‑150 (MQ=255)
                                            aaggaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCAcc                                                                                                          >  1:307834/1‑150 (MQ=255)
                                            aaggaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCAcc                                                                                                          >  2:1692727/1‑150 (MQ=255)
                                             aggaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCa                                                                                                         <  2:3116663/150‑1 (MQ=255)
                                             aggaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCa                                                                                                         >  1:1553858/1‑150 (MQ=255)
                                                  gTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCg                                                                                                    <  1:4553390/150‑1 (MQ=255)
                                                  gTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCg                                                                                                    <  2:2532786/150‑1 (MQ=255)
                                                     ccGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCa                                                                                                 <  2:4348587/150‑1 (MQ=255)
                                                      cGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAg                                                                                                >  2:3684842/1‑150 (MQ=255)
                                                      cGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAg                                                                                                >  2:250812/1‑150 (MQ=255)
                                                      cGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAg                                                                                                >  1:2002453/1‑150 (MQ=255)
                                                       ggAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGa                                                                                               >  2:4505456/1‑150 (MQ=255)
                                                           ccATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTc                                                                                           <  1:4251679/150‑1 (MQ=255)
                                                           ccATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTc                                                                                           <  1:3732385/150‑1 (MQ=255)
                                                            cATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCa                                                                                          <  1:3326527/150‑1 (MQ=255)
                                                               cTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGcc                                                                                       >  2:615828/1‑150 (MQ=255)
                                                                 gcgcGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCt                                                                                     <  2:2790051/150‑1 (MQ=255)
                                                                 gcgcGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCt                                                                                     <  2:4993365/150‑1 (MQ=255)
                                                                     ggcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGaa                                                                                 >  2:687065/1‑150 (MQ=255)
                                                                      gcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAAt                                                                                <  2:3775640/150‑1 (MQ=255)
                                                                       cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                               <  1:562583/150‑1 (MQ=255)
                                                                       cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                               <  1:3487353/150‑1 (MQ=255)
                                                                       cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                               >  2:4795146/1‑150 (MQ=255)
                                                                       cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                               >  2:288764/1‑150 (MQ=255)
                                                                        gTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGt                                                                              >  2:1855342/1‑150 (MQ=255)
                                                                        gTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGt                                                                              <  2:3857173/150‑1 (MQ=255)
                                                                        gTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGt                                                                              >  1:2417257/1‑150 (MQ=255)
                                                                         tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGtt                                                                             >  1:2055158/1‑150 (MQ=255)
                                                                         tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTACTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGtt                                                                             <  1:3936939/150‑1 (MQ=255)
                                                                           ctacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTaa                                                                           >  2:1223553/1‑150 (MQ=255)
                                                                            tacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAc                                                                          <  1:1080533/150‑1 (MQ=255)
                                                                                  aaTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAg                                                                    <  2:5172858/150‑1 (MQ=255)
                                                                                  aaTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAg                                                                    <  1:4381281/150‑1 (MQ=255)
                                                                                  aaTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAg                                                                    >  2:4743740/1‑150 (MQ=255)
                                                                                   aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGc                                                                   <  1:2513535/150‑1 (MQ=255)
                                                                                   aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGc                                                                   >  1:2103907/1‑150 (MQ=255)
                                                                                      ttttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGcacc                                                                >  2:4310138/1‑150 (MQ=255)
                                                                                        ttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGcaccac                                                              >  1:1471020/1‑150 (MQ=255)
                                                                                         tgtgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcg                                                             >  2:3323768/1‑150 (MQ=255)
                                                                                         tgtgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcg                                                             <  1:4027005/150‑1 (MQ=255)
                                                                                           tgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCATCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTATCCGCCAGCACCAcgcc                                                           <  1:1754979/150‑1 (MQ=255)
                                                                                           tgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgcc                                                           <  2:840501/150‑1 (MQ=255)
                                                                                            gAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgccg                                                          <  2:4987942/150‑1 (MQ=255)
                                                                                               ggCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCa                                                       >  1:2292426/1‑150 (MQ=255)
                                                                                                gCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCaa                                                      <  2:1412981/150‑1 (MQ=255)
                                                                                                   gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                   >  2:4606421/1‑150 (MQ=255)
                                                                                                    tttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCAc                                                  <  2:293847/150‑1 (MQ=255)
                                                                                                     ttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACg                                                 >  2:4526641/1‑150 (MQ=255)
                                                                                                       cAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGt                                               >  2:2908433/1‑150 (MQ=255)
                                                                                                         gACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCg                                             <  1:3018869/150‑1 (MQ=255)
                                                                                                          aCCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCgg                                            <  1:3398538/150‑1 (MQ=255)
                                                                                                              ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                        >  1:2385238/1‑150 (MQ=255)
                                                                                                              ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                        <  2:4683396/150‑1 (MQ=255)
                                                                                                              ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                        <  2:2124304/150‑1 (MQ=255)
                                                                                                               gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACt                                       <  2:3463665/150‑1 (MQ=255)
                                                                                                                aTTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTg                                      <  2:4376514/150‑1 (MQ=255)
                                                                                                                aTTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTg                                      >  1:1586497/1‑150 (MQ=255)
                                                                                                                aTTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTg                                      >  1:1586652/1‑150 (MQ=255)
                                                                                                                   cAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGGAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGaa                                   <  2:5052495/150‑1 (MQ=255)
                                                                                                                    aGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                  >  1:521635/1‑150 (MQ=255)
                                                                                                                    aGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                  >  1:3001820/1‑150 (MQ=255)
                                                                                                                       ttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGa                               <  1:1552149/150‑1 (MQ=255)
                                                                                                                       ttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGa                               <  1:3162152/150‑1 (MQ=255)
                                                                                                                               ggCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTg                       <  2:157778/150‑1 (MQ=255)
                                                                                                                                 cTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                     <  1:3133277/150‑1 (MQ=255)
                                                                                                                                    catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGACTGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                  >  1:4669368/1‑150 (MQ=255)
                                                                                                                                    catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                  >  1:2882263/1‑150 (MQ=255)
                                                                                                                                      tcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGt                >  2:4907318/1‑150 (MQ=255)
                                                                                                                                        atcCCACTTCGTGATTTGGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttt              <  2:5075593/150‑1 (MQ=255)
                                                                                                                                          cccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttttt            <  2:2685405/150‑1 (MQ=255)
                                                                                                                                           ccACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTc           <  1:3655851/150‑1 (MQ=255)
                                                                                                                                            cACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTct          <  2:492384/150‑2 (MQ=255)
                                                                                                                                             aCTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGt         >  1:5071278/1‑150 (MQ=255)
                                                                                                                                             aCTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGt         <  2:1438747/150‑1 (MQ=255)
                                                                                                                                              cttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTc        >  1:3199031/1‑150 (MQ=255)
                                                                                                                                                tcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAg      >  1:1115748/1‑150 (MQ=255)
                                                                                                                                                 cttTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa     >  2:359637/1‑150 (MQ=255)
                                                                                                                                                 cttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa     <  2:188163/150‑1 (MQ=255)
                                                                                                                                                  ttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGAc    <  2:2917316/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTTCCACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAG  >  CP009273/3904339‑3904636

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: