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,9710GA94.5% 1024.6 / 35.2 331E552K (GAG→AAG) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base A (154/157);  total (164/167)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 336 total aligned reads displayed.

GCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCCAG  >  CP009273/1478824‑1479117
                                                                                                                                                   |                                                                                                                                                  
gcgcGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAaaag                                                                                                                                                  <  1:742907/150‑1 (MQ=255)
gcgcGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAaaag                                                                                                                                                  >  1:387433/1‑150 (MQ=255)
 cgcgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGt                                                                                                                                                 <  2:596550/150‑1 (MQ=255)
  gcgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTc                                                                                                                                                >  2:2086351/1‑150 (MQ=255)
    gTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTg                                                                                                                                              <  2:4564305/150‑1 (MQ=255)
     tATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGa                                                                                                                                             <  2:4302982/150‑1 (MQ=255)
     tATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGa                                                                                                                                             <  2:3833175/150‑1 (MQ=255)
        ggTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACtt                                                                                                                                          <  1:5062434/150‑1 (MQ=255)
         gTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACttt                                                                                                                                         <  2:4504594/150‑1 (MQ=255)
            cTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtctc                                                                                                                                      >  1:2057596/1‑150 (MQ=255)
             tGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcg                                                                                                                                     >  1:3747412/1‑150 (MQ=255)
               gAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                   >  2:637822/1‑150 (MQ=255)
               gAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                   <  2:4031325/150‑1 (MQ=255)
                aaGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGt                                                                                                                                  <  2:3467854/150‑1 (MQ=255)
                 aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGtt                                                                                                                                 <  2:4684137/150‑1 (MQ=255)
                  gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGttt                                                                                                                                <  2:1290525/150‑1 (MQ=255)
                  gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGttt                                                                                                                                >  1:4335343/1‑150 (MQ=255)
                   cgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtg                                                                                                                               >  1:2690521/1‑150 (MQ=255)
                      aaaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGa                                                                                                                            >  2:3951109/1‑150 (MQ=255)
                         aaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATc                                                                                                                         <  1:2511676/150‑1 (MQ=255)
                         aaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATc                                                                                                                         >  2:2644814/1‑150 (MQ=255)
                         aaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATc                                                                                                                         <  1:1763226/150‑1 (MQ=255)
                              ggCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTgg                                                                                                                    >  1:1678350/1‑150 (MQ=255)
                              ggCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTgg                                                                                                                    >  1:1968453/1‑150 (MQ=255)
                                cTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGaa                                                                                                                  >  1:2341824/1‑150 (MQ=255)
                                  gcgtgcgtGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAAtt                                                                                                                <  1:799116/150‑1 (MQ=255)
                                   cgtgcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTa                                                                                                               >  2:4999042/1‑150 (MQ=255)
                                       cgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATTGACAAACAGCAGGCCTCGGACGTAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTGTCGCGGGTGTGAATCTGTGGAATTagagt                                                                                                           >  2:4828946/1‑146 (MQ=255)
                                       cgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCtt                                                                                                           >  2:4801365/1‑150 (MQ=255)
                                         tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTgg                                                                                                         <  1:2711866/150‑1 (MQ=255)
                                          gAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGc                                                                                                        >  1:2483191/1‑150 (MQ=255)
                                           aGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCg                                                                                                       <  1:3162449/150‑1 (MQ=255)
                                            ggCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACtt                                                                                                                                          <  1:4827389/114‑1 (MQ=255)
                                            ggCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGa                                                                                                      >  1:2286417/1‑150 (MQ=255)
                                             gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGag                                                                                                     >  1:4860495/1‑150 (MQ=255)
                                             gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGag                                                                                                     <  2:1530568/150‑1 (MQ=255)
                                              cGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagc                                                                                                    <  2:4179115/150‑1 (MQ=255)
                                              cGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagc                                                                                                    >  2:187322/1‑150 (MQ=255)
                                                atgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagcag                                                                                                  >  2:3499524/1‑150 (MQ=255)
                                                atgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagcag                                                                                                  <  2:355917/150‑1 (MQ=255)
                                                       tCACGTCCGCGCTCTCCAGTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAgg                                                                                           <  2:2997264/150‑1 (MQ=255)
                                                        cACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgc                                                                                          >  1:2971897/1‑150 (MQ=255)
                                                          cGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgcgc                                                                                        <  2:3129769/150‑1 (MQ=255)
                                                           gTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCt                                                                                       >  1:1784945/1‑150 (MQ=255)
                                                           gTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCt                                                                                       >  1:146425/1‑150 (MQ=255)
                                                               gcgcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttct                                                                                   <  2:2889697/150‑1 (MQ=255)
                                                               gcgcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttct                                                                                   >  1:815865/1‑150 (MQ=255)
                                                                cgcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttctt                                                                                  <  2:2755151/150‑1 (MQ=255)
                                                                  ctctCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTcc                                                                                >  2:4458567/1‑150 (MQ=255)
                                                                  ctctCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTcc                                                                                <  1:101936/150‑1 (MQ=255)
                                                                       cATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGc                                                                           >  1:4747918/1‑150 (MQ=255)
                                                                        aTTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGcc                                                                          >  1:4197124/1‑150 (MQ=255)
                                                                          tCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCtt                                                                        >  2:2508830/1‑150 (MQ=255)
                                                                             ggATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCg                                                                     >  1:4846278/1‑150 (MQ=255)
                                                                              gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGt                                                                    >  2:2328100/1‑150 (MQ=255)
                                                                                tCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCg                                                                  >  1:1299344/1‑150 (MQ=255)
                                                                                   gcgcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCt                                                               <  2:1699240/150‑1 (MQ=255)
                                                                                     gcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCtgt                                                             <  1:66473/150‑1 (MQ=255)
                                                                                        aaCGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTc                                                          >  1:3768174/1‑150 (MQ=255)
                                                                                         aCGTCCGATGGACAAACAGCAGGCATCGGACGAAACACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCg                                                         <  2:3687427/150‑1 (MQ=255)
                                                                                           gTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTa                                                       >  1:4101551/1‑150 (MQ=255)
                                                                                            tCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTAc                                                      <  2:4637018/150‑1 (MQ=255)
                                                                                             ccGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAACAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTAcc                                                     <  2:3366020/150‑1 (MQ=255)
                                                                                                aTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGAt                                                  >  1:2125025/1‑150 (MQ=255)
                                                                                                 tGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGAtt                                                 <  2:1619737/150‑1 (MQ=255)
                                                                                                     cAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                             <  1:3693926/150‑1 (MQ=255)
                                                                                                      aaaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtctc                                            <  1:532845/150‑1 (MQ=255)
                                                                                                         cagcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAAc                                         <  1:2794000/150‑1 (MQ=255)
                                                                                                          agcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                        <  1:4999042/150‑1 (MQ=255)
                                                                                                           gcagTCATCGTACTAAAAACATCGTCGCTTCCACGAGAAATAGTCTGACTTTCTCGCGTTTGTGAATCTGTTGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTATTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTa                                       <  1:2886491/150‑1 (MQ=255)
                                                                                                             agGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGGTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATc                                     >  2:1808154/1‑150 (MQ=255)
                                                                                                             agGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCTACTATc                                     >  2:4747358/1‑150 (MQ=255)
                                                                                                              ggCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCt                                    <  2:1148779/150‑1 (MQ=255)
                                                                                                                cATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgc                                  >  1:4897477/1‑150 (MQ=255)
                                                                                                                 aTCGGACGAACAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcg                                 <  2:1005062/150‑1 (MQ=255)
                                                                                                                     gACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtg                             >  2:2219491/1‑150 (MQ=255)
                                                                                                                     gACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtg                             <  1:4549060/150‑1 (MQ=255)
                                                                                                                       cGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcg                           >  2:5135554/1‑150 (MQ=255)
                                                                                                                       cGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcg                           >  1:1190415/1‑150 (MQ=255)
                                                                                                                       cGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcg                           <  2:863694/150‑1 (MQ=255)
                                                                                                                        gAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcgc                          >  1:4627300/1‑150 (MQ=255)
                                                                                                                         aaaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcgcg                         >  2:1299777/1‑150 (MQ=255)
                                                                                                                          aaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTCCCGATTATCTCAACTATCTGCGCGTGCGCGa                        >  2:3764531/1‑150 (MQ=255)
                                                                                                                           aaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGaa                       >  1:4314229/1‑150 (MQ=255)
                                                                                                                            aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAAt                      <  1:3599680/150‑1 (MQ=255)
                                                                                                                             aCATCGTCGCTTACACGACAAAAATTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATATTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCAGTGTCGTACCGATTATATAAACTATATGCGCGTGCGCGAATg                     >  1:3040397/1‑150 (MQ=255)
                                                                                                                              cATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATgg                    <  1:3951217/150‑1 (MQ=255)
                                                                                                                                tcgtcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCa                  <  1:3297788/150‑1 (MQ=255)
                                                                                                                                 cgtcgCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAg                 <  2:2254360/150‑1 (MQ=255)
                                                                                                                                 cgtcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAg                 >  1:1916401/1‑150 (MQ=255)
                                                                                                                                   tcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGa               >  2:941992/1‑150 (MQ=255)
                                                                                                                                    cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGat              >  1:962420/1‑150 (MQ=255)
                                                                                                                                    cgCATCCACGACAAAAAGTCTGACTTTCTCGCGTATGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTACAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGat              >  1:1057873/1‑150 (MQ=255)
                                                                                                                                     gCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGata             <  1:4517715/150‑1 (MQ=255)
                                                                                                                                      cTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGatat            >  1:1551520/1‑150 (MQ=255)
                                                                                                                                       ttCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATc           >  1:4400156/1‑150 (MQ=255)
                                                                                                                                         ccACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTa         <  2:4418991/150‑1 (MQ=255)
                                                                                                                                         ccACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTa         >  1:887942/1‑150 (MQ=255)
                                                                                                                                              aCAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACAccc    >  1:3404427/1‑150 (MQ=255)
                                                                                                                                               cAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCCa   <  1:298943/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCCAG  >  CP009273/1478824‑1479117

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: ‑