Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,810,619 G→T 13.7% intergenic (+234/+24) katE → / ← chbG catalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,810,6190GT13.7% 763.7 / 67.5 227intergenic (+234/+24)katE/chbGcatalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base G (105/90);  new base T (11/20);  total (116/111)
Fisher's exact test for biased strand distribution p-value = 8.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 239 total aligned reads displayed.

GCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGATAGTATT  >  CP009273/1810470‑1810768
                                                                                                                                                     |                                                                                                                                                     
gCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcg                                                                                                                                                       <  2:250507/150‑1 (MQ=255)
  cTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgcc                                                                                                                                                     <  2:3618283/150‑1 (MQ=255)
  cTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgcc                                                                                                                                                     >  2:4221515/1‑150 (MQ=255)
   tACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgccg                                                                                                                                                    >  2:1676838/1‑150 (MQ=255)
   tACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgccg                                                                                                                                                    <  1:1005631/150‑1 (MQ=255)
      aaaaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCAt                                                                                                                                                 >  1:2889871/1‑150 (MQ=255)
       aaaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATc                                                                                                                                                >  1:4406353/1‑150 (MQ=255)
        aaCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCt                                                                                                                                               >  1:871522/1‑150 (MQ=255)
         aCTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTg                                                                                                                                              <  2:1343440/150‑1 (MQ=255)
           tgtgCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGc                                                                                                                                            >  1:4167072/1‑150 (MQ=255)
           tgtgCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGc                                                                                                                                            >  1:2602929/1‑150 (MQ=255)
                 aaTTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACACCGCATCTGGCAAACAt                                                                                                                                      <  1:2723204/150‑1 (MQ=255)
                   ttCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTTCGTTAGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCCGATGCAAGGTTCACGCCGCATCTGGCAAACATcc                                                                                                                                    >  2:3983040/1‑150 (MQ=255)
                    tCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCt                                                                                                                                   <  1:2212778/150‑1 (MQ=255)
                    tCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCt                                                                                                                                   <  1:1574680/150‑1 (MQ=255)
                     cAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTc                                                                                                                                  <  2:2297134/150‑1 (MQ=255)
                        tatatTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACt                                                                                                                               >  2:1903943/1‑150 (MQ=255)
                            ttGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTaca                                                                                                                           >  1:111326/1‑150 (MQ=255)
                               cAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACAt                                                                                                                        >  2:3216863/1‑150 (MQ=255)
                                 ggAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAAGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATcc                                                                                                                      <  2:830542/150‑1 (MQ=255)
                                  gAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATccc                                                                                                                     >  2:1972861/1‑150 (MQ=255)
                                   aaaCACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCg                                                                                                                    >  1:2661542/1‑150 (MQ=255)
                                      cacGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATa                                                                                                                 <  1:2048637/150‑1 (MQ=255)
                                        cGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAAc                                                                                                               <  2:802375/150‑1 (MQ=255)
                                             gCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTcccc                                                                                                          <  1:760521/150‑1 (MQ=255)
                                              ccTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCa                                                                                                         >  2:3131558/1‑150 (MQ=255)
                                               cTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCACACCTCCTCACTTACACATCCCGATAACTCCCCaa                                                                                                        >  1:3744500/1‑150 (MQ=255)
                                                tGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAAc                                                                                                       >  1:2311463/1‑150 (MQ=255)
                                                  aTAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCg                                                                                                     >  1:2295547/1‑150 (MQ=255)
                                                       cgaagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAAc                                                                                                <  1:3270020/150‑1 (MQ=255)
                                                        gaagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAAcc                                                                                               <  2:603568/150‑1 (MQ=255)
                                                          agcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATCACTCCCCAACCGATAACCAc                                                                                             >  2:266362/1‑150 (MQ=255)
                                                          agcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCAc                                                                                             >  2:3330981/1‑150 (MQ=255)
                                                               aTCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGa                                                                                        >  1:3579415/1‑150 (MQ=255)
                                                               aTCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGa                                                                                        <  2:4109523/150‑1 (MQ=255)
                                                                  aGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCg                                                                                     <  2:1273449/150‑1 (MQ=255)
                                                                  aGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCg                                                                                     <  1:2740795/150‑1 (MQ=255)
                                                                   ggCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGa                                                                                    >  2:4554593/1‑150 (MQ=255)
                                                                      aGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAg                                                                                 <  2:847354/150‑1 (MQ=255)
                                                                      aGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAg                                                                                 >  1:3942211/1‑150 (MQ=255)
                                                                         tttGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCAc                                                                              <  1:2260508/150‑1 (MQ=255)
                                                                            gCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCtt                                                                           >  2:1777484/1‑150 (MQ=255)
                                                                            gCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCtt                                                                           <  1:696879/150‑1 (MQ=255)
                                                                               tttGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCa                                                                        >  2:4483630/1‑150 (MQ=255)
                                                                               tttGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCa                                                                        <  2:3001221/150‑1 (MQ=255)
                                                                                ttGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCaa                                                                       <  2:3907141/150‑1 (MQ=255)
                                                                                     agcagTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCACCCGATAACCACGCTGAGCGATAGCACCTTTCAacgacg                                                                  >  2:1665014/1‑150 (MQ=255)
                                                                                     agcagTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAacgacg                                                                  <  2:2972860/150‑1 (MQ=255)
                                                                                     agcagTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAacgacg                                                                  >  2:4255747/1‑150 (MQ=255)
                                                                                          tctcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGAt                                                             >  1:954203/1‑150 (MQ=255)
                                                                                            tcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGt                                                           >  1:1629072/1‑150 (MQ=255)
                                                                                              aaGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCa                                                         >  1:1372373/1‑150 (MQ=255)
                                                                                              aaGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCa                                                         >  1:1372363/1‑150 (MQ=255)
                                                                                               aGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCaa                                                        >  1:4559605/1‑150 (MQ=255)
                                                                                                 cggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAaca                                                      >  2:2311356/1‑150 (MQ=255)
                                                                                                 cggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAaca                                                      >  2:991479/1‑150 (MQ=255)
                                                                                                  ggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAacac                                                     >  2:2062688/1‑150 (MQ=255)
                                                                                                  ggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAacac                                                     >  1:1788793/1‑150 (MQ=255)
                                                                                                  ggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCACGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTACAACGACGCTGATGTCAacac                                                     >  1:930084/1‑150 (MQ=255)
                                                                                                      cAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATcc                                                 <  2:4103309/150‑1 (MQ=255)
                                                                                                       aGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCa                                                >  1:1889094/1‑150 (MQ=255)
                                                                                                        gTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAg                                               <  2:4406428/150‑1 (MQ=255)
                                                                                                        gTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAg                                               >  1:1207021/1‑150 (MQ=255)
                                                                                                          tACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCt                                             <  1:3167407/150‑1 (MQ=255)
                                                                                                           aCGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTc                                            >  2:2900941/1‑150 (MQ=255)
                                                                                                             gccgccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCg                                          <  1:1698526/150‑1 (MQ=255)
                                                                                                             gccgccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCg                                          >  1:2660680/1‑150 (MQ=255)
                                                                                                              ccgccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGt                                         >  2:2344354/1‑150 (MQ=255)
                                                                                                              ccgccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGt                                         >  1:2591116/1‑150 (MQ=255)
                                                                                                               cgccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGtt                                        <  2:3305455/150‑1 (MQ=255)
                                                                                                                gccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTa                                       <  1:4131517/150‑1 (MQ=255)
                                                                                                                 ccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTaa                                      >  2:3626383/1‑150 (MQ=255)
                                                                                                                    ttGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCg                                   <  1:3980899/150‑1 (MQ=255)
                                                                                                                     tGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGt                                  >  1:382882/1‑150 (MQ=255)
                                                                                                                       tAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTgg                                >  2:154108/1‑150 (MQ=255)
                                                                                                                        aGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTggg                               <  1:1255606/150‑1 (MQ=255)
                                                                                                                        aGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTggg                               <  1:2620009/150‑1 (MQ=255)
                                                                                                                         ggAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATACCTCCCCACCCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGa                              <  2:1131340/150‑1 (MQ=255)
                                                                                                                          gAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGaa                             <  1:3397866/150‑1 (MQ=255)
                                                                                                                           aaTTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACTTCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGaaa                            <  1:1068625/150‑1 (MQ=255)
                                                                                                                           aaTTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGaaa                            >  1:3376524/1‑150 (MQ=255)
                                                                                                                           aaTTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGaaa                            >  1:712494/1‑150 (MQ=255)
                                                                                                                           aaTTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTGACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGaaa                            <  2:1177487/150‑1 (MQ=255)
                                                                                                                            aTTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAAc                           <  1:3694155/150‑1 (MQ=255)
                                                                                                                                aTCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTa                       >  2:1808202/1‑150 (MQ=255)
                                                                                                                                aTCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTa                       >  2:3780649/1‑150 (MQ=255)
                                                                                                                                aTCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCCCCGTTAAGCGTGGGAAACAGTa                       <  1:3419447/150‑1 (MQ=255)
                                                                                                                                 tCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTaa                      <  2:2394134/150‑1 (MQ=255)
                                                                                                                                  cGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAg                     <  2:4457828/150‑1 (MQ=255)
                                                                                                                                  cGCCGGATGCAAGGGTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAg                     <  1:480080/150‑1 (MQ=255)
                                                                                                                                   gCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGc                    <  2:4260998/150‑1 (MQ=255)
                                                                                                                                    ccGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCa                   >  2:3319603/1‑150 (MQ=255)
                                                                                                                                     cGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCAc                  >  1:4252629/1‑150 (MQ=255)
                                                                                                                                      ggATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCAct                 <  1:3459027/150‑1 (MQ=255)
                                                                                                                                       gATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCActc                <  1:4545092/150‑1 (MQ=255)
                                                                                                                                        aTGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCActct               >  2:907340/1‑150 (MQ=255)
                                                                                                                                          gCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGa             >  1:2746856/1‑150 (MQ=255)
                                                                                                                                           cAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGAc            >  2:205691/1‑150 (MQ=255)
                                                                                                                                                 tCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGATAg      <  1:1085842/150‑1 (MQ=255)
                                                                                                                                                     gccgcATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGATAGtatt  <  2:2661532/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGATAGTATT  >  CP009273/1810470‑1810768

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: