Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,070,248 +GA intergenic (+48/+94) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2461.G95.0% 277.1 / 8.5 100intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base G (48/47);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01
*NC_0029471,070,2462.A95.0% 280.4 / 9.5 100intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base A (48/47);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 101 total aligned reads displayed.

GGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAA  >  NC_002947/1070122‑1070360
                                                                                                                             ||                                                                                                                  
ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                      <  1:1013482/125‑1 (MQ=255)
    aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gaga                                                                                                                <  1:1013485/125‑3 (MQ=255)
    aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGATGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gaga                                                                                                                <  1:1013484/125‑3 (MQ=255)
      gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatc                                                                                                              <  1:1013488/125‑5 (MQ=255)
      gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatc                                                                                                              >  1:1013630/1‑121 (MQ=255)
        tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑Gagatcg                                                                                                             >  1:1013691/1‑119 (MQ=255)
             ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGcgcc                                                                                                         >  1:1013904/1‑125 (MQ=255)
              cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAgcgc                                                                                                          >  1:1013684/1‑123 (MQ=255)
               ggCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGcgccgc                                                                                                       >  1:1013606/1‑125 (MQ=255)
               ggCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGcgccg                                                                                                        >  1:1013829/1‑124 (MQ=255)
                  aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCcgcgcg                                                                                                    <  1:1013497/125‑1 (MQ=255)
                  aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCcgcgc                                                                                                     <  1:1013500/124‑1 (MQ=255)
                    ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                  <  1:1013503/125‑1 (MQ=255)
                    ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGG‑TTTCGGCGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                  <  1:1013504/124‑1 (MQ=39)
                     gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggc                                                                                                 >  1:1013626/1‑125 (MQ=255)
                     gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                  >  1:1013695/1‑124 (MQ=255)
                     gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                  >  1:1013673/1‑124 (MQ=255)
                       tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                               <  1:1013505/125‑1 (MQ=255)
                       tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                               <  1:1013506/125‑1 (MQ=255)
                        cGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                               <  1:1013507/124‑1 (MQ=255)
                         gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                             <  1:1013510/125‑1 (MQ=255)
                           tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                            >  1:1013736/1‑124 (MQ=255)
                           tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                           <  1:1013512/125‑1 (MQ=255)
                             tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCg                                                                                         >  1:1013747/1‑125 (MQ=255)
                                 gtgggtggCTTTCAAGCTGAAGAAGGCTGGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGTCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGTGGCGCTCGAtc                                                                                      <  1:1013516/124‑1 (MQ=38)
                                 gtgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                     >  1:1013777/1‑125 (MQ=255)
                                       ggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                               <  1:1013520/125‑1 (MQ=255)
                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                             >  1:1013719/1‑124 (MQ=255)
                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                             >  1:1013714/1‑124 (MQ=255)
                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                            <  1:1013521/125‑1 (MQ=255)
                                           ttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                             >  1:1013866/1‑123 (MQ=255)
                                           ttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGCTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                            >  1:1013690/1‑124 (MQ=37)
                                             cAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGa                                                                         <  1:1013523/125‑1 (MQ=255)
                                              aaGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaa                                                                        >  1:1013858/1‑125 (MQ=255)
                                               aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                       >  1:1013770/1‑125 (MQ=255)
                                               aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                       <  1:1013524/125‑1 (MQ=255)
                                               aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaa                                                                        <  1:1013525/124‑1 (MQ=255)
                                                  tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                    >  1:1013636/1‑125 (MQ=255)
                                                     agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                 >  1:1013825/1‑125 (MQ=255)
                                                      gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                <  1:1013530/125‑1 (MQ=255)
                                                       aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                  >  1:1013868/1‑122 (MQ=255)
                                                         ggCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                              >  1:1013707/1‑124 (MQ=255)
                                                         ggCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGc                                                             <  1:1013531/125‑1 (MQ=255)
                                                          gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                            >  1:1013887/1‑125 (MQ=255)
                                                             cGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                          <  1:1013537/124‑1 (MQ=255)
                                                             cGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                          <  1:1013536/124‑1 (MQ=255)
                                                              gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                           <  1:1013538/122‑1 (MQ=255)
                                                                cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                      <  1:1013542/125‑1 (MQ=255)
                                                                 aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                      >  1:1013688/1‑124 (MQ=255)
                                                                 aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                     <  1:1013544/125‑1 (MQ=255)
                                                                 aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                     >  1:1013692/1‑125 (MQ=255)
                                                                  ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                    <  1:1013545/125‑1 (MQ=255)
                                                                       aGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                               <  1:1013548/125‑1 (MQ=255)
                                                                       aGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                               <  1:1013549/125‑1 (MQ=255)
                                                                        gggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                              >  1:1013853/1‑125 (MQ=255)
                                                                         ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                   <  1:1013552/87‑1 (MQ=25)
                                                                         ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGccc                                             <  1:1013553/125‑1 (MQ=255)
                                                                         ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGccc                                             <  1:1013554/125‑1 (MQ=255)
                                                                          gCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                              >  1:1013903/1‑123 (MQ=255)
                                                                             gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGa                                          <  1:1013559/124‑1 (MQ=255)
                                                                              aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                        >  1:1013694/1‑125 (MQ=255)
                                                                                  ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                    <  1:1013564/125‑1 (MQ=255)
                                                                                   gTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGCCCCTGACGCGTTc                                   >  1:1013872/1‑125 (MQ=37)
                                                                                    tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                  <  1:1013565/125‑1 (MQ=255)
                                                                                    tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                  >  1:1013959/1‑125 (MQ=255)
                                                                                      tGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                                <  1:1013570/125‑1 (MQ=255)
                                                                                       gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCAACTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                               <  1:1013571/125‑1 (MQ=37)
                                                                                        ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                              <  1:1013572/125‑1 (MQ=255)
                                                                                         tCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                             <  1:1013574/125‑1 (MQ=255)
                                                                                             aaagggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                              <  1:1013575/117‑1 (MQ=37)
                                                                                             aaGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                          <  1:1013577/124‑1 (MQ=255)
                                                                                              aGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                        >  1:1013801/1‑125 (MQ=255)
                                                                                              aGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                        >  1:1013857/1‑125 (MQ=255)
                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                               >  1:1013662/1‑101 (MQ=255)
                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                               <  1:1013579/117‑1 (MQ=255)
                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                            <  1:1013578/120‑1 (MQ=255)
                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                       >  1:1013669/1‑125 (MQ=255)
                                                                                                gggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                      >  1:1013787/1‑125 (MQ=255)
                                                                                                gggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCAATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                      >  1:1013908/1‑125 (MQ=38)
                                                                                                 ggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                     >  1:1013811/1‑125 (MQ=255)
                                                                                                 ggTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcc                       >  1:1013826/1‑123 (MQ=255)
                                                                                                 ggTATGCCGTTAGGCTTGCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                     <  1:1013582/125‑1 (MQ=37)
                                                                                                  gTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                     >  1:1013812/1‑124 (MQ=255)
                                                                                                  gTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                     >  1:1013788/1‑124 (MQ=255)
                                                                                                  gTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                     >  1:1013961/1‑124 (MQ=255)
                                                                                                  gTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                    <  1:1013583/125‑1 (MQ=255)
                                                                                                   tATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTc                   >  1:1013748/1‑125 (MQ=255)
                                                                                                        cGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGTGGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc              <  1:1013586/125‑1 (MQ=37)
                                                                                                          ttAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctg            <  1:1013587/125‑1 (MQ=255)
                                                                                                          ttAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct             >  1:1013852/1‑124 (MQ=255)
                                                                                                           tAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc           >  1:1014041/1‑125 (MQ=255)
                                                                                                            atgTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc          >  1:1013873/3‑125 (MQ=255)
                                                                                                             ggTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCt         <  1:1013589/125‑1 (MQ=255)
                                                                                                              gTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc           <  1:1013591/122‑1 (MQ=255)
                                                                                                              gTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc        <  1:1013592/125‑1 (MQ=255)
                                                                                                              gTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc        <  1:1013590/125‑1 (MQ=255)
                                                                                                               ttttCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc       >  1:1013717/1‑125 (MQ=255)
                                                                                                               ttttCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc       >  1:1014021/1‑125 (MQ=255)
                                                                                                               ttttCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc       >  1:1014036/1‑125 (MQ=255)
                                                                                                                    ggCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa  >  1:1013849/1‑125 (MQ=255)
                                                                                                                             ||                                                                                                                  
GGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAA  >  NC_002947/1070122‑1070360

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑