Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,904,487 G→T 21.7% 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,4870GT21.7% 761.8 / 170.4 290I133I (ATC→ATApstSperiplasmic phosphate binding protein, high‑affinity
Reads supporting (aligned to +/- strand):  ref base G (117/110);  new base T (29/34);  total (146/144)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 295 total aligned reads displayed.

TTTTTCCACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAG  >  CP009273/3904338‑3904636
                                                                                                                                                     |                                                                                                                                                     
  tttCCACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGAt                                                                                                                                                     <  1:2145988/150‑1 (MQ=255)
   ttCCACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGAtt                                                                                                                                                    <  1:3505460/150‑1 (MQ=255)
              gTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTa                                                                                                                                         <  2:2384389/150‑1 (MQ=255)
              gTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTGTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTa                                                                                                                                         >  2:4009886/1‑150 (MQ=255)
                tCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGa                                                                                                                                       >  1:4598443/1‑150 (MQ=255)
                tCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGa                                                                                                                                       >  1:1667456/1‑150 (MQ=255)
                 cACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCCTCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCCGGTAGAt                                                                                                                                      >  1:4503775/1‑150 (MQ=255)
                 cACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGAt                                                                                                                                      <  1:3744844/150‑1 (MQ=255)
                  aCTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATg                                                                                                                                     <  2:696091/150‑1 (MQ=255)
                     ttCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCg                                                                                                                                  >  1:1993577/1‑150 (MQ=255)
                     ttCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCg                                                                                                                                  <  1:255594/150‑1 (MQ=255)
                     ttCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCg                                                                                                                                  <  2:2739084/150‑1 (MQ=255)
                     ttCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCg                                                                                                                                  <  2:3916606/150‑1 (MQ=255)
                     ttCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCg                                                                                                                                  >  1:2873941/1‑150 (MQ=255)
                      tCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGc                                                                                                                                 <  2:3289575/150‑1 (MQ=255)
                           aGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAgg                                                                                                                            <  1:3004198/150‑1 (MQ=255)
                                 cTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTa                                                                                                                      >  2:2128086/1‑150 (MQ=255)
                                 cTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTa                                                                                                                      >  2:25806/1‑150 (MQ=255)
                                 cTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTa                                                                                                                      <  1:1262514/150‑1 (MQ=255)
                                  tGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTAc                                                                                                                     <  2:4164667/150‑1 (MQ=255)
                                   ggTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTAcc                                                                                                                    <  1:1722806/150‑1 (MQ=255)
                                   ggTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTAcc                                                                                                                    <  2:73220/150‑1 (MQ=255)
                                    gTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCa                                                                                                                   <  1:3877512/150‑1 (MQ=255)
                                    gTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCa                                                                                                                   >  1:3084694/1‑150 (MQ=255)
                                      gAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATc                                                                                                                 >  1:2724621/1‑150 (MQ=255)
                                       aaGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATcc                                                                                                                <  1:1205384/150‑1 (MQ=255)
                                           cGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCa                                                                                                            <  1:2352760/150‑1 (MQ=255)
                                             aaggaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCAcc                                                                                                          <  2:1297505/150‑1 (MQ=255)
                                               ggaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAg                                                                                                        <  1:1506919/150‑1 (MQ=255)
                                                gaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGt                                                                                                       <  2:474737/150‑1 (MQ=255)
                                                gaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGt                                                                                                       <  1:3526268/150‑1 (MQ=255)
                                                gaagTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGt                                                                                                       <  2:900664/150‑1 (MQ=255)
                                                    tCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGc                                                                                                   >  1:1101656/1‑150 (MQ=255)
                                                    tCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGc                                                                                                   >  1:3382446/1‑150 (MQ=255)
                                                        ggAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGa                                                                                               >  1:874204/1‑150 (MQ=255)
                                                           gCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACtt                                                                                            <  1:1429185/150‑1 (MQ=255)
                                                            ccATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTc                                                                                           <  2:1667456/150‑1 (MQ=255)
                                                            ccATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTc                                                                                           <  2:3597184/150‑1 (MQ=255)
                                                             cATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCa                                                                                          <  2:3822238/150‑1 (MQ=255)
                                                              aTCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAg                                                                                         <  2:3760945/150‑1 (MQ=255)
                                                              aTCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAg                                                                                         >  1:4165738/1‑150 (MQ=255)
                                                                cTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGcc                                                                                       >  1:2602937/1‑150 (MQ=255)
                                                                   cgcgGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTg                                                                                    >  1:1529/1‑150 (MQ=255)
                                                                   cgcgGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTg                                                                                    >  1:1584183/1‑150 (MQ=255)
                                                                   cgcgGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTg                                                                                    >  2:785097/1‑150 (MQ=255)
                                                                       gcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAAt                                                                                <  2:3045044/150‑1 (MQ=255)
                                                                        cgTACTACAGCAATGTTTTGTGAAGGCCGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                               >  2:1549565/1‑150 (MQ=255)
                                                                        cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATg                                                                               <  2:4377615/150‑1 (MQ=255)
                                                                           actacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTa                                                                            <  2:629435/150‑1 (MQ=255)
                                                                           actacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTa                                                                            >  1:2515699/1‑150 (MQ=255)
                                                                             tacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAc                                                                          >  1:3427692/1‑150 (MQ=255)
                                                                              acAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTTCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAcc                                                                         >  2:1699639/1‑150 (MQ=255)
                                                                              acAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAAcc                                                                         <  1:1003392/150‑1 (MQ=255)
                                                                                    aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGc                                                                   >  1:3614488/1‑150 (MQ=255)
                                                                                        tttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGcacca                                                               <  2:4056374/150‑1 (MQ=255)
                                                                                        tttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGcacca                                                               <  2:4405699/150‑1 (MQ=255)
                                                                                        tttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGcacca                                                               >  2:2868637/1‑150 (MQ=255)
                                                                                            tgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgcc                                                           >  2:3659902/1‑150 (MQ=255)
                                                                                            tgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgcc                                                           <  1:545722/150‑1 (MQ=255)
                                                                                            tgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgcc                                                           <  1:1817117/150‑1 (MQ=255)
                                                                                              aaGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgccgc                                                         >  2:1318631/1‑150 (MQ=255)
                                                                                              aaGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCAcgccgc                                                         >  1:2796076/1‑150 (MQ=255)
                                                                                               aGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACgccgcc                                                        >  1:4232361/1‑150 (MQ=255)
                                                                                                   aGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATc                                                    >  2:2349216/1‑150 (MQ=255)
                                                                                                    gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                   <  2:2254486/150‑1 (MQ=255)
                                                                                                    gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCAATTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                   <  1:2894585/150‑1 (MQ=255)
                                                                                                      ttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACg                                                 <  2:477426/150‑1 (MQ=255)
                                                                                                      ttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACg                                                 <  2:1527199/150‑1 (MQ=255)
                                                                                                       tCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACgg                                                >  1:4005336/1‑150 (MQ=255)
                                                                                                        cAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGt                                               >  1:4353170/1‑150 (MQ=255)
                                                                                                            cccGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCggg                                           >  2:3192899/1‑150 (MQ=255)
                                                                                                               ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                        >  2:1772205/1‑150 (MQ=255)
                                                                                                                gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACt                                       >  1:925760/1‑150 (MQ=255)
                                                                                                                 aTTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTg                                      <  1:761428/150‑1 (MQ=255)
                                                                                                                     aGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                  <  2:349474/150‑1 (MQ=255)
                                                                                                                        ttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGa                               <  2:4587820/150‑1 (MQ=255)
                                                                                                                        ttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGa                               <  1:504812/150‑1 (MQ=255)
                                                                                                                         tGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAc                              >  2:2693340/1‑150 (MQ=255)
                                                                                                                         tGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAc                              >  2:3329441/1‑150 (MQ=255)
                                                                                                                           gCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCt                            <  1:1265782/150‑1 (MQ=255)
                                                                                                                            cGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCtt                           >  1:2443206/1‑150 (MQ=255)
                                                                                                                             gATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTc                          >  2:1224154/1‑150 (MQ=255)
                                                                                                                              aTGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTcc                         >  1:3591307/1‑150 (MQ=255)
                                                                                                                               tGGCTTCATCATCCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAc                              >  1:1385640/1‑144 (MQ=255)
                                                                                                                                  cTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                     >  1:4054983/1‑150 (MQ=255)
                                                                                                                                  cTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                     >  1:4542070/1‑150 (MQ=255)
                                                                                                                                   ttCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                    >  2:1138210/1‑150 (MQ=255)
                                                                                                                                   ttCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                    >  2:1397691/1‑150 (MQ=255)
                                                                                                                                   ttCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                    >  2:511466/1‑150 (MQ=255)
                                                                                                                                     catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                  >  1:1561820/1‑150 (MQ=255)
                                                                                                                                       tcatcCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGt                <  2:1374776/150‑1 (MQ=255)
                                                                                                                                        catcCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGtt               >  1:3330824/1‑150 (MQ=255)
                                                                                                                                         atcCCACTTCTTTATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttt              <  2:1965044/150‑1 (MQ=255)
                                                                                                                                             cACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGTATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCg          <  1:4194540/150‑1 (MQ=255)
                                                                                                                                                ttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCa       <  2:1869483/150‑1 (MQ=255)
                                                                                                                                                ttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCa       >  1:1339971/1‑150 (MQ=255)
                                                                                                                                                  cttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa     <  2:2796076/150‑1 (MQ=255)
                                                                                                                                                     gATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCATGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg  <  2:2014790/150‑1 (MQ=255)
                                                                                                                                                     gATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCATGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg  <  2:2485312/150‑1 (MQ=255)
                                                                                                                                                     gATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg  >  1:1713804/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTTTCCACTCTTCGTTCACTTTCGCCAGGTAGCTGGTGAAGACGAAGGAAGTCCCGGAGCCATCTGCGCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAG  >  CP009273/3904338‑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: