Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,573,623 C→A 5.9% 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 5.9% 1079.7 / 29.1 285pseudogene (3010/3567 nt)hsdRpseudogene, premature stop codon, host restriction; endonuclease R
Reads supporting (aligned to +/- strand):  ref base C (140/128);  new base A (9/8);  total (149/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 286 total aligned reads displayed.

ATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGC  >  CP009273/4573477‑4573772
                                                                                                                                                  |                                                                                                                                                     
aTCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAg                                                                                                                                                    <  2:1365856/150‑1 (MQ=255)
 tCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAgg                                                                                                                                                   <  1:1073931/150‑1 (MQ=255)
  cgcgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGa                                                                                                                                                  >  1:1065895/1‑150 (MQ=255)
  cgcgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGa                                                                                                                                                  >  1:3036880/1‑150 (MQ=255)
   gcgcgGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGaa                                                                                                                                                 <  1:4413777/150‑1 (MQ=255)
       gggCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATc                                                                                                                                             <  2:566455/150‑1 (MQ=255)
       gggCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATc                                                                                                                                             <  2:1322413/150‑1 (MQ=255)
        ggCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCg                                                                                                                                            <  2:4513522/150‑1 (MQ=255)
           gATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGgcg                                                                                                                                         >  1:1379356/1‑150 (MQ=255)
           gATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATATTAAACTCACCACTTCATCGTCGATATCGAGGAAGATCGgcg                                                                                                                                         >  2:3901797/1‑150 (MQ=255)
            aTTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGACGATCGgcgc                                                                                                                                        >  2:2554746/1‑150 (MQ=255)
             ttAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGca                                                                                                                                       >  2:755569/1‑150 (MQ=255)
              taataaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcat                                                                                                                                      <  1:2758195/150‑1 (MQ=255)
               aagaaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGCGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcctc                                                                                                                                     <  2:3065506/147‑1 (MQ=255)
                  aaCTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCAtt                                                                                                                                  >  1:3253358/1‑150 (MQ=255)
                      gcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAgg                                                                                                                              <  1:3264097/150‑1 (MQ=255)
                      gcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAgg                                                                                                                              >  1:355065/1‑150 (MQ=255)
                      gcctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGTGGAAGATCGGCGCATCATTCAgg                                                                                                                              <  1:4457117/150‑1 (MQ=255)
                       cctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgt                                                                                                                             >  2:3686574/1‑150 (MQ=255)
                        ctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgtt                                                                                                                            <  2:2009967/150‑1 (MQ=255)
                        ctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgtt                                                                                                                            <  1:1170676/150‑1 (MQ=255)
                        ctgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgtt                                                                                                                            >  2:3444004/1‑150 (MQ=255)
                         tgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttg                                                                                                                           >  1:759253/1‑150 (MQ=255)
                         tgcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttg                                                                                                                           >  1:601882/1‑150 (MQ=255)
                            aaTGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGttgttg                                                                                                                        >  1:1024147/1‑150 (MQ=255)
                             aTGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGa                                                                                                                       >  1:333799/1‑150 (MQ=255)
                                ccGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGt                                                                                                                    >  1:878194/1‑150 (MQ=255)
                                 cGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTc                                                                                                                   <  1:549658/150‑1 (MQ=255)
                                  ggTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTcc                                                                                                                  >  1:2949387/1‑150 (MQ=255)
                                  ggTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTcc                                                                                                                  >  1:667088/1‑150 (MQ=255)
                                           ttCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGc                                                                                                         <  2:2364545/150‑1 (MQ=255)
                                            tCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCt                                                                                                        <  2:3790418/150‑1 (MQ=255)
                                              ggggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCttt                                                                                                      <  2:2782670/150‑1 (MQ=255)
                                              ggggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCttt                                                                                                      >  1:3612925/1‑150 (MQ=255)
                                              ggggAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGAGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCttt                                                                                                      <  2:272784/150‑1 (MQ=255)
                                                 gAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTcc                                                                                                   <  1:2207225/150‑1 (MQ=255)
                                                        gCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACggg                                                                                            <  2:50246/150‑1 (MQ=255)
                                                           ccAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGa                                                                                         >  1:3114182/1‑150 (MQ=255)
                                                             aGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATa                                                                                       >  2:305127/1‑150 (MQ=255)
                                                             aGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATa                                                                                       <  1:1148086/150‑1 (MQ=255)
                                                              gCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATag                                                                                      <  2:2809616/150‑2 (MQ=255)
                                                               cGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATaaa                                                                                     <  2:1049897/150‑1 (MQ=255)
                                                                 aGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGc                                                                                   <  2:3448009/150‑1 (MQ=255)
                                                                    cAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAg                                                                                >  2:1682094/1‑150 (MQ=255)
                                                                    cAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAg                                                                                >  2:2404870/1‑150 (MQ=255)
                                                                       aGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGTAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCa                                                                             >  1:2173215/1‑150 (MQ=255)
                                                                        ggCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAg                                                                            <  2:188945/150‑1 (MQ=255)
                                                                           ttCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGttt                                                                         >  2:3346417/1‑150 (MQ=255)
                                                                           ttCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGttt                                                                         >  2:1069423/1‑150 (MQ=255)
                                                                              gAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGtt                                                                      <  1:4103138/150‑1 (MQ=255)
                                                                               aGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTa                                                                     <  2:1874846/150‑1 (MQ=255)
                                                                                  aaaTCCTGCGGAGTGTCGTAATCACCGTACAGCGATGTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATACAGCCAGGCAGTTTGTTaatg                                                                  >  2:4407963/1‑148 (MQ=255)
                                                                                    aTCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGAc                                                                >  1:3707407/1‑150 (MQ=255)
                                                                                     tCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACt                                                               <  2:1453472/150‑1 (MQ=255)
                                                                                      ccTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACtt                                                              <  1:4314867/150‑1 (MQ=255)
                                                                                        tGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCg                                                            >  1:185840/1‑150 (MQ=255)
                                                                                        tGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCg                                                            <  1:3185393/150‑1 (MQ=255)
                                                                                        tGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCg                                                            <  1:3185449/150‑1 (MQ=255)
                                                                                             cgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCCTCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCt                                                       <  2:1161563/150‑1 (MQ=255)
                                                                                             cgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGACGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCt                                                       >  2:224889/1‑150 (MQ=255)
                                                                                             cgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGACGAAATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGATTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCt                                                       >  2:2724958/1‑150 (MQ=255)
                                                                                              gtgtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTc                                                      <  2:2957324/150‑1 (MQ=255)
                                                                                               tgtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTcc                                                     >  1:1534817/1‑150 (MQ=255)
                                                                                                gtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCa                                                    >  2:1765984/1‑150 (MQ=255)
                                                                                                  cGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGt                                                  <  1:960031/150‑1 (MQ=255)
                                                                                                   gTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTg                                                 <  1:554973/150‑1 (MQ=255)
                                                                                                    tAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGc                                                <  1:2358681/150‑1 (MQ=255)
                                                                                                    caaTCACCGTACAGCGACTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGc                                                <  1:2359274/149‑1 (MQ=255)
                                                                                                     aaTCACCGTACAGCGATTTTACACTCCCCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCAATCAGGGTGTTGATGTCCGTTTTCAGCTTTACCAGACGGGCGAAAAAGCCAGGCAGATTGTTAAAGACTTCGGCGCTCCAGTGCg                                               >  1:2570790/1‑150 (MQ=255)
                                                                                                           cGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCtt                                         <  1:4329658/150‑1 (MQ=255)
                                                                                                            gTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCttt                                        <  2:2319896/150‑1 (MQ=255)
                                                                                                              aCAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCttttt                                      >  1:4037755/1‑150 (MQ=255)
                                                                                                               cAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTc                                     >  2:698612/1‑150 (MQ=255)
                                                                                                               cAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTc                                     >  2:3301385/1‑150 (MQ=255)
                                                                                                                 gCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTccc                                   >  2:2360710/1‑150 (MQ=255)
                                                                                                                   gATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGc                                 >  1:997267/1‑150 (MQ=255)
                                                                                                                   gATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGc                                 >  1:1995875/1‑150 (MQ=255)
                                                                                                                   gATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGc                                 >  2:539453/1‑150 (MQ=255)
                                                                                                                     ttttACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAg                               >  1:3124358/1‑150 (MQ=255)
                                                                                                                      tttACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAgg                              >  2:2285788/1‑150 (MQ=255)
                                                                                                                       ttACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGgc                             >  1:1612893/1‑150 (MQ=255)
                                                                                                                       ttACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGgc                             >  1:2328538/1‑150 (MQ=255)
                                                                                                                        tACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGgcg                            >  2:2080900/1‑150 (MQ=255)
                                                                                                                        tACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGgcg                            >  1:2889328/1‑150 (MQ=255)
                                                                                                                             tCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAgg                       >  1:39225/1‑150 (MQ=255)
                                                                                                                                ccacTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGa                    <  1:3427609/150‑1 (MQ=255)
                                                                                                                                ccacTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGa                    <  1:283837/150‑1 (MQ=255)
                                                                                                                                  acTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAg                  >  1:1762017/1‑150 (MQ=255)
                                                                                                                                  acTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAg                  >  1:493773/1‑150 (MQ=255)
                                                                                                                                   cTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGc                 >  1:1909037/1‑150 (MQ=255)
                                                                                                                                     tCATCGTCGATATAGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCg               >  1:2067258/1‑150 (MQ=255)
                                                                                                                                     acctcgtcgATATCGAGGAAGATCGGCGCATCATTCAGGGTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCg               <  2:4230342/147‑1 (MQ=255)
                                                                                                                                          gtcgATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTaaa          <  2:249619/150‑1 (MQ=255)
                                                                                                                                          gtcgATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTaaa          >  2:1541818/1‑150 (MQ=255)
                                                                                                                                             gATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGtt       >  2:1390702/1‑150 (MQ=255)
                                                                                                                                              ttatCGAGGAAGATCGGCGCATCATTCAGGTTTTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTc      <  1:3901797/149‑1 (MQ=255)
                                                                                                                                                atCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCAc    <  1:2361140/150‑1 (MQ=255)
                                                                                                                                                atCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCAc    >  2:471883/1‑150 (MQ=255)
                                                                                                                                                atCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCAc    <  1:1788108/150‑1 (MQ=255)
                                                                                                                                                atCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCAc    <  2:665276/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGC  >  CP009273/4573477‑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: