Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 150,452 C→A 24.6% E488* (GAA→TAA)  htrE ← putative outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273150,4520CA24.6% 868.3 / 214.2 376E488* (GAA→TAA) htrEputative outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base C (142/141);  new base A (52/41);  total (194/182)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 382 total aligned reads displayed.

TCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATG  >  CP009273/150303‑150601
                                                                                                                                                     |                                                                                                                                                     
tCCGGAAGCACAGTACTCGGACCAACTTCCTGATAGATACAATGATCCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGGACCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTc                                                                                                                                                       <  1:1090166/150‑1 (MQ=255)
 ccGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCg                                                                                                                                                      <  2:2390091/150‑1 (MQ=255)
 ccGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCg                                                                                                                                                      <  2:3832600/150‑1 (MQ=255)
  cGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                                                                                                                                     <  2:1012413/150‑1 (MQ=255)
     aaGCACAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                                                                                                                                  <  1:3727162/150‑1 (MQ=255)
       gCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGt                                                                                                                                                >  1:2486637/1‑150 (MQ=255)
            gTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCAc                                                                                                                                           <  1:4003155/150‑1 (MQ=255)
               aTCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtc                                                                                                                                        >  1:1735552/1‑150 (MQ=255)
                 cGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcat                                                                                                                                      <  1:1951534/150‑1 (MQ=255)
                  ggACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc                                                                                                                                     >  2:4272327/1‑150 (MQ=255)
                   gACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatca                                                                                                                                    <  1:4268647/150‑1 (MQ=255)
                     ccAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATATCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAAt                                                                                                                                  <  2:4197954/150‑1 (MQ=255)
                       aaCTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa                                                                                                                                >  1:1425036/1‑150 (MQ=255)
                       aaCTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa                                                                                                                                >  2:5241544/1‑150 (MQ=255)
                            ccGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt                                                                                                                           <  1:4771009/150‑1 (MQ=255)
                              ggAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAg                                                                                                                         <  2:2090288/150‑1 (MQ=255)
                                 aaGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCTAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGc                                                                                                                      >  2:3050274/1‑150 (MQ=255)
                                 aaGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGc                                                                                                                      >  2:1909541/1‑150 (MQ=255)
                                 aaGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGc                                                                                                                      >  2:512624/1‑150 (MQ=255)
                                 aaGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGc                                                                                                                      >  2:2640427/1‑150 (MQ=255)
                                 aaGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGc                                                                                                                      >  2:4609882/1‑150 (MQ=255)
                                        aaTGAACCGTAATCTTTTTTCTCAAATTACAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCAGCATTa                                                                                                               >  1:1001048/1‑150 (MQ=255)
                                         aTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTaa                                                                                                              <  1:1974572/150‑1 (MQ=255)
                                         aTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTaa                                                                                                              >  2:247366/1‑150 (MQ=255)
                                         aTGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTaa                                                                                                              >  1:3026720/1‑150 (MQ=255)
                                          tgcaCCGTAATCTTTTTTCTCAAATTTCATCGGTTGGTTAATACTGACCGGAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAg                                                                                                             <  2:2303229/147‑1 (MQ=255)
                                             aCCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGAcc                                                                                                          >  2:4229722/1‑150 (MQ=255)
                                               cGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCaa                                                                                                        <  2:2640158/150‑1 (MQ=255)
                                               cGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTGAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCaa                                                                                                        >  2:2730194/1‑150 (MQ=255)
                                                gTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTAACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAg                                                                                                       >  1:505465/1‑150 (MQ=255)
                                                gTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAg                                                                                                       >  2:4459129/1‑150 (MQ=255)
                                                 tAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAgg                                                                                                      <  1:5281736/150‑1 (MQ=255)
                                                 tAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAgg                                                                                                      >  1:69161/1‑150 (MQ=255)
                                                   aTCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTa                                                                                                    >  1:2645495/1‑150 (MQ=255)
                                                    tCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTaa                                                                                                   >  1:3204729/1‑150 (MQ=255)
                                                     cTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAAt                                                                                                  <  1:214556/150‑1 (MQ=255)
                                                      ttttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAAtt                                                                                                 <  2:3160727/150‑1 (MQ=255)
                                                      ttttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAAtt                                                                                                 <  2:151195/150‑1 (MQ=255)
                                                        ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCt                                                                                               >  1:4701324/1‑150 (MQ=255)
                                                         tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATATCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCtg                                                                                              >  2:2827936/1‑150 (MQ=255)
                                                             tcAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTc                                                                                          <  2:3596068/150‑1 (MQ=255)
                                                                 aTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCGAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAAt                                                                                      >  1:2079063/1‑150 (MQ=255)
                                                                 aTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAAt                                                                                      >  2:3507129/1‑150 (MQ=255)
                                                                  tttCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATa                                                                                     <  2:2926378/150‑1 (MQ=255)
                                                                     cAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCg                                                                                  <  1:890385/150‑1 (MQ=255)
                                                                         ggttggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAg                                                                              >  1:4430470/1‑150 (MQ=255)
                                                                          gttggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAgg                                                                             <  2:3569707/150‑1 (MQ=255)
                                                                           ttggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGc                                                                            <  2:4896141/150‑1 (MQ=255)
                                                                            tggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGcc                                                                           >  1:1498123/1‑150 (MQ=255)
                                                                            tggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGAGCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGcc                                                                           >  1:4343425/1‑150 (MQ=255)
                                                                             ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCcg                                                                          <  2:3193/150‑1 (MQ=255)
                                                                             ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCcg                                                                          <  1:2280746/150‑1 (MQ=255)
                                                                                 aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGat                                                                      <  2:2358666/150‑1 (MQ=255)
                                                                                 aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGat                                                                      >  1:4099804/1‑150 (MQ=255)
                                                                                  aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGata                                                                     <  2:3966072/150‑1 (MQ=255)
                                                                                    aCTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATAtt                                                                   <  2:831204/150‑1 (MQ=255)
                                                                                      tGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCa                                                                 <  2:4136445/150‑1 (MQ=255)
                                                                                      tGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCa                                                                 >  1:1905724/1‑150 (MQ=255)
                                                                                       gACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTTATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAg                                                                >  1:4554027/1‑150 (MQ=255)
                                                                                       gACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAg                                                                <  2:3715563/150‑1 (MQ=255)
                                                                                             aaCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGt                                                          >  2:1193011/1‑150 (MQ=255)
                                                                                                 tGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACtt                                                      <  2:5083058/150‑1 (MQ=255)
                                                                                                  gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTg                                                     >  1:29448/1‑150 (MQ=255)
                                                                                                   aTTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGt                                                    >  2:1637886/1‑150 (MQ=255)
                                                                                                   aTTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGt                                                    >  2:55253/1‑150 (MQ=255)
                                                                                                    tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGtt                                                   <  1:4902200/150‑1 (MQ=255)
                                                                                                      tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTttc                                                 >  1:3490492/1‑150 (MQ=255)
                                                                                                        tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTttctt                                               <  2:3255040/150‑1 (MQ=255)
                                                                                                         cATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTttcttc                                              >  1:5102755/1‑150 (MQ=255)
                                                                                                             cGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATGAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAAt                                          <  2:2199569/150‑1 (MQ=255)
                                                                                                              gTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATa                                         >  2:4989464/1‑150 (MQ=255)
                                                                                                                  gTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACtt                                     <  2:1791441/150‑1 (MQ=255)
                                                                                                                   tAATTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTg                                    <  1:1274185/150‑1 (MQ=255)
                                                                                                                    aaTTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGAGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAAATTGt                                   >  2:3417941/1‑150 (MQ=255)
                                                                                                                    aaTTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGt                                   <  1:3288765/150‑1 (MQ=255)
                                                                                                                     aTTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGtt                                  >  1:3248705/1‑150 (MQ=255)
                                                                                                                     aTTACGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGtt                                  >  2:672114/1‑150 (MQ=255)
                                                                                                                       tACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTcc                                >  1:3608992/1‑150 (MQ=255)
                                                                                                                        aCGCATGGATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCTATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCa                               <  2:1139469/150‑1 (MQ=255)
                                                                                                                           cATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGa                            >  1:585844/1‑150 (MQ=255)
                                                                                                                              gtaTTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTTTTCCAGGAAac                         >  2:2071507/3‑150 (MQ=255)
                                                                                                                              ggATTTCGGTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAac                         >  1:1710099/1‑150 (MQ=255)
                                                                                                                                aTTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAacac                       >  2:4975266/1‑150 (MQ=255)
                                                                                                                                aTTTCGGTTCAAGATCTTGTTCGGGATGTATCACTTCATCAATTAGAGTTAGAGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGAGAAGTACTTGATACTTCGAATAACTTGTTCCAGGAAacac                       >  1:3544316/1‑150 (MQ=255)
                                                                                                                                     ggTTCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATaa                  >  2:2503681/1‑150 (MQ=255)
                                                                                                                                        tCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACtt               <  2:4807053/150‑1 (MQ=255)
                                                                                                                                        tCAAGATCTTGTTAGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACtt               <  1:3848868/150‑1 (MQ=255)
                                                                                                                                         cAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCTCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACttt              <  2:4675814/150‑1 (MQ=255)
                                                                                                                                         cAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACttt              <  1:1064628/150‑1 (MQ=255)
                                                                                                                                         cAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACttt              <  2:1425036/150‑1 (MQ=255)
                                                                                                                                          aaGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTg             >  1:3401918/1‑150 (MQ=255)
                                                                                                                                          aaGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTg             <  2:3300325/150‑1 (MQ=255)
                                                                                                                                             aTCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGccc          >  1:485968/1‑150 (MQ=255)
                                                                                                                                               cTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCt        <  1:4246135/150‑1 (MQ=255)
                                                                                                                                               cTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCt        <  2:1075483/150‑1 (MQ=255)
                                                                                                                                               cTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCt        <  1:2817339/150‑1 (MQ=255)
                                                                                                                                                  gTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGt     >  1:1674523/1‑150 (MQ=255)
                                                                                                                                                  gTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGt     >  1:1148428/1‑150 (MQ=255)
                                                                                                                                                  gTTCGGGATGTTTCACTTCATCAATTAGAGTAAGAGCATCATTAAGACAAAGGTAATTCTGTGTCGAATAGCGATAGGCCGAGATATTCAGTGAAGTACTTGTTTCTTCGAAAAACTTGTTCCAGGAAACACGATAAATTTGCCCCTGGt     >  2:2582274/1‑150 (MQ=255)
                                                                                                                                                     cGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATg  <  1:2730318/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCGGAAGCCCAGTAATCGGACCAACTTCCGGAAAGATAAAATGAACCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCATTAAGACCAAGGTAATTCTGTGTCGAATAGCGATAGGCCGCGATATTCAGTGAAGTACTTGTTTCTTCGAATAACTTGTTCCAGGAAACACGATAACTTTGCCCCTGGTATG  >  CP009273/150303‑150601

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: