Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,478,971 G→A 100% E552K (GAG→AAG)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,478,9710GA99.3% 1282.3 / ‑1.3 402E552K (GAG→AAG) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (201/196);  total (201/201)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 405 total aligned reads displayed.

CGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCCAG  >  CP009273/1478825‑1479117
                                                                                                                                                  |                                                                                                                                                  
  cgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCt                                                                                                                                               >  1:2381377/1‑150 (MQ=255)
   gTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTg                                                                                                                                              >  2:2789182/1‑150 (MQ=255)
   gTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTg                                                                                                                                              >  1:1420744/1‑150 (MQ=255)
    tATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGa                                                                                                                                             >  1:4065670/1‑150 (MQ=255)
     aTGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGAc                                                                                                                                            >  1:978551/1‑150 (MQ=255)
      tGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACt                                                                                                                                           <  1:347794/150‑1 (MQ=255)
       ggTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACtt                                                                                                                                          <  1:3762727/150‑1 (MQ=255)
         tGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtc                                                                                                                                        <  1:3178915/150‑1 (MQ=255)
         tGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtc                                                                                                                                        >  2:4866845/1‑150 (MQ=255)
          gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtct                                                                                                                                       <  1:1642852/150‑1 (MQ=255)
          gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtct                                                                                                                                       <  1:152276/150‑1 (MQ=255)
          gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtct                                                                                                                                       <  2:3658752/150‑1 (MQ=255)
          gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtct                                                                                                                                       <  2:2422309/150‑1 (MQ=255)
           cTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtctc                                                                                                                                      >  2:1965351/1‑150 (MQ=255)
              gAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                   >  2:1342010/1‑150 (MQ=255)
               aaGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGt                                                                                                                                  <  2:4643913/150‑1 (MQ=255)
                aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGtt                                                                                                                                 <  1:5317376/150‑1 (MQ=255)
                aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGtt                                                                                                                                 <  1:4398128/150‑1 (MQ=255)
                    cAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTAGCGTTtgtg                                                                                                                             <  1:1351449/150‑1 (MQ=255)
                       aaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAAt                                                                                                                          <  1:1824725/150‑1 (MQ=255)
                        aaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATc                                                                                                                         <  2:192049/150‑1 (MQ=255)
                          cATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtg                                                                                                                       >  1:517593/1‑150 (MQ=255)
                              gCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGa                                                                                                                   <  1:3350355/150‑1 (MQ=255)
                                   gtgcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTAt                                                                                                              <  2:1674026/150‑1 (MQ=255)
                                   gtgcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTAt                                                                                                              <  1:128103/150‑1 (MQ=255)
                                     gcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCt                                                                                                            <  1:3462672/150‑1 (MQ=255)
                                      cgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCtt                                                                                                           >  1:2587616/1‑150 (MQ=255)
                                      cgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCtt                                                                                                           <  2:1719871/150‑1 (MQ=255)
                                        tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTgg                                                                                                         >  1:2117828/1‑150 (MQ=255)
                                        tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTGTGTGAATCTGTGGAATTATCTTgg                                                                                                         >  2:422287/1‑150 (MQ=255)
                                          aGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCg                                                                                                       <  1:109883/150‑1 (MQ=255)
                                              gatgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGACTTATCTTGGCGagca                                                                                                   >  2:109797/1‑150 (MQ=255)
                                              gatgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                                                                                   >  2:3486863/1‑150 (MQ=255)
                                                tgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagcagc                                                                                                 <  2:4557009/150‑1 (MQ=255)
                                                  atTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTTTCTTGGCGAGCAGCaa                                                                                               >  1:3590009/1‑150 (MQ=255)
                                                   ttATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCaaa                                                                                              <  1:4162209/150‑1 (MQ=255)
                                                      tCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTAACTTGGCGAGCAGCAAAAgg                                                                                           <  1:4913067/150‑1 (MQ=255)
                                                           tCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCtt                                                                                      >  2:220829/1‑150 (MQ=255)
                                                             cgcgCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACAACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttc                                                                                    <  2:28363/150‑1 (MQ=255)
                                                               cgcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATATGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttctt                                                                                  <  2:1260398/150‑1 (MQ=255)
                                                               cgcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACAGCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGGGGGGGAATCTGTGGACGTATCTTGGCGAGCAGCAATAGGCGCTttctt                                                                                  >  2:857179/1‑150 (MQ=255)
                                                                gcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttcttc                                                                                 <  1:5257235/150‑1 (MQ=255)
                                                                gatatcaaTTCAGGTTCCGCGATTTAGTCCGATGGAAATATAGAAGGCATCGGAAGTTAAAATTCGTGTCTTCCAATATTATAATTCTGACTTTCTCGCGTTTGTGATTCTGTGGATTTATCTTGGCGTGCTGGATTAGGCGCTttcttc                                                                                 <  2:4451322/143‑1 (MQ=255)
                                                                      cATTCCGGCTCCGCGCGAACGTCCGATGGACAAACCGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGc                                                                           <  2:1649602/150‑1 (MQ=255)
                                                                      cATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGc                                                                           >  1:873059/1‑150 (MQ=255)
                                                                       aTTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGcc                                                                          >  2:4406892/1‑150 (MQ=255)
                                                                         tCAGGATCCGCTCGAACGTCCGATTGACAAACAACAGGAATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGAAAAAGGCGCTTTCTTCCAACGCCtt                                                                        <  1:3315504/150‑1 (MQ=255)
                                                                         tCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCGTTCTTCCAACGCCtt                                                                        <  2:1185310/150‑1 (MQ=255)
                                                                         tCAGGATCCGCGCGAACGTCCGATGGACAAACAACAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCtt                                                                        <  1:4972001/150‑1 (MQ=255)
                                                                          cAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTc                                                                       >  1:1125114/1‑150 (MQ=255)
                                                                             gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGt                                                                    >  2:4001749/1‑150 (MQ=255)
                                                                                 cgcgcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGcc                                                                >  1:2473596/1‑150 (MQ=255)
                                                                                   cgcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTATCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCtg                                                              <  2:4754319/150‑1 (MQ=255)
                                                                                    gcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCtgt                                                             <  1:1931418/150‑1 (MQ=255)
                                                                                     cgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGGTTGTGAATCTGTGGAAGTATCTTGGCGAGCAGCAAAAGGCGCTTTCGTCCAACGCCGTCCGTCGCCtgtg                                                            >  1:2042389/1‑150 (MQ=255)
                                                                                      gAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTgtgt                                                           >  2:3208867/1‑150 (MQ=255)
                                                                                       aaCGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTc                                                          <  2:2278211/150‑1 (MQ=255)
                                                                                       aaCGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTc                                                          <  2:4104829/150‑1 (MQ=255)
                                                                                       aaCGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGGGTc                                                          <  1:2291793/150‑1 (MQ=255)
                                                                                           tCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTAc                                                      <  2:2518965/150‑1 (MQ=255)
                                                                                            ccGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTAcc                                                     >  2:2340252/1‑150 (MQ=255)
                                                                                             cGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCg                                                    >  2:381911/1‑150 (MQ=255)
                                                                                               aTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGAt                                                  >  1:1820292/1‑150 (MQ=255)
                                                                                                 ggACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTa                                                <  1:520065/150‑1 (MQ=255)
                                                                                                 ggACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGGGAATCTGTGGAATTATCGTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGGGTCGTACCGAttc                                                >  1:604861/1‑149 (MQ=255)
                                                                                                  gACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAt                                               <  1:2181143/150‑1 (MQ=255)
                                                                                                   aCAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGGCGCCTGTGTCGTACCGATTAtc                                              >  2:4266188/1‑150 (MQ=255)
                                                                                                     aaaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtctc                                            >  2:7939/1‑150 (MQ=255)
                                                                                                     aaaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtctc                                            >  2:1447777/1‑150 (MQ=255)
                                                                                                      aaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTTTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCa                                           >  1:430663/1‑150 (MQ=255)
                                                                                                        cagcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAAc                                         <  1:1680473/150‑1 (MQ=255)
                                                                                                          gcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTa                                       <  2:5013848/150‑1 (MQ=255)
                                                                                                             ggCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCt                                    >  1:2948937/1‑150 (MQ=255)
                                                                                                             ggCATCGGACGAAAAACATCGTCGCTTACACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTTGAATTATATTGGCGATCAGCAAAAGGCGCTTTATTACAACGCCTTCAGTCGCCTGTGTCGTACCGATTATCTCAACTATCt                                    >  2:1393904/1‑150 (MQ=255)
                                                                                                              gCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtg                                   <  1:325552/150‑1 (MQ=255)
                                                                                                                 tCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgc                                >  2:2954124/1‑150 (MQ=255)
                                                                                                                  cGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcg                               >  1:3668909/1‑150 (MQ=255)
                                                                                                                  cGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTACGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcg                               >  1:3941401/1‑150 (MQ=255)
                                                                                                                    gACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGGTTGGGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTACGTCGCCTGTGTCGTAACGATTATATAAACTATCttcgcgtg                             >  1:3600241/1‑150 (MQ=255)
                                                                                                                     aCGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgc                            <  2:3052435/150‑1 (MQ=255)
                                                                                                                       gAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcgc                          >  2:4845966/1‑150 (MQ=255)
                                                                                                                        aaaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcgcg                         >  2:5182621/1‑150 (MQ=255)
                                                                                                                          aaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGaa                       >  1:708803/1‑150 (MQ=255)
                                                                                                                            aCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATg                     >  2:3381920/1‑150 (MQ=255)
                                                                                                                             cATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCGGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATgg                    <  2:2611730/150‑1 (MQ=255)
                                                                                                                                 gtcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg                >  1:1059864/1‑150 (MQ=255)
                                                                                                                                  tcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGa               >  1:3049386/1‑150 (MQ=255)
                                                                                                                                  tcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGa               >  1:3860362/1‑150 (MQ=255)
                                                                                                                                   cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGat              >  1:4124706/1‑150 (MQ=255)
                                                                                                                                    gCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGata             <  2:2117828/150‑1 (MQ=255)
                                                                                                                                    gCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGata             >  2:1520499/1‑150 (MQ=255)
                                                                                                                                     cTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGatat            >  1:5044220/1‑150 (MQ=255)
                                                                                                                                      ttCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGACTTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATc           >  1:3130631/1‑150 (MQ=255)
                                                                                                                                      ttCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGCTATc           >  2:639518/1‑150 (MQ=255)
                                                                                                                                         cACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTac        <  1:2789182/150‑1 (MQ=255)
                                                                                                                                          aCGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTaca       >  2:1051269/1‑150 (MQ=255)
                                                                                                                                           cGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTacac      <  1:81498/150‑1 (MQ=255)
                                                                                                                                            gACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACAcc     >  1:4833065/1‑150 (MQ=255)
                                                                                                                                             aCAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACAccc    <  1:4641471/150‑1 (MQ=255)
                                                                                                                                               aaaaAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCCAg  >  1:2874781/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCCAG  >  CP009273/1478825‑1479117

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑