Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,508,782 C→T 5.5% intergenic (+38/‑60) xanP → / → tusD xanthine permease/sulfur transfer protein complex subunit TusD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,508,7820CT 5.5% 827.8 / 24.3 263intergenic (+38/‑60)xanP/tusDxanthine permease/sulfur transfer protein complex subunit TusD
Reads supporting (aligned to +/- strand):  ref base C (137/110);  new base T (6/8);  total (145/118)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 267 total aligned reads displayed.

GCGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGGG  >  NC_002947/4508634‑4508931
                                                                                                                                                    |                                                                                                                                                     
  gcACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGcacc                                                                                                                                                    >  1:1048337/1‑149 (MQ=255)
    acaTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAg                                                                                                                                                  >  2:3561504/1‑150 (MQ=255)
      aTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcg                                                                                                                                                >  2:2662393/1‑150 (MQ=255)
      aTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTGGCACaaaag                                                                                                                                                >  1:881408/1‑147 (MQ=255)
      aTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCATGCACAAccg                                                                                                                                                >  1:2194364/1‑150 (MQ=255)
       tGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgg                                                                                                                                               <  1:368800/150‑1 (MQ=255)
         ccGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggc                                                                                                                                             >  2:5512150/1‑150 (MQ=255)
          cggcggCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcg                                                                                                                                            >  2:6792164/1‑150 (MQ=255)
             cggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGACCCTGGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcgg                                                                                                                                         <  2:2462528/150‑1 (MQ=255)
             cggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgcgg                                                                                                                                         <  1:7190826/150‑1 (MQ=255)
             cggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcgg                                                                                                                                         <  2:6211856/150‑1 (MQ=255)
                 gATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCtt                                                                                                                                     >  1:6657944/1‑150 (MQ=255)
                    tCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCttttt                                                                                                                                  <  2:1666968/150‑1 (MQ=255)
                      gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTg                                                                                                                                >  2:4791923/1‑150 (MQ=255)
                      gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTg                                                                                                                                <  2:2227542/150‑1 (MQ=255)
                      gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTg                                                                                                                                >  1:2895150/1‑150 (MQ=255)
                       cAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGc                                                                                                                               >  1:4324378/1‑150 (MQ=255)
                        aaTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCt                                                                                                                              <  2:3067048/150‑1 (MQ=255)
                         aTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCATTTTTGCtt                                                                                                                             <  1:3935941/150‑1 (MQ=255)
                            tgCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTAt                                                                                                                          <  1:2012326/150‑1 (MQ=255)
                             gCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATc                                                                                                                         <  2:4133087/150‑1 (MQ=255)
                              cTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCa                                                                                                                        >  2:353204/1‑150 (MQ=255)
                               tGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGACCGCGCCCGCTTGCACAATCGGGCGCGGGGTTTTTTGCTTTATAAt                                                                                                                       >  1:593809/1‑150 (MQ=255)
                                 gAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATgg                                                                                                                     <  2:2765550/150‑1 (MQ=255)
                                 gAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATgg                                                                                                                     <  1:6961577/150‑1 (MQ=255)
                                  aGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCATTTTTGCTTTATCATGGc                                                                                                                    <  2:20857/150‑1 (MQ=255)
                                      gggTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCa                                                                                                                <  1:4107563/150‑1 (MQ=255)
                                       ggTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAt                                                                                                               >  1:3072282/1‑150 (MQ=255)
                                             gCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCttt                                                                                                         <  2:5268006/150‑1 (MQ=255)
                                               tACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGc                                                                                                       >  2:1766220/1‑150 (MQ=255)
                                               tACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTGCACAACCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAATCCTTTGc                                                                                                       >  1:161920/1‑150 (MQ=255)
                                                aCCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                      >  2:1984826/1‑150 (MQ=255)
                                                 ccGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCAt                                                                                                     <  1:4036006/150‑1 (MQ=255)
                                                   ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   <  2:877399/150‑1 (MQ=255)
                                                   ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   <  1:4667727/150‑1 (MQ=255)
                                                    gggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTTCACAAGCGGGCGCGGGCTTTTTTGGTTTATCATGGCACCAATCCTTTGCATGAg                                                                                                  >  1:6322477/1‑150 (MQ=255)
                                                      gggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGtt                                                                                                <  1:2061815/150‑1 (MQ=255)
                                                         aTTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGc                                                                                             <  2:6389756/150‑1 (MQ=255)
                                                         aTTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGc                                                                                             <  1:974070/150‑1 (MQ=255)
                                                                 gcTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAAt                                                                                     <  2:2134364/150‑1 (MQ=255)
                                                                 gcTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAAt                                                                                     <  2:6881163/150‑1 (MQ=255)
                                                                   tGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTc                                                                                   <  1:6126602/150‑1 (MQ=255)
                                                                   tGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTACTTTGCATGAGTTATCCATTAAATTc                                                                                   >  1:1530039/1‑150 (MQ=255)
                                                                         tGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTGGCTATc                                                                             >  2:3739829/1‑150 (MQ=255)
                                                                           aaCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATcgc                                                                           >  2:5949132/1‑150 (MQ=255)
                                                                           aaCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTATTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATcgc                                                                           <  1:5389982/150‑1 (MQ=255)
                                                                            aCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATcgcg                                                                          <  1:2881988/150‑1 (MQ=255)
                                                                             cTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGGAATTCGCTATCGCgg                                                                         >  2:5813703/1‑150 (MQ=255)
                                                                             cTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCgg                                                                         >  2:5673835/1‑150 (MQ=255)
                                                                              tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt                                                                        <  1:1931270/150‑1 (MQ=255)
                                                                              tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt                                                                        <  1:562951/150‑1 (MQ=255)
                                                                              tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCACCAGCCGGCGCGGGCTGTTTTGCTTTATCATGGCACCCTACCTTTGCAGGAGTTAGCCATGAAATAAGCGATCGCGgt                                                                        >  2:1511122/1‑150 (MQ=255)
                                                                              tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCCCCCCCGTGCACAAGCCGGGCCGGGCTTTTTTGCTTTTTAAAGGCACCATTACTTTGGATGAGTTATCCCAGAAAATCGCAATCGAggg                                                                        >  1:1012814/1‑149 (MQ=255)
                                                                              tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCCCACCCCCGGGCACAACCCGGCCCGGGCGTTTTTGGTTTATAAAGGCACCAAGCCAATGCCTGGGGTAGCCCAGAAAATAGCTATCGCggg                                                                        >  1:3160083/1‑149 (MQ=255)
                                                                                gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGGTTTTTTGGTTTATCATGGCACCAATACTTTGCATGAGTTTGCCCTGAAATTGGCGATAGCGgtgt                                                                      >  1:6849559/1‑150 (MQ=255)
                                                                                 ttCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTATTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGtt                                                                     <  1:6472126/150‑1 (MQ=255)
                                                                                   cTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtct                                                                   >  2:5031981/1‑150 (MQ=255)
                                                                                   cTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtct                                                                   >  2:7421692/1‑150 (MQ=255)
                                                                                    ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTACTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtctc                                                                  >  2:3140686/1‑150 (MQ=255)
                                                                                     tGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   >  2:6128642/1‑116 (MQ=255)
                                                                                      gCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTccc                                                                <  1:3258770/150‑1 (MQ=255)
                                                                                      gCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTccc                                                                <  1:5224148/150‑1 (MQ=255)
                                                                                        cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgg                                                                                                                                               <  2:6183603/69‑1 (MQ=255)
                                                                                        cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgg                                                                                                                                               >  1:6183603/1‑69 (MQ=255)
                                                                                        cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTAGCTATAGCGGTGTTCTCCCCg                                                              >  1:1139126/1‑150 (MQ=255)
                                                                                        cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCGTTGCATGTGTTAGCCCTGAAAATCGCTATCGCGGTGTTCTCCCCg                                                              >  2:645736/1‑150 (MQ=255)
                                                                                             ggCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTGTTTTGCGTTATCATGGCACCCTTCCTGTGCATGAGGTAGCCATGACATTCGCTCTCGCGGTGTTCTCCCCGGCCCa                                                         >  2:3898780/1‑150 (MQ=255)
                                                                                                aaGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCAATCGCGGTGTTCTCCCCGGCCCAAgc                                                      >  1:1170139/1‑150 (MQ=255)
                                                                                                aaGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCGTTGCATGAGTTAGCCCTGAAATTCGCTATCGCGGTGTTCTCCCCCGCCCCTgc                                                      >  2:968472/1‑150 (MQ=255)
                                                                                                    aGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGccc                                                  <  2:5613791/150‑1 (MQ=255)
                                                                                                          cTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCAAGAAATTCGCAATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTcg                                            >  2:7200084/1‑150 (MQ=255)
                                                                                                           ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCAAGCGCCCTCTTcgc                                           >  2:5345323/1‑150 (MQ=255)
                                                                                                              aaCCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAc                                        >  2:2905525/1‑150 (MQ=255)
                                                                                                              aaCCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTGTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTCGCGCGGTGTTCTCCCCGGCCCATGCGCCCGCGTCGCGag                                        >  2:7383761/1‑149 (MQ=255)
                                                                                                                ccTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAcgc                                      <  2:5765764/150‑1 (MQ=255)
                                                                                                                 cTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTTCTTTATAATGGCACCATTACTTTTCATGAGTTAGCCATGAAATTAGCAATCGCGGTGTTCTCCCCGGCCCATTCGCCCTCTTCGCGAcgcg                                     >  2:367487/1‑150 (MQ=255)
                                                                                                                  ttGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACgcgc                                    <  2:3336615/150‑1 (MQ=255)
                                                                                                                     agaCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATTAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCt                                 <  1:3608913/150‑1 (MQ=255)
                                                                                                                       aCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAAgcgcagggc                               <  1:7280063/150‑6 (MQ=255)
                                                                                                                            ctGAGAGGAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTCTCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCGATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTc                          <  2:3980213/149‑1 (MQ=255)
                                                                                                                               aaaaaaGCCCGCGCCCGCTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGcc                       >  1:482151/6‑150 (MQ=255)
                                                                                                                                  gTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAAGGCAACAATCCATTGCATGAGTTATCCAAGAAATTCGCGATCGCGGTGTTCTCCCCCGCCCATGCGCCCACTTAGCGACGCGCCCATCGCTTCGCCGAg                    >  2:650214/1‑150 (MQ=255)
                                                                                                                                   tAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAgg                   <  1:1285289/150‑1 (MQ=255)
                                                                                                                                     gcccgcgcccgcTTGCACAAGCGGGCGCGGGCATTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCg                 <  1:6726297/150‑1 (MQ=255)
                                                                                                                                       ccgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt               <  2:482151/150‑1 (MQ=255)
                                                                                                                                       ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCggccagca               >  1:6572126/1‑143 (MQ=255)
                                                                                                                                       ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt               >  1:5752900/1‑150 (MQ=255)
                                                                                                                                       ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTAGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt               >  1:796970/1‑150 (MQ=255)
                                                                                                                                        cgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              >  1:5362227/1‑150 (MQ=255)
                                                                                                                                         gcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCGCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGc             >  2:4656451/1‑150 (MQ=255)
                                                                                                                                         gcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGc             >  2:2408467/1‑150 (MQ=255)
                                                                                                                                           gcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTg           >  1:5785269/1‑150 (MQ=255)
                                                                                                                                           gcccgcTTGCACAAGCGGGCGCGGGCTGTTGTGCTTTATCATGGCACCCGTCCGTTGCCGGGGTTAGCCATGACCGTCGCTATCGCGGTGTTCTCCCCGGCCCCTGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGGGCTg           >  2:2806749/1‑150 (MQ=255)
                                                                                                                                             ccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCCTGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGc         >  2:5171568/1‑150 (MQ=255)
                                                                                                                                             ccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGc         >  1:1684649/1‑150 (MQ=255)
                                                                                                                                                 ttGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggc     >  2:1929980/1‑150 (MQ=255)
                                                                                                                                                  tGCACAAGCGGGCGCGGGCTTATTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcg    <  1:4937894/150‑1 (MQ=255)
                                                                                                                                                   gCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg   <  2:1206996/150‑1 (MQ=255)
                                                                                                                                                   gCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg   >  2:4364312/1‑150 (MQ=255)
                                                                                                                                                    cacaAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCggg  >  2:1847638/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGGG  >  NC_002947/4508634‑4508931

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: