Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,810,619 G→T 38.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,6190GT38.7% 304.8 / 242.1 229intergenic (+234/+24)katE/chbGcatalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base G (53/86);  new base T (45/44);  total (99/130)
Fisher's exact test for biased strand distribution p-value = 7.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGATAGTATT  >  CP009273/1810470‑1810768
                                                                                                                                                     |                                                                                                                                                     
gCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcg                                                                                                                                                       >  2:2176062/1‑150 (MQ=255)
  cTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCActcc                                                                                                                                                     <  2:2681571/150‑1 (MQ=255)
  cTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgcc                                                                                                                                                     >  2:2083471/1‑150 (MQ=255)
    aCAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCAcgccgc                                                                                                                                                   >  2:2205984/1‑150 (MQ=255)
     cAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCa                                                                                                                                                  <  2:4167572/150‑1 (MQ=255)
          cTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTgg                                                                                                                                             >  2:112172/1‑150 (MQ=255)
          cTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTgg                                                                                                                                             <  1:199726/150‑1 (MQ=255)
                aaaTTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACa                                                                                                                                       <  1:1166448/150‑1 (MQ=255)
                aaaTTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACa                                                                                                                                       <  2:1461120/150‑1 (MQ=255)
                 aaTTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACAt                                                                                                                                      >  2:3887422/1‑150 (MQ=255)
                  aTTCAATATATTGCAGGAAACACGTAGGCCTGATACGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACCTc                                                                                                                                     <  2:3225658/150‑1 (MQ=255)
                  aTTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATc                                                                                                                                     <  2:2990536/150‑1 (MQ=255)
                   ttCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATcc                                                                                                                                    >  1:357072/1‑150 (MQ=255)
                    tCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCt                                                                                                                                   >  2:1909818/1‑150 (MQ=255)
                        tatatTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACt                                                                                                                               <  2:1634900/150‑1 (MQ=255)
                          tatTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTa                                                                                                                             >  2:3432785/1‑150 (MQ=255)
                            ttGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGATTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTaca                                                                                                                           <  2:2492315/150‑1 (MQ=255)
                             tGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTacac                                                                                                                          >  1:3455189/1‑150 (MQ=255)
                             tGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTacac                                                                                                                          >  1:4130865/1‑150 (MQ=255)
                              gCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTAcaca                                                                                                                         >  1:2617781/1‑150 (MQ=255)
                                   aaaCACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCg                                                                                                                    <  2:1770915/150‑1 (MQ=255)
                                    aaCACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGa                                                                                                                   <  2:2014102/150‑1 (MQ=255)
                                       acGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATaa                                                                                                                <  1:4438971/150‑1 (MQ=255)
                                               cTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCaa                                                                                                        >  1:3974849/1‑150 (MQ=255)
                                                tGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAAc                                                                                                       <  2:683492/150‑1 (MQ=255)
                                                    aagcgaagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGAt                                                                                                   >  2:3226664/1‑150 (MQ=255)
                                                      gcgaagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATaa                                                                                                 <  2:655450/150‑1 (MQ=255)
                                                        gaagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAAcc                                                                                               <  2:862226/150‑1 (MQ=255)
                                                        gaagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAAcc                                                                                               >  1:3759514/1‑150 (MQ=255)
                                                         aagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCa                                                                                              <  2:3346080/150‑1 (MQ=255)
                                                         aagcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCa                                                                                              >  2:977377/1‑150 (MQ=255)
                                                          agcgCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCAc                                                                                             <  1:880306/150‑1 (MQ=255)
                                                              cATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTg                                                                                         <  2:1006891/150‑1 (MQ=255)
                                                                tCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAg                                                                                       <  2:4519611/150‑1 (MQ=255)
                                                                 cAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGc                                                                                      <  2:2265374/150‑1 (MQ=255)
                                                                 cAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGc                                                                                      <  1:4512002/150‑1 (MQ=255)
                                                                  aGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCg                                                                                     >  1:1272422/1‑150 (MQ=255)
                                                                   ggCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGa                                                                                    >  2:2962523/1‑150 (MQ=255)
                                                                      aGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAg                                                                                 >  2:2690223/1‑150 (MQ=255)
                                                                       gTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGc                                                                                >  2:3912377/1‑150 (MQ=255)
                                                                         tttGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCAc                                                                              <  2:3722/150‑1 (MQ=255)
                                                                         tttGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCAc                                                                              <  2:4461439/150‑1 (MQ=255)
                                                                          ttGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCAcc                                                                             <  2:4368085/150‑1 (MQ=255)
                                                                             cGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGTTAGCACCttt                                                                          <  2:717545/150‑1 (MQ=255)
                                                                             cGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCttt                                                                          <  2:567864/150‑1 (MQ=255)
                                                                                 tGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAac                                                                      <  1:633638/150‑1 (MQ=255)
                                                                                  gTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAacg                                                                     >  1:3802890/1‑150 (MQ=255)
                                                                                    catcagTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAacgac                                                                   <  1:2042656/147‑1 (MQ=255)
                                                                                     agcagTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAacgacg                                                                  >  2:4103412/1‑150 (MQ=255)
                                                                                       cagTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCt                                                                >  1:2697989/1‑150 (MQ=255)
                                                                                        agTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACCCTGAGCGATAGCACCTTTCAACGACGCTg                                                               >  2:722221/1‑150 (MQ=255)
                                                                                        agTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTg                                                               <  1:2599043/150‑1 (MQ=255)
                                                                                          tctcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGAt                                                             <  2:3545847/150‑1 (MQ=255)
                                                                                            tcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGt                                                           >  1:1383159/1‑150 (MQ=255)
                                                                                            tcAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGt                                                           <  2:4089083/150‑1 (MQ=255)
                                                                                             cAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCGACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTc                                                          >  1:1598634/1‑150 (MQ=255)
                                                                                             cAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTc                                                          >  2:971615/1‑150 (MQ=255)
                                                                                                gcggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACACCCCACCCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAac                                                       <  2:2354048/150‑1 (MQ=255)
                                                                                                  ggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAacac                                                     >  1:1950022/1‑150 (MQ=255)
                                                                                                  ggcgCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAacac                                                     >  1:3474942/1‑150 (MQ=255)
                                                                                                    cgcAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACAt                                                   >  2:3949644/1‑150 (MQ=255)
                                                                                                      cAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATcc                                                 <  1:1320717/150‑1 (MQ=255)
                                                                                                      cAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATcc                                                 <  1:4354605/150‑1 (MQ=255)
                                                                                                       aGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCa                                                >  1:2436381/1‑150 (MQ=255)
                                                                                                       aGTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCa                                                >  1:2045086/1‑150 (MQ=255)
                                                                                                        gTTACGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAg                                               <  2:560712/150‑1 (MQ=255)
                                                                                                           aCGCCGCCTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTc                                            >  1:4119567/1‑150 (MQ=255)
                                                                                                               cgccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACCTCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGtt                                        <  2:2137362/150‑1 (MQ=255)
                                                                                                                gccTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTa                                       >  1:1432626/1‑150 (MQ=255)
                                                                                                                  cTTTGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAg                                     <  1:1735776/150‑1 (MQ=255)
                                                                                                                   tttGTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGc                                    <  2:1039932/150‑1 (MQ=255)
                                                                                                                    ttGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCg                                   <  1:1016259/150‑1 (MQ=255)
                                                                                                                     tGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCACCCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCGGATGTCAACACATCCAGCTCCGTTAAGCGt                                  <  1:2222256/150‑1 (MQ=255)
                                                                                                                     tGTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTGACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGt                                  <  1:4602000/150‑1 (MQ=255)
                                                                                                                      gTAGGAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTg                                 <  1:3425100/150‑1 (MQ=255)
                                                                                                                      gTAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTg                                 >  1:2009687/1‑150 (MQ=255)
                                                                                                                       tAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTgg                                >  1:1162856/1‑150 (MQ=255)
                                                                                                                       tAGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTgg                                <  2:2142551/150‑1 (MQ=255)
                                                                                                                        aGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTggg                               >  2:1217450/1‑150 (MQ=255)
                                                                                                                        aGGAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTggg                               <  2:137786/150‑1 (MQ=255)
                                                                                                                          gAATTAATCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGaa                             <  1:1341739/150‑1 (MQ=255)
                                                                                                                          gAATTAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACGCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGaa                             <  2:3630530/150‑1 (MQ=255)
                                                                                                                             ttAATCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACa                          <  1:1847648/150‑1 (MQ=255)
                                                                                                                               aaTCGCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGATGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGt                        <  1:1530686/150‑1 (MQ=255)
                                                                                                                                 tCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTaa                      <  1:3556612/150‑1 (MQ=255)
                                                                                                                                 tCGCCGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTaa                      <  1:3823192/150‑1 (MQ=255)
                                                                                                                                   gCCGGATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGc                    <  2:184713/150‑1 (MQ=255)
                                                                                                                                    ccGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCa                   <  2:2380690/150‑1 (MQ=255)
                                                                                                                                    ccGGATGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCa                   >  2:1359346/1‑150 (MQ=255)
                                                                                                                                     cGGATGCAAGGTTAACTCCGCATCTGCCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGAGATAGCAAATATCAACGAAGCTGAAGTCAACACATCAAGCTCCGTTAAGCGTGGGAAACATAAAGCAc                  >  2:1566076/1‑150 (MQ=255)
                                                                                                                                      ggATGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCAct                 >  1:1554217/1‑150 (MQ=255)
                                                                                                                                        aTGCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCActct               >  1:448742/1‑150 (MQ=255)
                                                                                                                                         tGCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTg              <  1:2946426/150‑1 (MQ=255)
                                                                                                                                          gCAAGGTTCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGa             <  1:516652/150‑1 (MQ=255)
                                                                                                                                          gCAAGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGa             <  2:1329153/150‑1 (MQ=255)
                                                                                                                                           cAAGGTACACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATCACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGAc            <  1:3142075/150‑1 (MQ=255)
                                                                                                                                            aaGGTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACg           >  2:2809117/1‑150 (MQ=255)
                                                                                                                                               gTTCACGCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATACCCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGAt        <  2:4219512/150‑1 (MQ=255)
                                                                                                                                                 tCACTCCGCATCTGGCAAACATCCTCACTTACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGATAg      <  1:237364/150‑1 (MQ=255)
                                                                                                                                                     gccgcATCTGGCAAACATCCTCACTCACACATCCCGATAACTCCCCAACCGATAACCACGCTGAGCGATAGCACCTTTCAACGACGCTGATGTCAACACATCCAGCTCCGTTAAGCGTGGGAAACAGTAAGCACTCTGACGGATAGtatt  <  2:554510/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: