Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 499,221 C→G 7.0% intergenic (+62/‑234) uhpA → / → PP_0411 two‑component system response regulator UhpA/polyamine ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947499,2210CG 7.0% 500.3 / 11.2 165intergenic (+62/‑234)uhpA/PP_0411two‑component system response regulator UhpA/polyamine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (54/95);  new base G (2/11);  total (59/106)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

ACCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGC  >  NC_002947/499073‑499370
                                                                                                                                                        |                                                                                                                                                     
aCCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACg                                                                                                                                                      <  2:2901630/150‑1 (MQ=255)
aCCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACg                                                                                                                                                      <  2:6362076/150‑1 (MQ=255)
aCCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCCCGAACg                                                                                                                                                      >  2:6311596/1‑150 (MQ=255)
 cccACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCCTGCCCGCGCACgg                                                                                                                                                     >  2:6431975/1‑150 (MQ=255)
 cccACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACgg                                                                                                                                                     <  2:3909924/150‑1 (MQ=255)
  ccACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGACGAAACC‑‑T‑‑GTGGGAGCGGGCATGCCCGCGAAGAggcc                                                                                                                                                  <  2:5113950/150‑1 (MQ=39)
  ccACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGACGAAACC‑‑T‑‑GTGGGAGCGGGCATGCCCGCGAAGAggcc                                                                                                                                                  <  2:447519/150‑1 (MQ=39)
   gcaagcgcccgaaaatcgggggtactgagagaatgccatcatccggatggaagatttgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgaaaccT‑‑GTGGGAGCGGGCATGCCCCCCAACAggc                                                                                                                                                   <  1:5248812‑M1/29‑1 (MQ=255)
   cacaAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTTT‑‑GTGGGAGCGGGCATGCCCGCGAACGggc                                                                                                                                                   <  1:3456839/150‑1 (MQ=255)
   cacaAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCCTGCCCGCGAACGggc                                                                                                                                                   >  1:383805/1‑150 (MQ=255)
   cacaAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGggc                                                                                                                                                   <  1:1035426/150‑1 (MQ=255)
   cacaAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAACC‑‑T‑‑GTGGGAGCGGGCATGCCCGCGAAGAggccg                                                                                                                                                 <  2:963645/150‑2 (MQ=39)
        ggCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGCACGggcctggc                                                                                                                                              >  1:1559088/1‑150 (MQ=255)
           aCGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGgcctggcctg                                                                                                                                           <  2:4517065/150‑1 (MQ=255)
            cGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCctgc                                                                                                                                          <  1:5584504/150‑1 (MQ=255)
             gCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCctgct                                                                                                                                         <  1:3874079/150‑1 (MQ=255)
             gCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCAGGCctgct                                                                                                                                         >  1:5121950/1‑150 (MQ=255)
                tGATGCAGAAGTTGAAAGTGACCTCGCTGGCGGAGCGGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAt                                                                                                                                      <  2:2458767/150‑1 (MQ=255)
                tGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGggcctg                                                                                                                                                <  1:4135213/140‑1 (MQ=255)
                tGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGggcctg                                                                                                                                                >  2:4135213/1‑140 (MQ=255)
                   tGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAtctc                                                                                                                                   <  1:301148/150‑1 (MQ=255)
                   tGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAtctc                                                                                                                                   <  1:2645982/150‑1 (MQ=255)
                    gCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGGAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCAGGCCTGCTGATctct                                                                                                                                  >  2:5248812/1‑150 (MQ=255)
                    gCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATctct                                                                                                                                  >  2:4974896/1‑150 (MQ=255)
                     cAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTa                                                                                                                                 <  2:3864807/150‑1 (MQ=255)
                      ggttactgagataatgccatcatcgggatggatgttttgtgttagggatcgccctacaaggatggcgggtttggaatgcccttgcgagccaaccacgaaaccT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTAt                                                                                                                                >  1:6149358‑M1/103‑150 (MQ=255)
                      aGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTAt                                                                                                                                <  2:6013531/150‑1 (MQ=255)
                       taatttGAAAGTGAACTCGCTGGCGGAGCTGGTTAATTACGCAATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGCCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATg                                                                                                                               <  1:5754817/146‑1 (MQ=255)
                        aaGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATgg                                                                                                                              <  2:1494365/150‑1 (MQ=255)
                        aaGTTGAAAGTGAACTCGCGGGCGGTGCTGGTGAAGTACGCCTTGGAGCACAAGCTGATCTGATTGCGACATCCGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGTGCGGGCAGGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATgg                                                                                                                              <  1:4653043/150‑1 (MQ=255)
                           ttGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGcgg                                                                                                                           >  2:4368876/1‑148 (MQ=255)
                           ttGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCCCGAACCGGCCAGGCCAGCTGATCGCTATGGCTg                                                                                                                           >  2:3034466/1‑150 (MQ=255)
                            tGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGa                                                                                                                          <  1:4285849/150‑1 (MQ=255)
                                  ttaaCTCGCTGGCGGAGCTGGTTAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGc                                                                                                                    <  1:5297065/148‑1 (MQ=255)
                                  tGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCGCAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGc                                                                                                                    <  2:4204695/150‑1 (MQ=255)
                                  tGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGc                                                                                                                    <  2:784970/150‑1 (MQ=255)
                                   atgccatcatccggatggaagatttgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgaaaccT‑‑GTGGGAGCGGGCATGCCCGCGAAGAG‑‑‑‑GCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTc                                                                                                              >  2:3484695‑M1/90‑150 (MQ=255)
                                         gCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTcg                                                                                                             <  1:2588762/150‑1 (MQ=255)
                                          cTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCAAGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCACTTcgc                                                                                                            >  2:6152596/1‑150 (MQ=255)
                                           atccggatggaagatttgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgaaaccT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTcgcg                                                                                                           >  1:2137058‑M1/82‑150 (MQ=255)
                                            ggCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGTCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          <  1:1106584/150‑1 (MQ=255)
                                            ggCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          <  1:3111609/150‑1 (MQ=255)
                                            ggCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          <  2:6360023/150‑1 (MQ=255)
                                             gCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCAGGCCCCCGAACGGGCCTGGCCTGGTGATCTCGAAGGCTGAGCCGGCCGCATAGCggg                                                                                                         >  1:6224685/1‑150 (MQ=255)
                                              cggatggaagatttgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgacgccTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          >  2:2522618‑M1/79‑150 (MQ=255)
                                              cggatggaagatttgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgaagccTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          >  2:5817055‑M1/79‑150 (MQ=255)
                                              cggatggaagatttgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgaaaccT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGc                                                                                                        >  1:3204653‑M1/79‑150 (MQ=255)
                                                gcgCTGGTGAAGTACGCCATGGAGCACACGCTGATCTGATTGCGACATAGGGGTGCCCGACACTTGCCGAAGCCTAT‑‑GGGGGAGCGGGCTTGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCAt                                                                                                      <  1:144834/148‑1 (MQ=255)
                                                    tGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGccc                                                                                                  <  1:5023832/150‑1 (MQ=255)
                                                      gTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGACGCCTAT‑‑GTGGGAGCGGGCAGGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGc                                                                                                <  1:2063086/150‑1 (MQ=255)
                                                      gTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGc                                                                                                <  2:5584312/150‑1 (MQ=255)
                                                       tGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          >  1:3457857/1‑139 (MQ=255)
                                                           tgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgaaaccT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          <  1:4390974‑M1/70‑1 (MQ=255)
                                                           tgttgtagggatcgccctacaaggatggcgggattggaatgcccttgcgagccaaccacgaaaccT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          >  2:4390974‑M1/66‑135 (MQ=255)
                                                           gTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCCCGAACGGGCCGGGGCGGGTGTTGTGTATGGGTGCGCCGGCCTCTTCGCGGGCAAGACCGCGcccc                                                                                           >  1:6726774/1‑149 (MQ=255)
                                                              cGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGGATTGCCGCGAACGGGGGTGGGCTTGTGATATATATTTGTTAGCCGGCCTATTAGCGGGCATTACCGCAGCCACAg                                                                                        >  1:3289832/1‑150 (MQ=255)
                                                              cGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCAAGCGCGCGAAAGGGCCAGGCCAGCGGATCTCGAAGGCTGAGCCGGCCTCCTCGCGGGCCTGCCCGCGCCCCCAg                                                                                        >  2:4331486/1‑150 (MQ=255)
                                                               gCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCCCGCACGGGCCGGGCCAGCGGGTCTCTAAGGCTGGGCCGGCCGCGTCGCGGGCCAGCCCCCAccccccgg                                                                                       >  2:1257137/1‑145 (MQ=255)
                                                                  aTGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCCTGGCCCCGCACGGGCCGGGCCGGGGGGTCTCAAATGCGGGGCCCGCCGCGTAGCGGGCCTGCCCGCGCCcccacggtt                                                                                    >  2:2857876/1‑145 (MQ=255)
                                                                         acaAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGGCTGCCCCCGAACCGGCCAGGGCGGCGGGTCTCTATTGCGGGGCCGGCCGCATCGGGGGGCTGCCCCCGCCCACACGTTTAGTGGtt                                                                             >  2:49530/1‑150 (MQ=255)
                                                                           aaGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCAGCTGATCTCTATGGCTGAGCCGGCCACGTCGCGGGCAAGCCCCCTCCCACAGGTTTCGTGGTTgg                                                                           >  2:3564014/1‑150 (MQ=255)
                                                                              ctgctctgatTGCGACATCCGGGGGCCCGAACCTGGCGGACGCCTTT‑‑GTGGGCGCGGGCATCCCCGCGACCGGGCCGGGCCTGCTGATCTCTATGGCGGCGCCGGCCTCTTCGGGGGCAGGCCCGCTCCCACAGGTTTCGTGGTTGGCTc                                                                        <  2:417104/150‑1 (MQ=255)
                                                                                       ttGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCGTCGCGGGCATGCCCGCACCCACACGTTTCGTGGTTGGCTCGCAACGGCa                                                               >  2:2827002/1‑150 (MQ=255)
                                                                                           gACATACGGGTGCCCGAAACGTGCCGAAGCCTTT‑‑GTGGGCGCGGGCATGCCCGCGAACGGGCCGGGCCTGCTGATCTCTTTGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTcc                                                           <  2:4973659/150‑1 (MQ=39)
                                                                                                   ggTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCCAACGGGCCTGGCCGGCTGCTCTCTATGGCTGAGCCGGCCGCGTCGCGGGGCTGCCCGCTCCCACACGTTTAGGGGGTGGCTCGCAAGGGCAGTCCCATCCCGc                                                   >  2:979384/1‑150 (MQ=255)
                                                                                                   ggTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCCCGAACCGGCCTGGCCGGCTGATCTCTATGGCTGAGCCGGCCTCGTCGCGGGCATGCCCGCGCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCAATCCAAAACCGc                                                   >  2:6702706/1‑150 (MQ=255)
                                                                                                     ccttgcgagccaaccacgaagccTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgg                                                                                                          <  1:3305385‑M1/72‑1 (MQ=255)
                                                                                                        ccGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCAGGCCCGCGAACGGGCCTGGCCTGCTGCTCTCGATGGCTGTGCCGGCCGCTTCGCGGGCCTGCCCGCGCCCCCACGTTTCGTGGTTGGGTCGCAAGGGCCTTCCAAAACCGCCAAcc                                              >  2:2632415/1‑150 (MQ=255)
                                                                                                          taacgttGCCGTCGCCTTT‑‑GTGGGCGCGGGCTGGCCCGCTACGGGGCCTGGCCTGCTGAGCTCTATGGCTGCGCCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCtt                                            <  1:6213631/145‑1 (MQ=37)
                                                                                                          gAAACTGGCCGAAGCCTTG‑‑GTGGGCGCGGGCCTGCCCGCGACCGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCACTTCGCGGGCAGGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCtt                                            <  2:1495011/150‑1 (MQ=255)
                                                                                                           caaCTTGCCGACGCCTTT‑‑GGGGGCGGGGGCAGGCCCGCCACGGGGCGGGGCCTGCTGATCTCTTGGGGTGCGCCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTg                                           <  2:1801454/149‑1 (MQ=25)
                                                                                                           aacctggCCGAAGCCTCG‑‑GTGGGAGCGGGCAGGCCCGGTACCGGGCCTGGCCTGCTGATCTATATGGCTGAGCCGGCCTCTGCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTg                                           <  2:1801391/144‑1 (MQ=38)
                                                                                                            accTTGCCGAAGCCTAT‑‑GTGGGTGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGt                                          <  2:3398654/148‑1 (MQ=255)
                                                                                                            aaCTTGCCGACGCCTTT‑‑GTGGGAGGGGGCAGGCCCGCGACCGGGCGGGGCGTGCTTATCTCTTGGGCTGCGCGGGCCTGTTCGGGGGCTTGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGt                                          <  2:3032029/150‑1 (MQ=255)
                                                                                                             aCTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCGGGCCGGCGGATCTCGATGGCGGAGCCCGCCACGTCGCCGGCAGGCCCGCTCCCACACGTTTCGTGGTTGGCTCGCAAGGGCAATCCAAACCCCCCAACCATgga                                         >  2:3375155/1‑148 (MQ=255)
                                                                                                              cTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCGATGGCGGAGCCGGCCTCGTCGCGGGCATGCCCGCTCCCCCAGGTGTCGGGGGTGGCGCGCAAGGGCCTTCCCCGCCCGCCATCCTTGTAg                                        >  2:5581800/1‑150 (MQ=255)
                                                                                                               tggCCGACGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGACCGGGCCTGGCCTGCTGATCTCTCTGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAgg                                       <  2:3248414/148‑1 (MQ=255)
                                                                                                                      cgCCTTT‑‑GTGGGCGGGGGCATGCCCGCGACGGGGCGTGGCCTGCTTATCTCTTTGGGTGCGCCGGCCTGTTCGGGGGCTTGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATcc                                <  2:1401404/149‑1 (MQ=255)
                                                                                                                         cTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTa                             <  2:6224685/150‑1 (MQ=255)
                                                                                                                           cT‑‑GTGGGAGCGGGCATGCCCGTGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTaca                           <  1:3367666‑M1/149‑1 (MQ=255)
                                                                                                                            t‑‑gtgGGGGCGGGGCCTGCCCCGCGAGCGGGCCGGCCGTGCTGATCTCTTTGGCTGCCCCGCCCTCGGCGGGGCCAGCCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTac                            <  1:4244090/150‑1 (MQ=255)
                                                                                                                               gtgGGAGCGGGCATGCCCGCCAACGGGCCTGGCCTGCGGTGCTCTATTGCAGAGCCCGCCTGCTGGCGGGCATGCCCGCTCCCACAGGTTTGGTGGTTTGGTCGCAAGGGCCATCCCATACCGCCAACCGTGTAGGGCGAACCCAacacc                         >  2:4807353/1‑148 (MQ=255)
                                                                                                                               ccgttcGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCGATGGCTGCGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCCTTCCCATCCCGCCAACCTTGGAGGGCGCTCCCTacaac                         <  2:5284317‑M1/144‑1 (MQ=255)
                                                                                                                                        GGGCAGGCCGGCTATGGGGCCGGGCGTGGTTATTTCTTTGGCTGAGCGGGCCTCTTCGGGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTc                 >  2:1214517‑M1/1‑150 (MQ=255)
                                                                                                                                         ggatTGCCCGCGACCGGGCCTGGCCTGCTGATCTCTTTGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTcc                >  2:1552405‑M1/5‑150 (MQ=255)
                                                                                                                                         GGCTTGCCCGCGAAGGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTcc                >  2:5333693‑M1/1‑150 (MQ=255)
                                                                                                                                                CCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCgg          <  2:3459606‑M1/150‑1 (MQ=255)
                                                                                                                                                CCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCgg          <  1:2575206‑M1/150‑1 (MQ=255)
                                                                                                                                                 CCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCCTTCCAATCCCGCCATCCGTGTAGGGCGATCCCTACAACAAATCTTCCATCCGga         <  1:3618862‑M1/150‑1 (MQ=255)
                                                                                                                                                 CCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGga         <  1:6198212‑M1/150‑1 (MQ=255)
                                                                                                                                                 CCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGga         <  2:5886630‑M1/150‑1 (MQ=255)
                                                                                                                                                 CCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGga         <  2:4547482‑M1/150‑1 (MQ=255)
                                                                                                                                                 CCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGga         <  2:3858498‑M1/150‑1 (MQ=255)
                                                                                                                                                  cgtgAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGgat        <  2:3367666/147‑1 (MQ=255)
                                                                                                                                                  CGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCCTTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCAGCCGgat        <  1:1877363‑M1/150‑1 (MQ=255)
                                                                                                                                                  CGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATACGgat        <  1:5912450‑M1/150‑1 (MQ=255)
                                                                                                                                                  CGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTAGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGAGCCCTACAACAAATATTCCAACCGgat        <  1:2285161‑M1/150‑1 (MQ=255)
                                                                                                                                                      AACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGatgatg    <  1:4588491‑M1/150‑1 (MQ=255)
                                                                                                                                                       ACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATgg   <  1:4970194‑M1/150‑1 (MQ=255)
                                                                                                                                                        cGGGCCTGGCGTTTTTTTCTATTTGGTTGAGCGGGCCTCTGCGGGGTCTTGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGc  <  2:4590266/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
ACCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTAT‑‑GTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGC  >  NC_002947/499073‑499370

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: