Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,610 C→G 8.2% intergenic (+91/‑117) 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,6100CG 8.2% 410.1 / 2.3 136intergenic (+91/‑117)pyk/PP_4302pyruvate kinase/urea transporter
Reads supporting (aligned to +/- strand):  ref base C (65/58);  new base G (5/8);  total (70/66)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

GCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGT  >  NC_002947/4898461‑4898746
                                                                                                                                                     |                                                                                                                                        
ggttgatgggggagcggtgggcctctttgggcctcggttttccgcggtcgggtagagcggtattcagggagggctttgccctccaaTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGc                                                                                                                                          >  1:2552220‑M1/87‑150 (MQ=255)
      tgggggagcggtgggcctctttgggcctcggttttccgcggtcgggtagagcggtattcagggagggctttgccctccaaTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCtt                                                                                                                                    >  2:5396255‑M1/81‑150 (MQ=255)
               ggtgggcctctttgggcctcggttttccgcggtcgggtagagcggtattcagggagggctttgccctccaaTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGtcgcgg                                                                                                                           >  1:3232144‑M1/72‑149 (MQ=255)
                    ggTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCCCAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGa                                                                                                                      >  2:4443104/1‑150 (MQ=255)
                      tACGACTACCACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAgg                                                                                                                    <  2:814537/150‑1 (MQ=255)
                      tACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAgg                                                                                                                    <  2:5127111/150‑1 (MQ=255)
                       aCGACTAACCCGCTGCGGATCGAGACTTGTATCTGAGATTGCCGGGGGGGCGTTGCGCCCGTTTCGCGACACACGGCCGCTCCCACCGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTTTGTCGCGATTGGAggg                                                                                                                   <  2:1255062/150‑1 (MQ=255)
                       aCGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGTCAAAACGCCCCGCCCCCAGGGCCCCCCGGGCCCAAGGGGCTGCAGTCCCAGGGGGGGCGGCCGTGTGTGGCGCTTgggggg                                                                                                                   >  1:1607600/1‑150 (MQ=255)
                          aCTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGGCACAAGGCCCCGCCCACACGGCCCCGCAGGCACAAAGGGAAGCCAACCCAGTGGGGGCGGCCATGTGTGGCGATTTGGGGGCaa                                                                                                                >  1:989792/1‑150 (MQ=255)
                            taccccGCTGCGGATCGAGACTTTGAGCTGAGTTTGCGGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCGCCCCCCGGGCCACGCTGGCACACGGGAATGCAGGCCCGGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAg                                                                                                              <  2:1936905/145‑1 (MQ=255)
                            tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGTCAAAATGCCGCTCCCCAAGGGGCACGCTGGCAGAAGGGGATTCAGTCCCTGTGGGTGGGGCCTTGTGTCGCGATTTGTGGGCAAAg                                                                                                              >  1:1240646/1‑150 (MQ=255)
                            tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGAATGGAGGGCAAAg                                                                                                              >  2:3094292/1‑150 (MQ=255)
                            tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGAAACAAGGCCGCACCCAAACGGCCAACCAGGCACCAGGGGATGCGGGCCCAGTGGGCGCGGCCGTGGGGCGGGGGTGGGGGGCCAAg                                                                                                              >  1:1446870/1‑150 (MQ=255)
                                  gCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGccctcc                                                                                                        >  2:4609571/1‑150 (MQ=255)
                                  gCTGCGGATCGAGAATTTGATCTGAGATGGCCGGGGGCGCGTTGCGCCCCTTTCGCGAAACACGGCCGCTCCCACCGGGCCACGCTGGCAAACGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGACGCGATTGGAGGGCAAAGccctcc                                                                                                        <  1:2191804/150‑1 (MQ=255)
                                          tCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACCCGGGCACAAGGGAAAGCAGTCCCGGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCAGAATAc                                                                                                >  2:5182080/1‑150 (MQ=39)
                                           cGAGACTTTTATCTGAGATTGCGGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACTGGGCCACGCTGGCACATGGGATTGCATTCCCTGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAcc                                                                                               <  1:1748259/150‑1 (MQ=255)
                                            gagaCTTTGATGTGTGTGTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACCGGGCCACGCTGGCACAAGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCg                                                                                              <  2:4241474/150‑1 (MQ=38)
                                             agaCTTAGATCTGAGAATGCCGGGGGCGCGTTGCGCCCCTAACGCGAAAAAAAGCCGCACCCACACGGCCACCCAGGCACCAGGGAATACAGTCCCGGTGGGCGCGGCCATGTGAACCGATTGGAGGGAAAATCCCACCCGGAAAAcccc                                                                                             >  2:500736/1‑148 (MQ=255)
                                               aCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctc                                                                                           <  2:1713658/150‑1 (MQ=255)
                                                cTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctct                                                                                          <  2:1316077/150‑1 (MQ=255)
                                                 tttGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTa                                                                                         >  2:4931921/1‑150 (MQ=255)
                                                 tttGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCCCTACCACACGGCCACGCTGGCACAAGGGAATGCAGTCCCAGTGGGCTCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAACGCTCTa                                                                                         >  1:4560176/1‑150 (MQ=255)
                                                 tgtGATCTGAGATAGCCGGGGGCGCGTTGCGCCCCTTTCTCGTCACAAGGCCGCTCCCAAACGGCCACGCTGGCACAAGGGAATTCAGTCACTGTGGGAGCGGCCTTGTGTCGCGAGTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTa                                                                                         >  2:2736708/3‑150 (MQ=38)
                                                      tCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCCCGCAGGCACAAGGGAATGCAGTCCCTGTGGGCGCCGCCGTTTGTCGCGAATGGAGGGCAAAGCCCTCCCTGAATACCGCGCTAccccc                                                                                    >  2:1471385/1‑148 (MQ=38)
                                                        tGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCACTACCCGAcc                                                                                  >  2:4559436/1‑150 (MQ=255)
                                                        tGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAAAGCAGTCCCTGTGGGAGCGGCCATGTGTCGCGATTGGAGGGCCAAGCCCACCCGGAATAACGCTCTACCCGAcc                                                                                  >  1:4659046/1‑150 (MQ=38)
                                                          agaTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACcgc                                                                                <  2:4333942/150‑1 (MQ=255)
                                                           gaTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAAAACCGCTCTACCCGACcgcg                                                                               >  2:3172846/1‑150 (MQ=255)
                                                            aTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCgg                                                                              <  2:5537249/150‑1 (MQ=255)
                                                            aTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCAAAGCCCTCCCGGAATACCGCTCTACCCGACCGCgg                                                                              >  2:4348724/1‑150 (MQ=255)
                                                                 cGGGGGCGGGTTGCGCCCCTTACGCGACACCCGGCCGCTCCCCCCGGGCCCCGCTGGCACATGGGAATGCAGTCCCTGTGGGCGGGGCCTTGTGCCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAc                                                                         <  1:4491714/150‑1 (MQ=255)
                                                                  gggggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGACTGCAGTCCCTGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCCAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAcc                                                                        >  2:2728337/1‑150 (MQ=255)
                                                                   ggggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCg                                                                       <  2:5627882/150‑1 (MQ=255)
                                                                   ggggCGCAACGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCAGGCACAAGGGAATGCAGTCCCAGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCAGAAAACCGCTCTAACCGACAGCGGAAAACCg                                                                       >  1:5606572/1‑150 (MQ=38)
                                                                    gggCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACCGGGCCACGCTGGCACATGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                      <  2:103519/150‑1 (MQ=255)
                                                                    gggCGCGTTGCGCCCCTTTCGCGACACACGCCCGCTCCCACTGGGCCACGCTGGAACAAGGGAATGCTGTCCCTGTGGGCGCGGCCTTGTTACGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                      <  1:2943988/150‑1 (MQ=255)
                                                                    gggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACACGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                      <  2:873176/150‑1 (MQ=255)
                                                                    gggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                      <  2:2302406/150‑1 (MQ=255)
                                                                    gggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                      <  2:1627752/150‑1 (MQ=255)
                                                                    gggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCGCCCACCGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                      <  1:1798159/150‑1 (MQ=255)
                                                                    gggCGCGTTGCCCCCCTTTCGCGACACACGGCCGCTCCCACAGGGCCACGCTGCCACAAGGGAATGATGTCCCTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                      <  1:4891522/150‑1 (MQ=255)
                                                                     ggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                     <  2:5095459/150‑1 (MQ=255)
                                                                     ggCGCGTTGCCCCCCTTTCGCCACACCCGGCCGCCCCCCCCGGGCCACGCTGGAACACGGTATTGCAGCCCCTGTGGGAGGGGCCTTGTGACGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                     <  1:2067478/150‑1 (MQ=255)
                                                                      gcgcGTTGCGCCCATTTCGCGACACACGGCCGCGCCCACCGGGCCACGCTGGCACATGGTAATGCAGTCATTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAgg                                                                    <  1:3116547/150‑1 (MQ=255)
                                                                        gcgTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCGGGCACAAGGGAATGCAGTCCCGGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCCACCGCGGAAAACCGAgggc                                                                  >  2:5280558/1‑148 (MQ=255)
                                                                            ggcgcCCTTTTCGCGCCACACGGCCGCCCCCCCCGGGCCCCGCTGGAACATGGTAATGCTGTCCGTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                              <  1:4756384/149‑1 (MQ=255)
                                                                            cgcgcCCCTTGCGAGACCCAAGCCCCCACCCCAAGGCCCACCCTGGCACACCGGATTGCAGCCCCTCGGGCAGCGGCCTTGTGTCGCGATTGGAGGGCAAACCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                              <  1:4880349/149‑1 (MQ=255)
                                                                             gcgcCCCTTTCGCGACACACGGCCGCTCCCACTGGGCCCCGTGGGCAAACGGGATTGCTGGCCTTTTGGGAGGGGCCTTGTTTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAag                                                             <  1:161506/150‑1 (MQ=255)
                                                                               gcCCTTTGCGCGACACCCGGCGGCTCCCACGGGGCCCGGGTGGCAAATGGTAAGGCTGGCCGTGTGGTGGGGGCTTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAagcg                                                           <  1:3031487/150‑3 (MQ=255)
                                                                                ccccTTTCGCGACACAATGCCGCTCCCACAGGGCCACGCTGGGACAAGGGAATGCAGTCCCTGTGGGTGCGGCCGTGTGTGGCGATTGGAGGGCAAAGCCCTCCCTGAATAACGCTATACCCCTCCGCGGAAAACCGAGGCCCAAATAgg                                                          >  2:3080654/1‑150 (MQ=255)
                                                                                  ccTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGTAGCGGCCTTGTGTCGCGATTGGAGGGCAAGGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGcc                                                        <  2:2275478/150‑1 (MQ=255)
                                                                                  ccTTTCGCGACACAAGGCCGCACCCACAGGGCCACGCGGGCACAAGGGAAAGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCACCCTGAATACCGCTCTACCCGAACGCGGAAAACCGAGGCCCAAAGAGGcc                                                        >  2:1078677/1‑150 (MQ=39)
                                                                                    tttCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                      >  1:4533373/1‑150 (MQ=255)
                                                                                          cACACAGGGCCGCGCCCAGCGGGCCACGCTGTCACTCGGTATTGCCGCCCCTGTGGGTGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTc                                                >  1:2351079‑M1/2‑150 (MQ=255)
                                                                                             ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGa                                                                                                                      >  2:295929‑M1/1‑77 (MQ=255)
                                                                                             ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGa                                                                                                                      <  1:295929‑M1/77‑1 (MQ=255)
                                                                                             ACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTcccc                                             <  2:680874‑M1/150‑1 (MQ=255)
                                                                                               acGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                           >  2:2949429‑M1/3‑150 (MQ=255)
                                                                                                 GGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTTTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                         >  2:2491545‑M1/1‑150 (MQ=255)
                                                                                                 GGCCGCTCCCACAGGGCCACGCTGGCAAAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                         >  2:2491519‑M1/1‑150 (MQ=255)
                                                                                                  cccGCTCCCACAGGGCCACGCTGGCACAAGGCAATGCAGTCCCTGTGGGAGCGGCCTTGTGACGCGATTGGAGGGCAAAGCCCTCCCAGAATACCGCTCTACCCGACCGCGGAAAACCGACGCCCAAAGAGGCCCACCGCTCCCCCATCa                                        <  2:1354359/149‑1 (MQ=255)
                                                                                                    CGCTCCCCCTGGGCCACGCGGGAACATGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                      >  2:2347904‑M1/1‑150 (MQ=255)
                                                                                                    CGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAc                                      <  1:1829968‑M1/150‑1 (MQ=255)
                                                                                                     GCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAcc                                     <  2:3242222‑M1/150‑1 (MQ=255)
                                                                                                       TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCaaccac                                   <  2:3740679‑M1/150‑3 (MQ=255)
                                                                                                       TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                   <  1:3422522‑M1/150‑1 (MQ=255)
                                                                                                       TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                   <  1:1032711‑M1/150‑1 (MQ=255)
                                                                                                         CCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCa                                 <  2:4283805‑M1/150‑1 (MQ=255)
                                                                                                         CCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCa                                 <  1:4197554‑M1/150‑1 (MQ=255)
                                                                                                          CACAGGGCCACGCTGGCACAAGGTAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                >  2:4714025‑M1/1‑150 (MQ=255)
                                                                                                          CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                <  1:4523766‑M1/150‑1 (MQ=255)
                                                                                                           ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCACAGAGGCCCACCGCTCCCCCATCAACCTCCact                               <  2:3530506‑M1/150‑1 (MQ=255)
                                                                                                           ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                               <  2:3785641‑M1/150‑1 (MQ=255)
                                                                                                           ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                >  2:77718‑M1/1‑149 (MQ=255)
                                                                                                           ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                >  2:77701‑M1/1‑149 (MQ=255)
                                                                                                           ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                <  1:77701‑M1/149‑1 (MQ=255)
                                                                                                           ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                <  1:77718‑M1/149‑1 (MQ=255)
                                                                                                             AGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactga                             >  2:5329319‑M1/1‑150 (MQ=255)
                                                                                                              GGGCCAGGCTGGCACATGGTTTTGCTGTCCGTGGGGGTGGGGCGTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgac                            >  1:514636‑M1/1‑150 (MQ=255)
                                                                                                              GGGCCACGCTGCCACAAGGTAAGTCTGTCCGTGTTGGAGGGGCTTTTTTTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgac                            >  1:193896‑M1/1‑150 (MQ=255)
                                                                                                                   aCGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGccc                       <  2:5111945/150‑1 (MQ=255)
                                                                                                                     gCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCa                     >  1:4437145/1‑150 (MQ=255)
                                                                                                                         gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGTCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                 >  2:1427029/1‑150 (MQ=255)
                                                                                                                         gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                 <  2:3101198/150‑1 (MQ=255)
                                                                                                                         gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                 <  2:5158180/150‑1 (MQ=255)
                                                                                                                           acaTGGGAATGCAGTCCGGGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca               <  1:1585372/150‑1 (MQ=255)
                                                                                                                           acaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca               >  1:3580381/1‑150 (MQ=255)
                                                                                                                           acaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca               >  1:3345608/1‑150 (MQ=255)
                                                                                                                            caAGGGAAGGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTacac              <  1:3184329/150‑1 (MQ=255)
                                                                                                                              aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa            >  2:4621665/1‑150 (MQ=255)
                                                                                                                              aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa            >  2:1822965/1‑150 (MQ=255)
                                                                                                                               gggAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaa           >  2:1957806/1‑150 (MQ=255)
                                                                                                                                ggAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaa          >  1:5026981/1‑150 (MQ=255)
                                                                                                                                  aaTGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaaaa        >  1:3553104/1‑150 (MQ=255)
                                                                                                                                     gCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAAttt     <  2:4456602/150‑1 (MQ=255)
                                                                                                                                     gCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACAGCCCCATGTACACCAAAAAttt     >  2:2540523/1‑150 (MQ=255)
                                                                                                                                      cAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc    <  2:1607600/150‑1 (MQ=255)
                                                                                                                                       aGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCg   >  1:3523408/1‑150 (MQ=255)
                                                                                                                                       aGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCg   >  1:4566004/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
GCGTGCCTTTCGGCAAGCCGGGTACGACTAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGT  >  NC_002947/4898461‑4898746

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: