Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,636 A→T 5.2% intergenic (+117/‑91) pyk → / → PP_4302 pyruvate kinase/urea transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,898,6360AT 5.2% 707.5 / 19.2 210intergenic (+117/‑91)pyk/PP_4302pyruvate kinase/urea transporter
Reads supporting (aligned to +/- strand):  ref base A (106/93);  new base T (8/3);  total (114/96)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 217 total aligned reads displayed.

CTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4898488‑4898784
                                                                                                                                                    |                                                                                                                                                    
cTAACCCGCTGCGTATCGAGAGTTTGATCTGTGTTTGCGGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGTAATGATGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCaaa                                                                                                                                                     <  2:1303326/150‑1 (MQ=38)
cTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCaaa                                                                                                                                                     >  2:4363716/1‑150 (MQ=255)
 tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAATGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTACCAGTGGGAGCGGCCGTGTGTCGCGGTTGGAGGGCAAAg                                                                                                                                                    >  1:5548632/1‑150 (MQ=255)
 tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAg                                                                                                                                                    >  2:510747/1‑150 (MQ=255)
 tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCCCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCAAAg                                                                                                                                                    >  2:2442266/1‑150 (MQ=255)
       gCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGGCACAAAGCCCCGCCCCCAGGGCCCCCCCGGGACAACGGGATGCAGTCCCAGTGGGGGCGGCCGTGTGTGGCGATTTGTGGGGAAAGcccacc                                                                                                                                              >  1:254239/1‑150 (MQ=255)
            ggATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGGCACAAGGCCCCTCCCACAGCGGCACCCGGGCCCAAGGGACAGCACTACCGGTGGGGGCGGCCATGTGTGGCGGTTGGAGGGGCAATcccacccagca                                                                                                                                         >  1:5185091/1‑146 (MQ=255)
            ggATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCAGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCGCCCGGaa                                                                                                                                         >  1:996950/1‑150 (MQ=255)
               tCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAc                                                                                                                                      <  2:3231428/150‑1 (MQ=255)
                    aCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctc                                                                                                                                 >  2:5499780/1‑150 (MQ=255)
                     cTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAACGCCCCGCCCACCGGGCCACCCGGGCCCCACGGCAAGCAGTCCCGGGGGGCGCGGCCGTGTGTCGCGCTTGGAGGGCAAACCCCACCCGGAAAAACGcaca                                                                                                                                >  1:2029786/1‑147 (MQ=255)
                         gATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCCCGCGGGCACAAGGGAATGCAGGCCCGGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCCAAGCCCGCCCTGAAAACCGCGCTAccc                                                                                                                            >  1:545387/1‑150 (MQ=38)
                                     cggggggCGCGTTGCGCCCCTTTCGCCACACACGCCCGCTCCCACCGGGCCCCGGTGGCACACGGGAATGCAGTCCCTGTGGGTGGGGCCTTGTGCCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGaaaa                                                                                                                <  1:1282827/148‑1 (MQ=255)
                                       aggggCGCAACGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCGCCCGGAATACCGCTCTACCCGACCGCGGAAAAcc                                                                                                              >  2:410695/2‑150 (MQ=39)
                                         gggCGCTTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACCGGGCCCCGCTGGAACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                                                            <  1:2612962/150‑1 (MQ=255)
                                          ggCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGATGGCCGGCCCTGTGGGTGCGGCCTTGTGCCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                           <  1:1550139/150‑1 (MQ=255)
                                          ggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                           <  2:5408204/150‑1 (MQ=255)
                                           gcgcTTGGCCCCCCTTGCGCCACACCCGCCCGCCCCCACCGGGCCACGCGGGCACACGGGATTGCTGACCCTGGGGGAGCGCCCTTTGGACGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAgg                                                                                                          <  1:3020727/150‑1 (MQ=255)
                                           gcgcGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGa                                                                                                                          <  2:5611806/134‑1 (MQ=255)
                                           gcgcGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGa                                                                                                                          >  1:5611806/1‑134 (MQ=255)
                                                       ccTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGcc                                                                                              >  2:2622941/1‑150 (MQ=255)
                                                         aatcgcgacacaaggccgctcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                                                            >  1:5509602‑M2/82‑150 (MQ=255)
                                                           tCGCTACACAATGCCGGTCCCACAGGCCCACGCTGTCACAAGGGATTTCAGTCCCTGTGGGTTCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGACTACCGCTCTTCCCGACAGCGTAAATTGGTGTCCCATAGAGGCCCTcc                                                                                          <  2:1518524/150‑1 (MQ=255)
                                                            CGCGACACACGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCg                                                                                         >  1:1055577‑M1/1‑150 (MQ=255)
                                                            CGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCg                                                                                         >  2:6803033‑M1/1‑150 (MQ=255)
                                                             gcgacacaaggccgctcccacagggactgcattcccgtgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGc                                                                                        >  1:837830‑M2/78‑150 (MQ=255)
                                                                  ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcccc                                                                                   <  1:3661883‑M1/150‑1 (MQ=255)
                                                                    AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  2:4783050‑M1/150‑1 (MQ=255)
                                                                    AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  1:442906‑M1/150‑1 (MQ=255)
                                                                    AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  2:754090‑M1/150‑1 (MQ=255)
                                                                    AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  2:6340526‑M1/150‑1 (MQ=255)
                                                                         CGCTCCCCCAGGGCCACGCTGGCACAAGGGAATGCCGTCCCTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                                                            >  2:642844‑M1/1‑150 (MQ=255)
                                                                         CGCTCCCACCGGGCCACGCTGGCACATGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                                                            >  2:4060079‑M1/1‑150 (MQ=255)
                                                                            aCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         <  1:6120102‑M1/149‑1 (MQ=255)
                                                                               CACAGGGCCGCGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATGGGAGGGCtttgccctccaatc                                                                                                                                          >  1:4523461‑M1/1‑78 (MQ=255)
                                                                               CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      <  2:5218009‑M1/150‑1 (MQ=255)
                                                                                cccGGGCCACGCTGGCACACGGGATTGCAGTCCGTGTGGTAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                                                                     >  2:1767370‑M1/4‑150 (MQ=255)
                                                                                ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                                                                     <  2:1708864‑M1/150‑1 (MQ=255)
                                                                                  AGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactga                                                                   >  1:925412‑M1/1‑150 (MQ=255)
                                                                                    GGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccct                                                                                                                                                >  1:1079029‑M1/1‑71 (MQ=255)
                                                                                      ccACGCTGGCACAAGGGTATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGc                                                               <  2:1858363/150‑1 (MQ=255)
                                                                                          gCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCa                                                           >  1:2178612/1‑150 (MQ=255)
                                                                                          gCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCa                                                           <  2:5873682/150‑1 (MQ=255)
                                                                                           cTGGCACATGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCAt                                                          <  2:2387546/150‑1 (MQ=255)
                                                                                           cTGCCACACGGGAATGCAGTCCCTGGGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCAt                                                          <  1:1899919/150‑1 (MQ=255)
                                                                                              gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       <  2:5279986/150‑1 (MQ=255)
                                                                                                acaAGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:6238417/150‑1 (MQ=255)
                                                                                                acaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     >  2:2359468/1‑150 (MQ=255)
                                                                                                  atgggATTGCAGTCCCTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAcc                                                   <  1:3080127/148‑1 (MQ=255)
                                                                                                   aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  <  2:2178612/150‑1 (MQ=255)
                                                                                                   aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  >  2:925153/1‑150 (MQ=255)
                                                                                                   aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  >  1:4797331/1‑150 (MQ=255)
                                                                                                      gAATGCAGTCCCTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaaa                                               <  2:1590812/150‑1 (MQ=255)
                                                                                                      gAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaaa                                               >  1:3977781/1‑150 (MQ=255)
                                                                                                          cgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAAttt                                           >  1:2104523‑M2/33‑150 (MQ=255)
                                                                                                           ctttCCCTGGGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                          <  1:2477343/147‑1 (MQ=255)
                                                                                                           cAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                          <  2:3977668/150‑1 (MQ=255)
                                                                                                           cAGTCCCTGTGGGAGCGGAATTTTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                          <  2:877821/150‑1 (MQ=255)
                                                                                                             gTCCGTGTGGGCGCGGCGTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGt                                        <  1:25697/150‑1 (MQ=255)
                                                                                                             gTCCCTGTGGTAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGt                                        <  2:3101541/150‑1 (MQ=255)
                                                                                                             gTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCACAAATTTCGt                                        >  1:4467467/1‑150 (MQ=255)
                                                                                                             gTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGt                                        >  1:6537054/1‑150 (MQ=255)
                                                                                                             gTCCATGTGGGAGCGGCCTTGTGTCGCGATAGGATGGCAAAGCCCTCCCTGACTACCGCTCTCCCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAGCAACAAATTCgg                                        >  2:6779874/1‑149 (MQ=38)
                                                                                                              tCCCTGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACGGACTGCCCCATGTACACCAAAAATTTCGTc                                       <  2:1121361/150‑1 (MQ=255)
                                                                                                              tCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTc                                       >  2:4559467/1‑150 (MQ=255)
                                                                                                              tCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTc                                       >  1:6649612/1‑150 (MQ=255)
                                                                                                                cctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCaa                                     >  2:1225070‑M2/27‑150 (MQ=255)
                                                                                                                 ctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAAc                                    >  2:4859026‑M2/26‑150 (MQ=255)
                                                                                                                 ctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAAc                                    <  2:3661883‑M2/125‑1 (MQ=255)
                                                                                                                     gggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGt                                >  2:3949394‑M2/22‑150 (MQ=255)
                                                                                                                            gccttgggtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGa                         <  1:4305520‑M2/136‑1 (MQ=255)
                                                                                                                             ccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGAc                        <  2:4889756‑M2/137‑1 (MQ=255)
                                                                                                                              cttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACt                       >  1:4598731‑M2/13‑150 (MQ=255)
                                                                                                                                tgtgtcgcgaTTGGAGGGCGAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  2:6293231‑M2/140‑1 (MQ=255)
                                                                                                                                tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  2:2029786‑M2/140‑1 (MQ=255)
                                                                                                                                tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  2:5187634‑M2/140‑1 (MQ=255)
                                                                                                                                tgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCCCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                     <  1:5499780‑M2/140‑1 (MQ=255)
                                                                                                                                    tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  2:3920371‑M2/7‑150 (MQ=255)
                                                                                                                                    tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  2:3361301‑M2/7‑150 (MQ=255)
                                                                                                                                    tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:1710606‑M2/7‑150 (MQ=255)
                                                                                                                                    tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:6619376‑M2/7‑150 (MQ=255)
                                                                                                                                    tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:293796‑M2/7‑150 (MQ=255)
                                                                                                                                     cgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGTCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCAc                <  1:1509109‑M2/145‑1 (MQ=255)
                                                                                                                                       ttcagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACgg              >  1:3321766/6‑150 (MQ=255)
                                                                                                                                        tcagggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGc             >  1:2286206/5‑150 (MQ=255)
                                                                                                                                        gaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGc             >  1:5101290‑M2/3‑150 (MQ=255)
                                                                                                                                        gaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGc             <  1:4910729‑M2/148‑1 (MQ=255)
                                                                                                                                         aTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCg            <  1:754090‑M2/149‑1 (MQ=255)
                                                                                                                                          agggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           >  1:2460531/3‑150 (MQ=255)
                                                                                                                                          agggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           >  2:4769353/3‑150 (MQ=255)
                                                                                                                                          agggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           >  2:6513801/3‑150 (MQ=255)
                                                                                                                                           gggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGtt          >  2:893555/2‑150 (MQ=255)
                                                                                                                                            ggAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTa         <  1:4168515/150‑1 (MQ=255)
                                                                                                                                              aGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACt       <  1:2359468/150‑1 (MQ=255)
                                                                                                                                               gggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTc      >  2:2108013/1‑150 (MQ=255)
                                                                                                                                               gggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTc      >  2:2786660/1‑150 (MQ=255)
                                                                                                                                                ggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGt                                >  1:3703575/1‑123 (MQ=255)
                                                                                                                                                ggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGt                                <  2:3703575/123‑1 (MQ=255)
                                                                                                                                                 gCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCaa    >  1:5803059/1‑150 (MQ=255)
                                                                                                                                                   aaaGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACg  >  1:276733/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4898488‑4898784

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: