Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,508,782 C→T 7.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 7.8% 522.2 / 25.2 179intergenic (+38/‑60)xanP/tusDxanthine permease/sulfur transfer protein complex subunit TusD
Reads supporting (aligned to +/- strand):  ref base C (103/62);  new base T (7/7);  total (110/69)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

Only 100 of 187 total aligned reads displayed.

GGCGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGGG  >  NC_002947/4508633‑4508931
                                                                                                                                                     |                                                                                                                                                     
 gcgcACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGca                                                                                                                                                      >  1:3682418/1‑150 (MQ=255)
   gcACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGcaca                                                                                                                                                    >  1:2427446/1‑150 (MQ=255)
      caTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgc                                                                                                                                                 <  1:4203201/150‑1 (MQ=255)
       aTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcg                                                                                                                                                <  2:3268309/150‑1 (MQ=255)
         gCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGTAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcggg                                                                                                                                              >  2:5118115/1‑150 (MQ=255)
             gcggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcg                                                                                                                                          >  1:4045097/1‑150 (MQ=255)
              cggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgcgg                                                                                                                                         <  1:4417682/150‑1 (MQ=255)
              cggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgcgg                                                                                                                                         <  1:2088303/150‑1 (MQ=255)
               ggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcggg                                                                                                                                        >  1:2188340/1‑150 (MQ=255)
                gcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcgggc                                                                                                                                       <  2:4409972/150‑1 (MQ=255)
                  gATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTGCACAAGCGGGCGCGGGCtt                                                                                                                                     >  1:5163787/1‑150 (MQ=255)
                    ttCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCtttt                                                                                                                                   >  1:2986187/1‑150 (MQ=255)
                       gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTg                                                                                                                                >  2:1171527/1‑150 (MQ=255)
                       gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTg                                                                                                                                >  2:3999801/1‑150 (MQ=255)
                        cAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgg                                                                                                                                               >  2:3930084/1‑134 (MQ=255)
                        cAATGGGCTGGAGTCGGGTGTGGCTACCGGGGGGTTTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGc                                                                                                                               <  1:575395/150‑1 (MQ=255)
                              gCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATc                                                                                                                         <  2:471611/150‑1 (MQ=255)
                               cTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACACGCGGGCGCGGGCTTTTTTGCTTTATCa                                                                                                                        >  2:3869068/1‑150 (MQ=255)
                                tGGATTAGGTTGTGGCTACAGGGGGTATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCAGCGACCGCTTGCACAAGCGGGCGCGGTCTTTTTTGCTTTATAAt                                                                                                                       >  2:1200964/1‑150 (MQ=255)
                                tGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAt                                                                                                                       <  1:4879653/150‑1 (MQ=255)
                                 ggAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTTCACAAGCGGGCGCGGGCTTTTTTTCTTTATAATg                                                                                                                      >  1:2049605/1‑150 (MQ=255)
                                  gAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATgg                                                                                                                     >  2:1119322/1‑150 (MQ=255)
                                     tCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCAc                                                                                                                  >  2:2871713/1‑150 (MQ=255)
                                       gggTGTGGCTACCGGGGGGATTTGTGCGCAGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCATGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAAGGCACCa                                                                                                                >  2:2798616/1‑150 (MQ=255)
                                        ggTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGGGCCCCCGCGGGCGCGGGCTTTTTTGGTTTATAATGGGACCAt                                                                                                               >  1:488571/1‑150 (MQ=255)
                                         gtgtGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCaat                                                                                                              >  1:2663031/1‑148 (MQ=255)
                                         gtgtGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAtt                                                                                                              >  1:884099/1‑150 (MQ=255)
                                          tgtgGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTATTTTGCATTATCATGGCACCATTc                                                                                                             <  2:3734356/150‑1 (MQ=255)
                                              gCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCttt                                                                                                         <  1:3965981/150‑1 (MQ=255)
                                              gCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCttt                                                                                                         <  2:674549/150‑1 (MQ=255)
                                              gCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCATATTTGCTTTATCATGGCAccattccatt                                                                                                         <  2:3886595/150‑1 (MQ=255)
                                               cTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTg                                                                                                        <  2:3217109/150‑1 (MQ=255)
                                               cTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTg                                                                                                        >  2:1979880/1‑150 (MQ=255)
                                                tACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAACCGGGCGCGGGCTGTTTTGCATTATAAATGCACCATTACATTga                                                                                                       >  1:3381212/1‑149 (MQ=255)
                                                  ccGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCAt                                                                                                     <  2:4913685/150‑1 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTTCTTTATAATGGCACCATTCCTTTTCAAGa                                                                                                   >  1:598379/1‑150 (MQ=255)
                                                    ggggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTGTGCATGa                                                                                                   >  2:4283138/1‑150 (MQ=255)
                                                      ggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGt                                                                                                 <  1:4394017/150‑1 (MQ=255)
                                                      ggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAACCGGGCCCGGGCTTGTTTGCTTTATCAAGGCACCAGTACTTTGCATTAGt                                                                                                 >  1:3298212/1‑150 (MQ=255)
                                                       gggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAATCCTTTGCATGAGtt                                                                                                >  1:2332896/1‑150 (MQ=255)
                                                        ggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGTAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTa                                                                                               <  1:5118115/150‑1 (MQ=255)
                                                        ggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTa                                                                                               <  1:5258136/150‑1 (MQ=255)
                                                          aTTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCCCGCCCCCGTGCACAAACGGGCGCGGGCGTTTTTGCTTTATAAATGCACCAATCCATTGCAAGAGTTAGc                                                                                             >  1:4714785/1‑150 (MQ=255)
                                                            ttGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATTGCACCAATAATTTTCATTAGTTAGCCa                                                                                           >  1:4883602/1‑150 (MQ=255)
                                                               tgCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTTCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATGGCACCATTACTTTGCATGAGTTAGCCATGa                                                                                        >  1:3574223/1‑150 (MQ=255)
                                                                  gcTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCATTCACAAACGGGCGCGGGCTTTTTTTGTTTATAATTGGACCAATACTTTTCATTTGTTAGGCATGAAAt                                                                                     >  1:1398158/1‑150 (MQ=255)
                                                                    tGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTc                                                                                   <  2:564843/150‑1 (MQ=255)
                                                                    tGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGGAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTc                                                                                   <  1:1983114/150‑1 (MQ=255)
                                                                         cTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATGGCACCATTCCGTTGCAGGAGGTAGCCATGAAATTCGCGAt                                                                              >  2:1916054/1‑150 (MQ=255)
                                                                            aaCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTAGCTATagc                                                                           >  1:5029436/1‑150 (MQ=255)
                                                                               tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATTGCACCATTACTTTGCATGAGTTATCCATGAAATTAGCTATCGCGgt                                                                        >  1:4501165/1‑150 (MQ=255)
                                                                               tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCCCCCCCGTGCACAAGCCGGCGCGGGCGTTTTTGGGTTATCAAGGCACCATTCCGTTGCCAGAGTTTGCCAAGAAAATAGCTATcgccgg                                                                        >  1:1073227/1‑147 (MQ=255)
                                                                                 gTTCTTGCCGGAGGGCCAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgtgt                                                                      >  1:4612737/1‑150 (MQ=255)
                                                                                 gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCAGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgtgt                                                                      <  2:4957329/150‑1 (MQ=255)
                                                                                    cTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGGTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCGATCGCGGTGTtct                                                                   >  1:2747518/1‑150 (MQ=255)
                                                                                     ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtctc                                                                  <  1:3922526/150‑1 (MQ=255)
                                                                                     ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtctc                                                                  >  1:3638613/1‑150 (MQ=255)
                                                                                     ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCCCGCCCCCGTGCACAACCGGGCGCGGGCTTTTTTTGTTTAAAAAGGCACCAATACATTTCATGAGTTAGCCCAGAAAATAGCAATAGCGGTGTtctc                                                                  >  1:758458/1‑150 (MQ=255)
                                                                                        ccGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGGTTTATCATGGGACCAATCCATTGCATGAGTTAGCCATGGAATTAGCTGTAGCGGGGTTCTcccc                                                               >  1:3807430/1‑150 (MQ=255)
                                                                                         cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCATGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATTGCACCAATCCATTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCg                                                              >  1:3641504/1‑150 (MQ=255)
                                                                                          ggAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCGTTTTTGCTTTATAATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGGGTTCTACCCgg                                                             >  1:4360123/1‑150 (MQ=255)
                                                                                           gAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCGTGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGc                                                            <  1:2637579/150‑1 (MQ=255)
                                                                                                cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCGTGCACAAGCGGGCGCGGGCTTGTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATg                                                       >  2:3561523/1‑150 (MQ=255)
                                                                                                cAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTTCACAAGCGGGCGCGGGGGTTTTTGGGTTATAAAGGCACCAAAACATTGCAAGAGTTAGCCAATAAATTAGCGATAGCGGTGTTAACCCCCGccccag                                                       >  1:1463910/1‑147 (MQ=255)
                                                                                                 aaGGAGCAGGCTTGAACCTTGAGACTGACTGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATgc                                                      <  1:5572499/150‑1 (MQ=255)
                                                                                                 aaGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATgc                                                      <  2:2653765/150‑1 (MQ=255)
                                                                                                 aaGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCCGCCCATgc                                                      >  1:3727456/1‑150 (MQ=255)
                                                                                                 aaGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCAGGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATgc                                                      >  2:2946060/1‑150 (MQ=255)
                                                                                                    gAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGcc                                                   <  2:884099/150‑1 (MQ=255)
                                                                                                        aGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCAATACTTTGCATGAGTTAGCCAAGAAATTCGCTATCGCGGTGTTATCCCCGGCCCATGCGCCcact                                               >  1:1827769/1‑150 (MQ=255)
                                                                                                            ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCACAACCGGGCGCGGGCATTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCCTGAAAATCGCGATCGCGGTGTTCACCCCGGCCCATGCGCCCTCGTcgc                                           >  2:2094398/1‑150 (MQ=255)
                                                                                                              gAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTGATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGa                                         <  2:1998446/150‑1 (MQ=255)
                                                                                                               aaCCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAc                                        <  1:114231/150‑1 (MQ=255)
                                                                                                               aaCCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAc                                        <  1:1863374/150‑1 (MQ=255)
                                                                                                                  cTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAATCGGGCGCGGGCTTTTTTGCTTTAAAATGGCACCATTCCTTTGCATGCGTTAGCCATGAAAATCGCTATCGCGGTGTTCAACCCCGCCCATGCGCCCTCTTCGCGAcgcg                                     >  2:127978/1‑150 (MQ=255)
                                                                                                                  cTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCGTTATCATGGCACCATTCCTTTTCATGAGTTATCCAAGATATTGGCTATAGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAcgcg                                     >  2:1591073/1‑150 (MQ=255)
                                                                                                                  cTTGAGACTGAATGAGAGTAGCCCCCGCCCGCATTCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCCTGGAATTCGCTATCGCGGTGGTCTCCCCGGCCCAAGCGCCCGCGTCGCGAcgcg                                     >  2:127924/1‑150 (MQ=255)
                                                                                                                   ttGAGACTGAATGAGAGTAGCCCCCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTTCTTTATAATGGCACCATTACTTTGCATGAGTTAGCCATGAAATTAGCTATCGCGGTGTTCTACCCGGCCCAAGCGCCCTCTTCGCGACgcgc                                    >  1:3152764/1‑150 (MQ=255)
                                                                                                                          tGAATGAGAGTAGCCCGCGCCCGCGTGCACAAGCCGGCGCGGGCTTTTTTTCTTTATAATGGCACCAATCCTTTGCAAGAGTTAGCCATTAAATTCGCTATCCCGGTGTTCTCCCCCGCCCATGCGCCCTCATCGCGACGCGCCCAgcgc                             >  2:521017/1‑150 (MQ=255)
                                                                                                                                       cccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCgg                >  2:309173/1‑150 (MQ=255)
                                                                                                                                       cccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCgg                >  2:245843/1‑150 (MQ=255)
                                                                                                                                       cccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCgg                >  1:3375349/1‑150 (MQ=255)
                                                                                                                                        ccgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt               >  1:2853438/1‑150 (MQ=255)
                                                                                                                                        ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATAATGGCACCATTCCTTTGAATGAGTTAGCCATGAAATACGCTATCGCGGTGTTCTCACCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt               >  2:724788/1‑150 (MQ=255)
                                                                                                                                         cgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                      >  2:2386316/1‑62 (MQ=255)
                                                                                                                                         cgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa                                                                                                      <  1:2386316/62‑1 (MQ=255)
                                                                                                                                         cgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              >  2:1441247/1‑150 (MQ=255)
                                                                                                                                          gcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGc             >  1:3188949/1‑150 (MQ=255)
                                                                                                                                          gcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGc             >  2:3644321/1‑150 (MQ=255)
                                                                                                                                              ccgcTTGCACAAGCGGGCGCGGGCTGTATTGCTTTATCATGGCACCATTCCTTTGCAAGAGTTATCCATGAAATTCGCAATCGCGGTGTACTCCCCGGACAATGCGCCCACTTCGCGACGCGCCCTGCGCTTCCCCGAGGCGGTGCTGGc         >  1:1705388/1‑150 (MQ=255)
                                                                                                                                               cgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGcc        >  1:5575259/1‑150 (MQ=255)
                                                                                                                                                gcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCgg                                                                         >  1:188963/1‑84 (MQ=255)
                                                                                                                                                gcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCgg                                                                         <  2:188963/84‑1 (MQ=255)
                                                                                                                                                gcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcg       >  1:3518065/1‑150 (MQ=255)
                                                                                                                                                gcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcg       <  1:2285780/150‑1 (MQ=255)
                                                                                                                                                   tGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcg    <  2:1398158/150‑1 (MQ=255)
                                                                                                                                                    gCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg   <  1:2983996/150‑1 (MQ=255)
                                                                                                                                                    gCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAGGGCACCATTCCTTGGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCGGCGCGACGCGCCCTGCGCGTCGCCGAGGCGGTGCTGGCcggcgg   <  2:2599612/150‑1 (MQ=255)
                                                                                                                                                     cacaAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCggg  <  1:1557019/150‑1 (MQ=255)
                                                                                                                                                     cacaAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCggg  >  2:869135/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 < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: