Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,508,782 C→T 5.8% 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.8% 763.7 / 22.9 245intergenic (+38/‑60)xanP/tusDxanthine permease/sulfur transfer protein complex subunit TusD
Reads supporting (aligned to +/- strand):  ref base C (112/118);  new base T (6/8);  total (119/126)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 257 total aligned reads displayed.

GGCGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGGG  >  NC_002947/4508633‑4508931
                                                                                                                                                     |                                                                                                                                                     
      caTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgc                                                                                                                                                 <  2:6999088/150‑1 (MQ=255)
          ccGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggc                                                                                                                                             <  2:3019119/150‑1 (MQ=255)
           cggcggCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcg                                                                                                                                            >  1:2743355/1‑150 (MQ=255)
            ggcggcGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgc                                                                                                                                           <  1:7383997/150‑1 (MQ=255)
            ggcggcGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgc                                                                                                                                           <  2:5566340/150‑1 (MQ=255)
             gcggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcg                                                                                                                                          <  2:1648079/150‑1 (MQ=255)
              cggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgcgg                                                                                                                                         <  1:3411731/150‑1 (MQ=255)
              cggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcgg                                                                                                                                         <  2:7892994/150‑1 (MQ=255)
                      cGCAATGTGCTGGCGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGACTGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCtttttt                                                                                                                                 <  1:7286319/150‑1 (MQ=255)
                       gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAACGGGCGCGGGCTTTTTTg                                                                                                                                >  1:3730923/1‑150 (MQ=255)
                         aaTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCt                                                                                                                              <  2:1262203/150‑1 (MQ=255)
                          aTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCtt                                                                                                                             <  1:8541321/150‑1 (MQ=255)
                          aTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCtt                                                                                                                             <  2:8383822/150‑1 (MQ=255)
                           tgtgCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCttt                                                                                                                            >  1:1935191/1‑150 (MQ=255)
                              gCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATc                                                                                                                         >  2:4105905/1‑150 (MQ=255)
                                  gAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATgg                                                                                                                     <  2:7966237/150‑1 (MQ=255)
                                  gAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATgg                                                                                                                     <  1:3620550/150‑1 (MQ=255)
                                  gAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAACCGGGCGCGGGCTTTTTTGCGTTATAATgg                                                                                                                     >  1:5592287/1‑150 (MQ=255)
                                   aGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGc                                                                                                                    >  2:2750017/1‑150 (MQ=255)
                                    gTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCa                                                                                                                   <  1:5845309/150‑1 (MQ=255)
                                     tCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTTTAATGGCAc                                                                                                                  >  1:6241561/1‑150 (MQ=255)
                                        ggTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAt                                                                                                               >  2:2508715/1‑150 (MQ=255)
                                          tgtgGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTc                                                                                                             <  1:3365590/150‑1 (MQ=255)
                                            tgGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTGTTTGCTTTATCATGGCACCATTCCt                                                                                                           >  2:3547101/1‑150 (MQ=255)
                                              gCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCttt                                                                                                         <  1:954829/150‑1 (MQ=255)
                                                tACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGc                                                                                                       >  2:2894117/1‑150 (MQ=255)
                                                tACCGGGGGGATTTGCGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGACTGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCCCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGc                                                                                                       >  2:2239447/1‑150 (MQ=255)
                                                 aCCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                      <  1:8140263/150‑1 (MQ=255)
                                                 aCCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                      <  2:2607368/150‑1 (MQ=255)
                                                 aCCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                      <  2:4579397/150‑1 (MQ=255)
                                                  ccGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAGTCCTTTGCAt                                                                                                     >  2:8181580/1‑150 (MQ=255)
                                                  ccGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATTGCAACATTACTTTGCAt                                                                                                     >  1:8460665/1‑150 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAAGTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGCGAGTAGCCCGCGCCCGGTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGTCACCATTCCTTTGCATGa                                                                                                   <  1:80702/150‑1 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   <  1:6008813/150‑1 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   <  1:4468949/150‑1 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   <  1:7144408/150‑1 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   <  1:4819339/150‑1 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGa                                                                                                   <  1:80380/150‑1 (MQ=255)
                                                    ggggggATGTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCGTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCGTTTTTGCTTTATCATGGCACCATTCCTTTGCATgg                                                                                                   >  2:1225290/1‑149 (MQ=255)
                                                     gggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATGGCACCATTCCTTTGCATGAg                                                                                                  >  2:6863868/1‑150 (MQ=255)
                                                        ggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTGTTTGCTTTATCATGGCACCATTCCTTTGCCTGAGGTa                                                                                               >  2:3222609/1‑150 (MQ=255)
                                                         gATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCATTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAg                                                                                              <  1:5449119/150‑1 (MQ=255)
                                                         gATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCAAAAACGGGCGCGGGCTTTATTTCTTTATCATGGAACCATTACTTTGCAAGAGTTAg                                                                                              >  2:811451/1‑150 (MQ=255)
                                                           tttGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGcc                                                                                            <  1:3926421/150‑1 (MQ=255)
                                                           tttGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGAAAGGAGCAGGCTTGAACCTTTAGACTGAATGAGAGTAGCCCGGGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGcc                                                                                            <  1:2174675/150‑1 (MQ=255)
                                                             tgtgCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCAt                                                                                          <  1:2849739/150‑1 (MQ=255)
                                                                 cgcTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCATTGCATGAGTTAGCCCTGtaa                                                                                      >  1:6156647/1‑150 (MQ=255)
                                                                    tGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTc                                                                                   >  1:7058547/1‑150 (MQ=255)
                                                                     ggTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCg                                                                                  <  2:1401804/150‑1 (MQ=255)
                                                                      gTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGc                                                                                 <  2:3325865/150‑1 (MQ=255)
                                                                      gTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAATCCTTTGCATGAGTTAGCCATGAAATTCGc                                                                                 >  2:6266557/1‑150 (MQ=255)
                                                                        gCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTa                                                                               >  2:6447550/1‑150 (MQ=255)
                                                                          tGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATc                                                                             >  2:2079563/1‑150 (MQ=255)
                                                                              cTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCgg                                                                         <  1:1610664/150‑1 (MQ=255)
                                                                               tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt                                                                        <  1:7167606/150‑1 (MQ=255)
                                                                               tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt                                                                        <  2:3607449/150‑1 (MQ=255)
                                                                                ggTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTTCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAGTACGTTTCATGAGTTAGCCATGAAAATAGCTATCGCgggg                                                                       >  1:4572947/1‑148 (MQ=255)
                                                                                 gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgtgt                                                                      <  1:6781077/150‑1 (MQ=255)
                                                                                 gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAACGGGCGCGGGCTTTTTTTCTTTATAAATGCACCAATCCATTTCGAGAGTTAGCCAAGAAAATAGCAATAGCGgtgt                                                                      >  1:115733/1‑150 (MQ=255)
                                                                                    cTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtct                                                                   >  2:7143633/1‑150 (MQ=255)
                                                                                         cGGAGTGCAAGGAGCAGGCTTGAACCTTGAGACTTAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCg                                                              <  1:3538341/150‑1 (MQ=255)
                                                                                         cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATTGCACCAATAAATTGCAAGAGTTAGCCAAGAAAATCGCTATCGCGGTGTTCTCCCCg                                                              >  1:3258950/1‑150 (MQ=255)
                                                                                         cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCACAAGCCGGCGCGGGCTTTTTTTCTTTATAATGGCACCATTCCATTGCATGAGTTAGCCATGAAATTAGCTATCGCGGGGTTCTCCCCg                                                              >  1:6945666/1‑150 (MQ=255)
                                                                                                cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCAt                                                        >  1:8558162/1‑149 (MQ=255)
                                                                                                cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCAt                                                        <  2:8558162/149‑1 (MQ=255)
                                                                                                cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCATTGCATGAGTTAGCCATTAAAGTCGCTATCGCGGTGTTCACCCCGGCCCATg                                                       >  2:4340412/1‑150 (MQ=255)
                                                                                                   ggAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGGACAAGCGGGCGCGGGGTTTTTTGGTTTATGATGGCACCATTACTTTTGATTTGTTATGCATGTAATTTGGTATTGCGGTGTTCTACCCGGCCCATgcgc                                                    >  2:852189/1‑150 (MQ=255)
                                                                                                    gAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGcc                                                   <  1:2215597/150‑1 (MQ=255)
                                                                                                    gAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCATTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGGTCTCCCCGGCCCATGCGcc                                                   >  2:5523351/1‑150 (MQ=255)
                                                                                                          gCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCATGCACAACCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAATCCTTTGCATGGGTTAGCCATGAAAATCGCGATCGCGGTGTTCTCCCCGGCCCATGCGcccgcgtc                                             >  1:5782102/1‑145 (MQ=255)
                                                                                                          gCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCATGCACAAACGGGCGCGGGCTTTTTTGGTTTATAAATGCACCAATACATTTCATGAGTTATCCAAGAAAATAGCTATAGCGGTGTTGTACCCCGCCCAAGCGcccacata                                             >  1:7564967/1‑145 (MQ=255)
                                                                                                          gCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATTGCACCATTCCTTTGCATGAGGTAGCCATGAAATTAGCTATCGCGGTGTTCTCCCCCGCCCAAGCGcccgcgtc                                             >  1:1204471/1‑145 (MQ=255)
                                                                                                                  cTTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTGCACAAGCCGGCGCGGGCTTTTTTGCGTTATAAAGGCACCAATACATTGCAAGAGTTAGCCAATAAATTCGCGATCGCGGTGTTATACCCCGCCCCTGCGCCCACGTCGCGAcccg                                     >  1:4149579/1‑150 (MQ=255)
                                                                                                                   ttGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGGTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTGGCTATCGCGGTGTTGTCCCCGGCCCATGCGCCCTCTTCGCGACgcgc                                    >  2:3470721/1‑150 (MQ=255)
                                                                                                                   ttGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCACAAAACGGCGCGGGCTTTTTTGCTTTATAATGGCACCATTACATTGCATGAGTTAGCCATTAAATTAGCAATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACgcgc                                    >  1:2063550/1‑150 (MQ=255)
                                                                                                                      agaCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCt                                 <  2:6035251/150‑1 (MQ=255)
                                                                                                                       gaCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTg                                <  1:3611066/150‑1 (MQ=255)
                                                                                                                       gaCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTg                                <  1:4938179/150‑1 (MQ=255)
                                                                                                                        aCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTgc                               <  1:4376050/150‑1 (MQ=255)
                                                                                                                         cTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTgcg                              <  2:4700365/150‑1 (MQ=255)
                                                                                                                         cTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAAGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTgcg                              >  2:3275941/1‑150 (MQ=255)
                                                                                                                          tGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAAATCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTgcgc                             >  2:1915673/1‑150 (MQ=255)
                                                                                                                                agagTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCATTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGcc                       <  2:5200941/150‑1 (MQ=255)
                                                                                                                                  cggaGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGa                     <  2:8361865/147‑1 (MQ=255)
                                                                                                                                  agTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGa                     <  2:7378234/150‑1 (MQ=255)
                                                                                                                                  agTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAAGCGCAGGGCGCGTCGcgaa                     <  1:8332192/150‑4 (MQ=255)
                                                                                                                                    tAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAgg                   <  2:2119972/150‑1 (MQ=255)
                                                                                                                                    tAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAgg                   <  1:3040687/150‑1 (MQ=255)
                                                                                                                                    aaGCCCGCGCCCGCTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAgg                   >  1:4989018/2‑150 (MQ=255)
                                                                                                                                      gcccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCg                 <  2:4887732/150‑1 (MQ=255)
                                                                                                                                       cccgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCgg                >  1:1393675/1‑150 (MQ=255)
                                                                                                                                        ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt               >  1:6370693/1‑150 (MQ=255)
                                                                                                                                         cgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              >  2:6299596/1‑150 (MQ=255)
                                                                                                                                         cgcgcccgcTTGCACACGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              <  1:5222644/150‑1 (MQ=255)
                                                                                                                                         cgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              >  2:4643532/1‑150 (MQ=255)
                                                                                                                                            gcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTg           <  1:2894117/150‑1 (MQ=255)
                                                                                                                                              ccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGc         >  2:6221748/1‑150 (MQ=255)
                                                                                                                                                    gCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg   <  2:3759854/150‑1 (MQ=255)
                                                                                                                                                    gCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg   <  2:5908759/150‑1 (MQ=255)
                                                                                                                                                     cacaAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCggg  >  2:8699287/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGGG  >  NC_002947/4508633‑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: