Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,573,623 C→A 7.5% pseudogene (3010/3567 nt) hsdR ← pseudogene, premature stop codon, host restriction; endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,573,6230CA 7.5% 1211.9 / 44.8 314pseudogene (3010/3567 nt)hsdRpseudogene, premature stop codon, host restriction; endonuclease R
Reads supporting (aligned to +/- strand):  ref base C (131/159);  new base A (9/15);  total (140/174)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

GAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGC  >  CP009273/4573474‑4573772
                                                                                                                                                     |                                                                                                                                                     
gagaTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATc                                                                                                                                                       <  2:5031113/150‑1 (MQ=255)
  gaTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGACCGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGa                                                                                                                                                     >  2:1971285/1‑150 (MQ=255)
   aTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAg                                                                                                                                                    >  2:3955329/1‑150 (MQ=255)
     cgcgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGa                                                                                                                                                  <  1:1932113/150‑1 (MQ=255)
      gcgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGaa                                                                                                                                                 <  1:1728202/150‑1 (MQ=255)
       cgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAg                                                                                                                                                >  1:2134713/1‑150 (MQ=255)
        gcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGa                                                                                                                                               <  2:1745620/150‑1 (MQ=255)
         cgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGAt                                                                                                                                              <  2:344425/150‑1 (MQ=255)
         cgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGAt                                                                                                                                              <  2:3132879/150‑1 (MQ=255)
          gggCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATc                                                                                                                                             <  1:3187219/150‑1 (MQ=255)
          gggCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATc                                                                                                                                             <  2:981501/150‑1 (MQ=255)
                ttAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGca                                                                                                                                       <  2:1047837/150‑1 (MQ=255)
                 taataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcat                                                                                                                                      >  2:906589/1‑150 (MQ=255)
                 taataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcat                                                                                                                                      >  2:3481519/1‑150 (MQ=255)
                 taataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcat                                                                                                                                      >  2:1310764/1‑150 (MQ=255)
                 taataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcat                                                                                                                                      <  1:3612037/150‑1 (MQ=255)
                  aataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcatc                                                                                                                                     <  1:4583418/150‑1 (MQ=255)
                  aataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcatc                                                                                                                                     >  1:1564379/1‑150 (MQ=255)
                   ataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcatca                                                                                                                                    >  2:436421/1‑150 (MQ=255)
                   ataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcatca                                                                                                                                    >  2:1249242/1‑150 (MQ=255)
                   ataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcatca                                                                                                                                    <  1:2771899/150‑1 (MQ=255)
                     aaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCAtt                                                                                                                                  <  2:4144021/150‑1 (MQ=255)
                       ctgcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCa                                                                                                                                <  1:3367715/150‑1 (MQ=255)
                        tgcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAg                                                                                                                               >  1:995722/1‑150 (MQ=255)
                          cctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgt                                                                                                                             >  2:880805/1‑150 (MQ=255)
                          cctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgt                                                                                                                             <  1:1526308/150‑1 (MQ=255)
                            tgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttg                                                                                                                           <  1:1031443/150‑1 (MQ=255)
                            tgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttg                                                                                                                           >  2:3186296/1‑150 (MQ=255)
                            tgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttg                                                                                                                           >  1:117553/1‑150 (MQ=255)
                             gcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttgt                                                                                                                          <  1:132238/150‑1 (MQ=255)
                             gcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttgt                                                                                                                          <  2:1718670/150‑1 (MQ=255)
                              cAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttgtt                                                                                                                         <  2:2471378/150‑1 (MQ=255)
                                aTGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGa                                                                                                                       >  1:2122463/1‑150 (MQ=255)
                                       ttGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTGCACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGt                                                                                                                >  1:1783162/1‑150 (MQ=255)
                                        tGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGtt                                                                                                               >  2:4183819/1‑150 (MQ=255)
                                          cgcgTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGtttt                                                                                                             <  2:2897415/150‑1 (MQ=255)
                                             gTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGACGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAg                                                                                                          >  1:2038794/1‑150 (MQ=255)
                                             gTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAg                                                                                                          >  1:3700556/1‑150 (MQ=255)
                                             gTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAg                                                                                                          <  1:2903425/150‑1 (MQ=255)
                                             gTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAg                                                                                                          <  2:1834392/150‑1 (MQ=255)
                                               tCGGGGAACCTTGCACCAGCGAGTCAAACGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATATTACAAACACAACTTCATCGTCGATATCGAGGAAGATCCGCGCATCATTCAGGATGTTGATTTCCGTTTTCAGca                                                                                                        >  2:492607/1‑149 (MQ=255)
                                                  gggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCtttt                                                                                                     >  2:4991083/1‑150 (MQ=255)
                                                   ggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTc                                                                                                    >  2:1801527/1‑150 (MQ=255)
                                                    gAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTcc                                                                                                   <  2:2708226/150‑1 (MQ=255)
                                                    gAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTcc                                                                                                   <  2:3911322/150‑1 (MQ=255)
                                                      aCGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAg                                                                                                 >  2:2546402/1‑150 (MQ=255)
                                                       cGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGa                                                                                                <  2:4162446/150‑1 (MQ=255)
                                                         ttGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACg                                                                                              >  2:214375/1‑150 (MQ=255)
                                                          tGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACgg                                                                                             >  2:1380204/1‑150 (MQ=255)
                                                            cACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGc                                                                                           <  1:1053532/150‑1 (MQ=255)
                                                             aCCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCg                                                                                          >  1:1426459/1‑150 (MQ=255)
                                                               cAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGAt                                                                                        <  2:580199/150‑1 (MQ=255)
                                                                aGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATa                                                                                       <  2:4017664/150‑1 (MQ=255)
                                                                aGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATa                                                                                       >  1:3760685/1‑150 (MQ=255)
                                                                 gCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATaa                                                                                      <  1:1371738/150‑1 (MQ=255)
                                                                    aGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGc                                                                                   >  1:4940287/1‑150 (MQ=255)
                                                                      tCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCa                                                                                 >  2:926312/1‑150 (MQ=255)
                                                                      tCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCa                                                                                 <  2:151191/150‑1 (MQ=255)
                                                                        aaaGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAgg                                                                               <  1:1960208/150‑1 (MQ=255)
                                                                          aGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGGTAAAGCCAGGCa                                                                             <  1:4250603/150‑1 (MQ=255)
                                                                              ttCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGttt                                                                         <  1:892639/150‑1 (MQ=255)
                                                                              ttCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGttt                                                                         <  1:4099254/150‑1 (MQ=255)
                                                                              ttCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGttt                                                                         <  1:874465/150‑1 (MQ=255)
                                                                                 gAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGtt                                                                      >  1:677173/1‑150 (MQ=255)
                                                                                    gAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTaaa                                                                   <  1:348047/150‑1 (MQ=255)
                                                                                    gAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTaaa                                                                   >  1:4502529/1‑150 (MQ=255)
                                                                                       aTCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGAc                                                                <  2:2295781/150‑1 (MQ=255)
                                                                                        tCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACt                                                               <  2:1928594/150‑1 (MQ=255)
                                                                                           tGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCg                                                            >  1:2280407/1‑150 (MQ=255)
                                                                                             cggcgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGgc                                                          >  2:115441/1‑150 (MQ=255)
                                                                                             cggcgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGgc                                                          <  1:4782873/150‑1 (MQ=255)
                                                                                                cgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCt                                                       >  2:2427883/1‑150 (MQ=255)
                                                                                                   gtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCa                                                    >  2:17900/1‑150 (MQ=255)
                                                                                                    tCGTACTCACCGTACAGCGATTTTACACTCACCACGTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAg                                                   <  2:4018613/150‑1 (MQ=255)
                                                                                                    tCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAg                                                   >  2:3306872/1‑150 (MQ=255)
                                                                                                       tAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGc                                                >  2:3743886/1‑150 (MQ=255)
                                                                                                       tAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGc                                                >  1:2814472/1‑150 (MQ=255)
                                                                                                        aaTCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCg                                               >  1:569875/1‑150 (MQ=255)
                                                                                                         aTCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCgg                                              >  2:3555587/1‑150 (MQ=255)
                                                                                                            aCCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGccc                                           <  1:4296143/150‑1 (MQ=255)
                                                                                                            aCCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGccc                                           <  1:2502453/150‑1 (MQ=255)
                                                                                                            aCCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGccc                                           <  2:428839/150‑1 (MQ=255)
                                                                                                               gTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCttt                                        >  1:4616052/1‑150 (MQ=255)
                                                                                                                  cAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTc                                     <  1:3536308/150‑1 (MQ=255)
                                                                                                                  cAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTc                                     <  1:803682/150‑1 (MQ=255)
                                                                                                                   aGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTcc                                    <  1:3955329/150‑1 (MQ=255)
                                                                                                                     cGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTccc                                   >  2:4575844/1‑149 (MQ=255)
                                                                                                                             cacTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGcgcg                          <  1:1131403/150‑1 (MQ=255)
                                                                                                                              acTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGa                         <  1:5016386/150‑1 (MQ=255)
                                                                                                                                tCACCACTTCCTCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAgg                       <  2:5000701/150‑1 (MQ=255)
                                                                                                                                 caccacTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGc                      <  2:4777336/150‑1 (MQ=255)
                                                                                                                                    cacTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGaa                   <  1:5135427/150‑1 (MQ=255)
                                                                                                                                    cacTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGaa                   >  1:4842129/1‑150 (MQ=255)
                                                                                                                                     acTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAg                  >  1:2676369/1‑150 (MQ=255)
                                                                                                                                     acTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAg                  >  1:1664887/1‑150 (MQ=255)
                                                                                                                                       ttCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGcc                >  2:1420557/1‑150 (MQ=255)
                                                                                                                                           tcgtcgATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTa            >  1:592653/1‑150 (MQ=255)
                                                                                                                                            cgtcgATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTaa           >  2:3861936/1‑150 (MQ=255)
                                                                                                                                                gATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGtt       >  2:2114694/1‑150 (MQ=255)
                                                                                                                                                  tatCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCa     <  2:628727/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGC  >  CP009273/4573474‑4573772

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: