Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,508,456 T→G 10.7% M438R (ATG→AGG)  PP_4839 → iron‑regulated membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,508,4560TG10.7% 304.0 / 2.6 141M438R (ATG→AGG) PP_4839iron‑regulated membrane protein
Reads supporting (aligned to +/- strand):  ref base T (53/64);  new base G (5/13);  total (58/83)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

GGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTA  >  NC_002947/5508308‑5508605
                                                                                                                                                    |                                                                                                                                                     
gggCTGGTGATGTGGTGGATGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAtg                                                                                                                                                      <  2:1612689/150‑1 (MQ=255)
gggCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAtg                                                                                                                                                      <  2:3888109/150‑1 (MQ=255)
   cTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTgg                                                                                                                                                   <  2:1542091/150‑1 (MQ=255)
    tGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTggg                                                                                                                                                  >  1:4141735/1‑150 (MQ=255)
    tGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCGGGTGGGGGTGTCGCTGGTGGTGATGTggg                                                                                                                                                  >  2:1364225/1‑150 (MQ=255)
    tGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCGGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTggg                                                                                                                                                  <  2:962241/150‑1 (MQ=255)
        gATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGTTGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGt                                                                                                                                              <  1:1045644/150‑1 (MQ=255)
        gATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGt                                                                                                                                              <  2:1488137/150‑1 (MQ=255)
         aTGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTg                                                                                                                                             <  2:1593528/150‑1 (MQ=255)
         aTGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTg                                                                                                                                             <  2:1004567/150‑1 (MQ=255)
           gtggtgGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGTCCTGCCACGCTGGATGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGa                                                                                                                                           <  1:285371/150‑1 (MQ=255)
           gtggtgGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGGGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGa                                                                                                                                           <  2:666143/150‑1 (MQ=255)
            tggtggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGAt                                                                                                                                          <  1:2139028/150‑1 (MQ=255)
             ggtggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATa                                                                                                                                         <  1:1027825/150‑1 (MQ=255)
             ggtggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCGGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATa                                                                                                                                         <  1:1956085/150‑1 (MQ=255)
              gtggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAg                                                                                                                                        >  1:3196642/1‑150 (MQ=255)
              gtggAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCGGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGGGGTGAGGGGGGGGGGGGATAg                                                                                                                                        >  2:2746202/1‑150 (MQ=255)
                  aaGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtg                                                                                                                                    <  2:2047433/150‑1 (MQ=255)
                  aaGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtg                                                                                                                                    >  1:2508600/1‑150 (MQ=255)
                  aaGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtg                                                                                                                                    >  1:4518886/1‑150 (MQ=255)
                  aaGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtg                                                                                                                                    <  1:2764073/150‑1 (MQ=255)
                  aaGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtg                                                                                                                                    <  2:3196642/150‑1 (MQ=255)
                   aGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtgg                                                                                                                                   <  1:674596/150‑1 (MQ=255)
                   aGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtgg                                                                                                                                   <  2:765523/150‑1 (MQ=255)
                    gcgcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggt                                                                                                                                  <  2:1011832/150‑1 (MQ=255)
                     cgcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggtg                                                                                                                                 <  1:4421581/150‑1 (MQ=255)
                      gcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggtg                                                                                                                                 <  1:3801929/149‑1 (MQ=255)
                      gcCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTtggtg                                                                                                                                 >  2:3801929/1‑149 (MQ=255)
                        cGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGgtggtg                                                                                                                              <  1:1161940/150‑1 (MQ=255)
                        cGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGgtggtg                                                                                                                              >  1:2051372/1‑150 (MQ=255)
                         gTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCGGGGGGTGATGGTCCCGCGGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGa                                                                                                                             <  1:1207888/150‑1 (MQ=255)
                          tCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGAGAGTTTGGTGGTGAg                                                                                                                            >  1:748014/1‑150 (MQ=255)
                            cTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCCGGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGGTCCCGCTGGTGGGGGGGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGgc                                                                                                                          <  1:1019160/150‑1 (MQ=255)
                            cTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCCGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCGGGTGGGGGTGTCGATGGTGGTGCGGTGGGGGGGGGATAGTTTGGTGGTGAGgc                                                                                                                          >  2:4000331/1‑150 (MQ=255)
                            cTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGAGAGTGTGGTGGTGAGgc                                                                                                                          >  2:4000315/1‑150 (MQ=255)
                            cTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGgc                                                                                                                          >  2:306980/1‑150 (MQ=255)
                               acggcggcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGGGGCCGTGATGCGGGCGCTGGGGGTGATGGTCCCGCTGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGcgcc                                                                                                                       <  1:2320520/148‑1 (MQ=255)
                               aaGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGcgcc                                                                                                                       >  1:2492092/1‑150 (MQ=255)
                                    ggcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCgg                                                                                                                  <  2:2591267/150‑1 (MQ=255)
                                    ggcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCgg                                                                                                                  <  2:2671727/150‑1 (MQ=255)
                                    ggcCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCgg                                                                                                                  <  2:872112/150‑1 (MQ=255)
                                    gcccTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTTATGTGGGGGGTGGATAGTTTGGTGGTGAGGCgcctcctg                                                                                                                  <  2:4505208/148‑3 (MQ=255)
                                       cTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTTATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTc                                                                                                               >  1:1552748/1‑150 (MQ=255)
                                       cTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTc                                                                                                               >  1:3828835/1‑150 (MQ=255)
                                       cTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGCGGGTGTCGATGGTGGTGATGTGGGGGGTGGAGAGTGTGGGGGTGAGGCGCCGCCGGGTc                                                                                                               >  2:1964658/1‑150 (MQ=255)
                                       cTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCGGGTGCGGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTc                                                                                                               >  2:4465739/1‑150 (MQ=255)
                                          ggCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTg                                                                                                            <  2:2454415/150‑1 (MQ=255)
                                          ggCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTg                                                                                                            <  2:2454360/150‑1 (MQ=255)
                                          ggCGGACCGCCTGTGCGTCACGCCCCGCCACGCTGGAAGGCGGCAGGGGCCGTGATGCGGGTGCGGGGGGTGAGGGCCCCGCGGGTGGGGGTGCCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTg                                                                                                            <  1:3040034/150‑1 (MQ=255)
                                                  gCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAgcgc                                                                                                    >  2:3460150/1‑150 (MQ=255)
                                                   ccTTCGCGCCATGACCTGCCACGCTGGACGGCGGCAGTGGCCGTGATGCTGGCGCTGGGGGTGATGGCCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCt                                                                                                   <  2:3435497/150‑1 (MQ=255)
                                                    cTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCtt                                                                                                  <  2:1731038/150‑1 (MQ=255)
                                                    cTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGGGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCtt                                                                                                  <  2:4094366/150‑1 (MQ=255)
                                                     tttGCGTCATGACCTGCCACGCTGGTAGGCGGCAGTGGCCGTGAGGCGGGTGCTGGGGGTGCTGTTCCCGCGGGTGGGGGTGTCGCTGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAgcacagg                                                                                                 <  2:4018347/150‑6 (MQ=255)
                                                           tCATGACCCGCCCCGCGGGAAGGCGGCAGTGGCCGCGCGGCGGGTGCTGGGGGGGATGGTCCCGCTGGTGGGGGTGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCtt                                                                                           <  2:3300303/150‑1 (MQ=255)
                                                            cATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTtg                                                                                          <  2:2492092/150‑1 (MQ=255)
                                                             cggccccgCCACGCGGGCCGGCGGCCGCGGCCGTGCGGCCGGTGCGGGGGGCGACGGGCCCGCGGGCGGGGGTGCCGCTGGCGGGGCTGCGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTtgt                                                                                         <  1:1735350/143‑1 (MQ=255)
                                                             aTGACCTGCCACGCTGGAAGCCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGACGTGGGGGGTGGATAGTATGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTtgt                                                                                         <  2:1169827/150‑1 (MQ=255)
                                                                  ccgCCCCGCTGGAAGGCGGCAGGGGCCGCGAGGCTGGTGCCGGGGGTGCTGTTCCCGCTGGGGGGGGGGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGc                                                                                    <  2:2068476/148‑1 (MQ=255)
                                                                  cTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGc                                                                                    >  1:602651/1‑150 (MQ=255)
                                                                  cTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGc                                                                                    >  1:1603430/1‑150 (MQ=255)
                                                                  cTGCCACGCTGCCCGGCGGCAGGGGCCGTGCCGCTGGTGCCGGGGGTGCTGGGCCCGCTGGTGGGGGCGGCGATGGCGGGGAGGCGGGGGGCGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGc                                                                                    <  1:1378126/150‑1 (MQ=255)
                                                                   tGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGcc                                                                                   >  1:582781/1‑150 (MQ=255)
                                                                    gCCACGCGGGCCGGCGGCCGGGGCCGTGCGGCGGGGGCGGGGGGGGAGGGTCCCGCGGGCGGGGGGGGCGATGGGGGTGCTGTGGGGGGTGGATAGTTGGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCt                                                                                  <  2:2929593/150‑1 (MQ=255)
                                                                     ccACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCtg                                                                                 >  1:2008032/1‑150 (MQ=255)
                                                                       aCGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCtgtg                                                                               >  2:2353324/1‑150 (MQ=255)
                                                                       aCGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCtgtg                                                                               >  1:804309/1‑150 (MQ=255)
                                                                        cGCTGGAAGGCGGCCGGGGCCGTGATGCGGGTGCGGGGGGGGAGGGGCCCGCGGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGc                                                                              <  2:2541405/150‑1 (MQ=255)
                                                                            ggccggcggcAGTGGCCGCGCTGCGGGTGCGGGGGGTGCCGGGCCCGCGGGGGGGGGGGGCGCTGGGGGTGAGGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGc                                                                          <  2:2937543/146‑1 (MQ=255)
                                                                            ggACGGCGGCAGTGGCCGTGCTGCGGGTGCGGGGGGTGCGGGGCCCGCGGGTGGGGGTGTCGAGGGTGGTGCGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGc                                                                          <  2:668522/150‑1 (MQ=255)
                                                                              ccggcggcAGTGGCCGTGCTGCTGGCGCGGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGCTGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCct                                                                        <  2:2517675/148‑1 (MQ=255)
                                                                                 gcggcAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCtt                                                                     >  1:3954584/1‑150 (MQ=255)
                                                                                     ctgTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTcgcg                                                                 <  2:3748868/148‑1 (MQ=255)
                                                                                     cAGCGGCCGAGTGGCTGGGGCGGCGGGTCCTGCGCCCGCGGCCGGGGGTGTCGACGGTGGCGATGCGGGCGGTGGATAGTTGGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTcgcg                                                                 <  1:3382374/150‑1 (MQ=255)
                                                                                      cggggCCGTGATGCGGGCGCCGGGGGCCCCGCCCCCGCGGGCGGGGGTGTCGAGGGTGGTGCGGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCgg                                                                <  1:3867679/147‑1 (MQ=255)
                                                                                      agcggCCGTGATGCTGGCGCCGGGGGTGCCGGGCCCGCTGGGGGGGGTGCCGATGGTGGTGCCGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCgg                                                                <  2:317626/147‑1 (MQ=255)
                                                                                      aGTGGCCGTGTTGCTGGTGCTGGGGTTGATGTTCCCGCTGGTGGGTGTGTCGATGTTGGTGATGGGTGGGGTGGATATTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCgg                                                                <  2:1213665/150‑1 (MQ=255)
                                                                                        tGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGt                                                              >  1:641501/1‑150 (MQ=255)
                                                                                            cGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                          >  1:3123542/1‑150 (MQ=255)
                                                                                            cGTGACGCTGGTGCTGGGGGTGCTGTTCCCGCCGGTGGGGGTGTCGACCGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                          <  2:2135989/150‑1 (MQ=255)
                                                                                              tgctgctggtgctggGGGCGACGGACCCGCTGGCGGGGGCGTCGATGGTGGCGAGGAGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAccc                                                        <  2:690860/147‑1 (MQ=255)
                                                                                               gATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCg                                                       >  2:3176437/1‑150 (MQ=255)
                                                                                                aTGCTGGGGCGGGGGGGGACGCGCCCGCGGGGCGGGGCGTCGAGGGCGGTGAGGTGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGc                                                      <  1:3931380/150‑1 (MQ=37)
                                                                                                    tggtgctggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGGGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCt                                                  >  2:4115112/1‑150 (MQ=255)
                                                                                                       tgctggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACg                                               >  2:2591098/1‑150 (MQ=255)
                                                                                                        gctggGGGTGCTGCGCCCGCGGGTGGGGGCGCCGCTGGCGGCGCAGTGGGGGGTGGATCGTTGGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGa                                              <  1:2739964/150‑1 (MQ=25)
                                                                                                          tggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCAACGAgg                                            >  1:2881441/1‑150 (MQ=255)
                                                                                                          tggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGGAAAACCCCGCCAAAGAgg                                            >  1:148421/1‑150 (MQ=255)
                                                                                                          ggggggTGCTGGTCCCGCGGGGGGGGGTGGCGCGGGGGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAgg                                            <  2:3650751/149‑1 (MQ=37)
                                                                                                           gggggTGATCGTCCCGCGGGTGGGGGCGGCGCCGGTGGTGATGTGGGGGGTGGAGAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAggg                                           <  2:3412584/150‑1 (MQ=38)
                                                                                                           gggggCGAGGTTCCCGCTGGTGGGGGTGTCGAGGGCGGTGAGGTGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAggg                                           <  2:3412566/150‑1 (MQ=38)
                                                                                                               gTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACcg                                       >  2:1634326/1‑150 (MQ=255)
                                                                                                                 gATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGAGGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCGCCAACCAGGGCCcccg                                     >  1:902593/1‑150 (MQ=255)
                                                                                                                   cccccccGCGGGCGGGGGCGCCGACGGCGGCGACGCGGCGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGtt                                   <  1:4065219/146‑1 (MQ=37)
                                                                                                                      tCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGGAAAACCCCACCAACCAGGGACCCCGTTGTa                                >  1:1930428/1‑150 (MQ=255)
                                                                                                                      gcccGCGGGGGGGGGGGTCGCTGGGGGGGCGGGGGGGGGTGGAGAGTTGGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCGAGGGACCGCGTTGTa                                <  2:2902817/149‑1 (MQ=37)
                                                                                                                       cccGCGGGTGGGGGTGTCGATGGTGGTGCTGGGGGGGGGGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAg                               <  2:3119574/150‑1 (MQ=255)
                                                                                                                              gtggGGGTCTCGATGGTGGTGATGTGGGGGGCGGAGAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCggg                        <  2:2367876/150‑1 (MQ=255)
                                                                                                                                         gATGGTGGTGATGTGGGGGGTGGCTAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCAAAGAGGGACCGCGGTGTAGGAGCGGGTGTGCCcgcgc             >  1:3143906/1‑149 (MQ=38)
                                                                                                                                                    tgtgGGGGCTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTCGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGCGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGccagca  <  2:4518886/150‑5 (MQ=14)
                                                                                                                                                    |                                                                                                                                                     
GGGCTGGTGATGTGGTGGAAGCGCCGTCCTGAAGGCGGCCTGGGCGTACCGCCTTTGCGTCATGACCTGCCACGCTGGAAGGCGGCAGTGGCCGTGATGCTGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTA  >  NC_002947/5508308‑5508605

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: