Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,878,998 G→C 7.5% intergenic (+42/+63) rpiA → / ← PP_5151 ribose‑5‑phosphate isomerase A/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,878,9980GC 7.5% 691.8 / 7.0 216intergenic (+42/+63)rpiA/PP_5151ribose‑5‑phosphate isomerase A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (67/128);  new base C (16/4);  total (83/133)
Fisher's exact test for biased strand distribution p-value = 1.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

ATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATCCCC  >  NC_002947/5878849‑5879145
                                                                                                                                                     |                                                                                                                                                   
aTCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCAAAAAGGTACCCCACCg                                                                                                                                                     >  1:1859462/1‑150 (MQ=255)
aTCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCAGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCACCTACAACGTACCCCCCCg                                                                                                                                                     >  1:3694278/1‑150 (MQ=255)
    aCGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCt                                                                                                                                                 <  2:1189923/150‑1 (MQ=255)
      gCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTgg                                                                                                                                               >  2:3241444/1‑150 (MQ=255)
        cATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGaa                                                                                                                                             <  2:4122806/150‑1 (MQ=255)
        cATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGaa                                                                                                                                             <  2:2369706/150‑1 (MQ=255)
          tcgtcgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAAcc                                                                                                                                           <  2:2230443/150‑1 (MQ=255)
              cgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                                       <  1:2274596/150‑1 (MQ=255)
                 cGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGa                                                                                                                                    <  1:988448/150‑1 (MQ=255)
                    ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTTAAACCCCGCCAAAAAAGTAACCCCCCCGCCCGGAACCCGGGCGCTcc                                                                                                                                 >  1:1961429/1‑150 (MQ=255)
                      tGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCt                                                                                                                               >  1:203131‑M2/1‑149 (MQ=255)
                       gACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtg                                                                                                                              >  2:399411‑M2/1‑148 (MQ=255)
                       gACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCAACGACAAGGTAACCCACCCGCCGGGTACCGTGGCGATaccag                                                                                                                              >  1:2557291/1‑148 (MQ=255)
                         ccAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgta                                                                                                                            >  1:2101058‑M2/1‑146 (MQ=255)
                         atgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTa                                                                                                                            >  2:46093‑M1/89‑150 (MQ=255)
                          tgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAg                                                                                                                           >  1:4164768‑M1/88‑150 (MQ=255)
                          cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCTTGTAAAACCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGTCGATCCCtgtag                                                                                                                           >  2:2955466‑M2/1‑145 (MQ=255)
                          cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCAACAAGGTACCCCACCGGCCTGGGACCGGGGCGATcccaggag                                                                                                                           >  2:1728327/1‑145 (MQ=255)
                           aaCGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCAACAAGGTACCCCACCGGCCTGGAACCCTGGCGATCCCggtagg                                                                                                                          <  1:3344372‑M2/150‑7 (MQ=255)
                            aCGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACACGCTCCAACAAAGTACCCCACCCGCCAGGAACCGTGGCGATACCagtagga                                                                                                                         <  1:2139878‑M2/150‑8 (MQ=255)
                              ggCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagc                                                                                                                       >  1:857401‑M2/1‑141 (MQ=255)
                                      cgcggcgcggCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttac                                                                                                               >  2:169410‑M2/1‑133 (MQ=255)
                                      cgcggcgcggCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttac                                                                                                               >  2:2031507‑M2/1‑133 (MQ=255)
                                      catttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTa                                                                                                                <  2:1303338‑M1/74‑1 (MQ=255)
                                      catttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTa                                                                                                                >  1:1303338‑M1/76‑149 (MQ=255)
                                        cggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttaccc                                                                                                             >  1:2961733‑M2/1‑131 (MQ=255)
                                        cggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttaccc                                                                                                             >  2:2107881‑M2/1‑131 (MQ=255)
                                          gcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgc                                                                                                           >  1:724292‑M2/1‑129 (MQ=255)
                                             cggcgggcGGATCTGCTGGTGCTGGGCACCGCTGACGGCGTGAAGAGCGTGATGGCCGAGTAACCCCCTTCGCGGGTAACCCCGCTCCTCCACGGTACCCCACCGGCCTTGAAGCGTGGCGATCCCtgtaggagcgggtttacccgcgaa                                                                                                        >  2:171814‑M2/1‑126 (MQ=255)
                                               gccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACa                                                                                                      <  1:1815532/150‑1 (MQ=255)
                                                ccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCAACAAGGTACCCCACCGGCCGGGAACCGTGGCGCTCCCTGGAGGAGCGGGTTTACCCGCGAACAg                                                                                                     >  2:1411256/1‑150 (MQ=255)
                                                  ggcggATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCGGGAACCGTGGCGCTCCCTGTAGGAGCGGGTTTACCCGCGAACAGGc                                                                                                   >  2:3424753/1‑150 (MQ=255)
                                                        tCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCCCTCCTACAAGGTACCCCACCCGCCTGGAACCGTGGCGATCCCAGTAGGAGCGGGTTTACCCGCGACCAGGCCCAGAc                                                                                             >  2:1369824/1‑150 (MQ=255)
                                                        tCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCTCCGACAAGGTACCCCACCGGCCGGGAACCGTGGCGATCCCAGTAGGCGCGGGGGTTCCCGCGAACAGGCCCAGAc                                                                                             >  1:2917387/1‑150 (MQ=255)
                                                        tCTGCTGCTGCTGGGCACCGCTAAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGAc                                                                                             <  2:622286/150‑1 (MQ=255)
                                                           gctggtgctgGGCACCGCTGATGGCGTTAAGCGCCTTAAGGCCGTGTAAAACCCTGCGCGGTTAAACCGGCTCCTACACGGTACCCCACGGGCCTGTAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAgg                                                                                          <  1:3557114/150‑1 (MQ=255)
                                                             tgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCg                                                                                        <  1:4032166/150‑1 (MQ=255)
                                                             tgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGGACCCCACCGGCCGGGAACCGTGGCGATCCCGGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCg                                                                                        >  2:4448122/1‑150 (MQ=255)
                                                              gctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTTAACCCCCTCCGACAAGGTACCCCACCGGCCGGGGACCGGGGCGGTCCCAGTAGGAGCGGGGTTACCCGCGAACAGGCCCGGACAGGCgc                                                                                       >  2:3031984/1‑149 (MQ=255)
                                                              gctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCTCCAACAAGGTACCCCACCCGCCTGGAACCCTGGCGATACCAGTAGGAGCGGGTTTAACCCCCAACAGGCCCAGACAGGCGa                                                                                       >  2:154092/1‑150 (MQ=255)
                                                              gctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCGCCAACCAGGTACCCCACCGGCCGGGAAACGAGGCGCTCCCGGTAGGAGCGGGGTTACCCCCGCACACGCCCTGACAGGCgc                                                                                       >  2:1884797/1‑149 (MQ=255)
                                                                 gctgGGCACCGCTGATGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                                    <  1:666626/150‑1 (MQ=255)
                                                                 gctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                                    <  2:541026/150‑1 (MQ=255)
                                                                 gctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                                    <  2:1727326/150‑1 (MQ=255)
                                                                 gctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                                    <  2:86566/150‑1 (MQ=255)
                                                                   tgGGCACCGCTGACGGCGTGACGAGCCTGACGGCCGAGTAAAACCCTGCGCGGTTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTTTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagt                                                                                  <  1:1581876/150‑1 (MQ=255)
                                                                   tgGGCACCGCTGAAGGCGTGACGAGCCTGAAGGCCGAGTAAACCCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagt                                                                                  <  1:1364028/150‑1 (MQ=255)
                                                                   tgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCCCGCCAAAAAGGTACCCCACCCGCCGGGAACCCGGGCGATCCCAGGAGGTGCGGGTTTAACCGCCAAAAGGCCCAGACAGGCGCtcagt                                                                                  >  1:1094207/1‑150 (MQ=255)
                                                                    gggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtc                                                                                 <  1:343011/150‑1 (MQ=255)
                                                                     ggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGTGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                                <  2:1355475/150‑1 (MQ=255)
                                                                     ggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                                <  1:4513599/150‑1 (MQ=255)
                                                                     ggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTGCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                                <  1:628605/150‑1 (MQ=255)
                                                                          cGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCTACAAGGTACCCCACCCGCCTGGAACCGTGGCGATCCCGGGTGGCGCGGGTTTACCCCCGCACAGGCCCAGACAGGCGCTCAGTCAGTCgg                                                                           >  2:261346/1‑150 (MQ=255)
                                                                          cGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCACCTACAAGGTACCCCACCCGCCAGGAACCCTGGCCGTTGCGGTAGGAGCGGGTTTACCCGCGAACAGGCCCAGACAGGCGATCAGTAACTCgg                                                                           >  1:4076065/1‑150 (MQ=255)
                                                                           gCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgc                                                                          <  2:2021968/150‑1 (MQ=255)
                                                                              gatggCGTGAAGAGCCTGAAGGCCGCGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgc                                                                       <  1:2528069/147‑1 (MQ=255)
                                                                              gAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGTTACCCCCCCGGCCTGGAACCGTGGCGATCCCTGTAGTAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgc                                                                       <  1:3612647/150‑1 (MQ=255)
                                                                               aaGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgct                                                                      <  2:430240/150‑1 (MQ=255)
                                                                                   cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGaa                                                                  <  2:3475600/150‑1 (MQ=255)
                                                                                   cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGaa                                                                  <  2:3036550/150‑1 (MQ=255)
                                                                                   cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGaa                                                                  <  2:2807758/150‑1 (MQ=255)
                                                                                    gTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGGACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAac                                                                 <  2:2397779/150‑1 (MQ=255)
                                                                                     ttacgcgCCTGACGGCCGAGTAAACCCCTTCGCGGGTAACCCCGCTCCAACACGGTACCCCCCCGGCCTGGAACCGTGGCGAGCCCTGTAGGCGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTCGAaca                                                                <  1:107676/144‑1 (MQ=255)
                                                                                       aaGAGCCTGACGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAAcaca                                                              <  2:724669/150‑1 (MQ=255)
                                                                                             cTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAAc                                                        <  1:2870815/150‑1 (MQ=255)
                                                                                              tGACGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCAATCAGTCAGTCGGCTGCTTGAACACATAGAACa                                                       <  1:521214/150‑1 (MQ=255)
                                                                                              tGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACa                                                       <  2:802187/150‑1 (MQ=255)
                                                                                              tGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACa                                                       <  2:3071802/150‑1 (MQ=255)
                                                                                                   gCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTc                                                  <  2:3745984/150‑1 (MQ=255)
                                                                                                     cGAGTAAAACCCTTCGCGGTTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCgg                                                <  1:3105967/150‑1 (MQ=255)
                                                                                                     cGAGTAAAACCCTTCGCGGGTAACCCCGCTCATACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGGGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCgg                                                <  1:943451/150‑1 (MQ=255)
                                                                                                         tAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCa                                            <  2:2153610/150‑1 (MQ=255)
                                                                                                          aaaaCCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCAc                                           <  1:4430093/150‑1 (MQ=255)
                                                                                                              cccTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGAc                                       <  1:1427629/150‑1 (MQ=255)
                                                                                                              cccTTCGCGGGTAAACCCCCTCCTACAAGGTAACCCACCCGCCGGGAACCGTGGCGATCCCTGTTGGATCGGGTTTACCCGCGAACAGGCCCAGACACGCGATCAGTCAGTCGGCGGCTTGAACACAAAGAACAGGGTCGGCGCACAGAc                                       >  2:1861796/1‑150 (MQ=255)
                                                                                                                 TTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTTTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAg                                    >  1:4264263‑M1/1‑150 (MQ=255)
                                                                                                                  TCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAgg                                   >  1:3181724‑M1/1‑150 (MQ=255)
                                                                                                                   cgggggTAAACCCGATCCTACATGGTACCCCACCGGCCTGGAACCGTGGCGATCCGTGTAGGAGGGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGt                                  <  1:4359555/147‑1 (MQ=255)
                                                                                                                   CGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGt                                  >  2:2557291‑M1/1‑150 (MQ=255)
                                                                                                                      GGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAgg                                                                                          >  2:3490963‑M1/1‑91 (MQ=255)
                                                                                                                      GGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAgg                                                                                          <  1:3490963‑M1/91‑1 (MQ=255)
                                                                                                                             CCCGCTCCTACAATGTACCCCACCGTCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGAAAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACAAGGTAAATGGTccc                        <  1:1701887‑M1/150‑1 (MQ=255)
                                                                                                                             CCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTccc                        <  1:2474985‑M1/150‑1 (MQ=255)
                                                                                                                             CCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTccc                        <  1:3845965‑M1/150‑1 (MQ=255)
                                                                                                                             CCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTccc                        <  1:821421‑M1/150‑1 (MQ=255)
                                                                                                                              ccGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacagg                       <  2:753307‑M2/150‑106 (MQ=255)
                                                                                                                              CCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTcccc                       <  1:2229019‑M1/150‑1 (MQ=255)
                                                                                                                              CCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTcccc                       <  1:3783007‑M1/150‑1 (MQ=255)
                                                                                                                               CGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCg                      <  1:2798624‑M1/150‑1 (MQ=255)
                                                                                                                                  tCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggggcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgt                   <  2:3704235‑M2/150‑110 (MQ=255)
                                                                                                                                     tACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttgg                <  2:243651‑M2/150‑113 (MQ=255)
                                                                                                                                     TACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCTGGTAAATGGTCCCCGCTTCAt                <  2:4125101‑M1/150‑1 (MQ=255)
                                                                                                                                      aCAAGGTACCCAAACGGACTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaccaggccgttggt               <  2:3534228‑M2/150‑114 (MQ=255)
                                                                                                                                      ACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATc               <  2:4506487‑M1/150‑1 (MQ=255)
                                                                                                                                      ACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATc               >  2:2292332‑M1/1‑150 (MQ=255)
                                                                                                                                        AAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGACCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCa             <  1:2624841‑M1/150‑1 (MQ=255)
                                                                                                                                               cccACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAAt      >  2:396248/1‑150 (MQ=255)
                                                                                                                                                ccACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATc     >  1:518567/1‑150 (MQ=255)
                                                                                                                                                  aCCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCACTccc   >  2:3885616/1‑150 (MQ=255)
                                                                                                                                                  aCCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATccc   >  2:2476841/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATCCCC  >  NC_002947/5878849‑5879145

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: