Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,009,711 G→A 6.0% intergenic (+71/‑101) PP_4412 → / → PP_4413 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,009,7110GA 6.0% 194.9 / 2.0 68intergenic (+71/‑101)PP_4412/PP_4413hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GAGGAGCAAGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCC  >  NC_002947/5009563‑5009859
                                                                                                                                                      |                                                                                                                                                     
gaggagCAAGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCg                                                                                                                                                       <  1:362263/149‑1 (MQ=255)
 aggagCAAGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGt                                                                                                                                                      <  3:144324/149‑1 (MQ=255)
       aaGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATggcg                                                                                                                                                >  6:197701/1‑149 (MQ=255)
          tgtTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCCTATGgcg‑gcg                                                                                                                                            >  8:123625/1‑149 (MQ=255)
                    tGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCt                                                                                                                                  >  6:197139/1‑149 (MQ=255)
                     gggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCt                                                                                                                                  >  7:164534/1‑148 (MQ=255)
                      ggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTa                                                                                                                                <  1:325844/149‑1 (MQ=255)
                      ggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGACGC‑C‑GCCATACGGCGAGCGGCTTggttca                                                                                                                                  <  1:47744/147‑5 (MQ=38)
                        cagcatgagGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGGGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCTGGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTAc                                                                                                                               <  1:92235/149‑1 (MQ=255)
                        cagcagcagGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACg                                                                                                                              <  2:109661/149‑1 (MQ=255)
                               agGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGAGTCCTGGGCGCGCCCTTCGTGCC‑GGCTGTCGTGAACCAAGCATCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCtac                                                                                                                        >  7:352752/1‑146 (MQ=255)
                                   ttCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATcc                                                                                                                   >  1:25576/1‑149 (MQ=255)
                                    tCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTCGGGCCTAAGGCG‑GCGTGATTGCCCTTACCGGCttgccta                                                                                                                    >  1:125600/1‑143 (MQ=255)
                                       gACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTcgc                                                                                                               <  6:80935/149‑1 (MQ=255)
                                         cGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGt                                                                                                             <  8:164534/149‑1 (MQ=255)
                                          gATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGGGGTCGTGAACCAAGCCCCAC‑GCCGTATGGCC‑GCGTCTTGGCCCTTACGGGCTTCCAACCTCGCGtc                                                                                                            >  4:192151/1‑149 (MQ=255)
                                                   cctcctCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTc                                                                                                   >  7:223186/1‑149 (MQ=255)
                                                    ctcctcCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCACGCGTCTCCGCCTTcc                                                                                                  >  3:240639/1‑149 (MQ=255)
                                                        tccAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGc                                                                                              >  4:298839/1‑149 (MQ=255)
                                                           aaaCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgc                                                                                           <  4:240639/149‑1 (MQ=255)
                                                            aaCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGgcg‑gcg                                                                                                                                            <  6:240246/99‑1 (MQ=255)
                                                            aaCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGgcg‑gcg                                                                                                                                            >  5:240246/1‑99 (MQ=255)
                                                             aCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAACAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCGC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgcgc                                                                                         <  5:197701/149‑1 (MQ=255)
                                                               ggggAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgcgcgc                                                                                       >  5:225086/1‑149 (MQ=255)
                                                               ggggAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTAGGGCC‑GGCGGTCGTGAACCAATCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTAAGGGCTTGCATACTAGCCGCGCCGCCTTCCTGGCTtcgcgc                                                                                       >  3:214860/1‑149 (MQ=255)
                                                                     ctctTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCt                                                                                 >  5:353859/1‑149 (MQ=255)
                                                                        ttctAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGGCTCCGCCTTCCTGGCTGCGCGCTCCGCGTgg                                                                              >  4:43005/1‑148 (MQ=255)
                                                                                  gCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAAc                                                                    >  6:380728/1‑149 (MQ=255)
                                                                                     aGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGt                                                                  >  2:122885/1‑148 (MQ=255)
                                                                                         ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCaaa                                                              <  3:64587/148‑1 (MQ=255)
                                                                                         ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCaaa                                                              <  3:192151/148‑1 (MQ=255)
                                                                                         ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAAt                                                             >  7:13230/1‑149 (MQ=255)
                                                                                          tgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGGTCTGCCCAACGTCAAAt                                                             >  1:349009/1‑148 (MQ=255)
                                                                                          tgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAAt                                                             >  6:350412/1‑148 (MQ=255)
                                                                                            cAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGTTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                           <  1:297911/148‑1 (MQ=255)
                                                                                            cAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                           <  1:170567/148‑1 (MQ=255)
                                                                                            cAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAg                                                          <  8:223186/149‑1 (MQ=255)
                                                                                             aaaaGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGc                                                         >  6:169458/1‑149 (MQ=255)
                                                                                               atgCGTCCTGGGCGGGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGTGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGt                                                       <  8:352752/147‑1 (MQ=255)
                                                                                               aaGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCCTCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                           >  8:329286/1‑145 (MQ=255)
                                                                                               aaGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                           <  7:329286/145‑1 (MQ=255)
                                                                                               aaGCGTCCTGGGCGCGCCCATCGGGCC‑GGCTGTGGTGAAACAACCCGCAC‑GCCGTATTGGG‑GCGTCGTTGCCCGTTAGGGCTTGCCTACTAGCGGGTCCGCCGTGCTGGCATCGCGCGCCGCGTGCAAGGGCCCACGTCAAATGAGCgg                                                       >  4:370784/1‑148 (MQ=255)
                                                                                                aGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGtc                                                      <  5:380728/149‑1 (MQ=255)
                                                                                                 gCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGGCTCCGCCTGCCTGGCTGCGCGCTCCGCTTGCTCTGCCCACCGTCAAATGCGCGtct                                                     >  6:192398/1‑149 (MQ=255)
                                                                                                      cTGGGCGCGCGCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTGTGCCCAACGTCAAATGAGCGtctttctt                                                <  3:266226/149‑1 (MQ=255)
                                                                                                      cTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGtctttctt                                                <  7:123625/149‑1 (MQ=255)
                                                                                                         ggCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTg                                               <  4:139119/147‑1 (MQ=255)
                                                                                                          gcgcgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGc                                             >  8:133350/1‑148 (MQ=255)
                                                                                                          gcgcgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCa                                            >  8:153118/1‑149 (MQ=255)
                                                                                                            gcgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATc                                          <  3:43005/149‑1 (MQ=255)
                                                                                                             cgcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCt                                         <  7:133350/149‑1 (MQ=255)
                                                                                                              gcCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTg                                        >  4:195759/1‑149 (MQ=255)
                                                                                                                 cTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCAt                                     <  3:298839/149‑1 (MQ=255)
                                                                                                                    cGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATg                                    <  2:25576/147‑1 (MQ=255)
                                                                                                                             cTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGtt                          <  4:27605/149‑1 (MQ=255)
                                                                                                                               gTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTaa                        >  7:322814/1‑149 (MQ=255)
                                                                                                                               gTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTaa                        >  7:330225/1‑149 (MQ=255)
                                                                                                                                tCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAAc                       <  5:197139/149‑1 (MQ=255)
                                                                                                                                   tGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCa                    >  7:270477/1‑149 (MQ=255)
                                                                                                                                       ccAAGACGC‑C‑GCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGTCCAACGTCAAATGAGCGTCTTTATTGGCATCTGCATGCCACCCTGTTAAATCActt                 >  2:336138/1‑148 (MQ=255)
                                                                                                                                       ccAAGACGC‑C‑GCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtc                >  4:266954/1‑149 (MQ=255)
                                                                                                                                        cAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtct               <  8:330225/149‑1 (MQ=255)
                                                                                                                                         aaGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtctt              >  7:172753/1‑149 (MQ=255)
                                                                                                                                         aaGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtctt              >  2:58445/1‑149 (MQ=255)
                                                                                                                                          aGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCttt             <  6:225086/149‑1 (MQ=255)
                                                                                                                                           gCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTg            <  1:58445/149‑1 (MQ=255)
                                                                                                                                           agacgcc‑GCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCttt             <  3:266954/143‑1 (MQ=255)
                                                                                                                                                   gCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTa     <  7:153118/149‑1 (MQ=255)
                                                                                                                                                   gCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTa     <  3:370784/149‑1 (MQ=255)
                                                                                                                                                      gTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGcc  >  2:69072/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GAGGAGCAAGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCC‑GGCTGTCGTGAACCAAGCCGCTC‑GCCGTATGGCG‑GCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCC  >  NC_002947/5009563‑5009859

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑