Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,009,718:1 +A 6.2% intergenic (+78/‑94) 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,7181.A 6.2% 252.7 / 9.7 64intergenic (+78/‑94)PP_4412/PP_4413hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (31/29);  new base A (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

AAGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTTCTCTC  >  NC_002947/5009570‑5009866
                                                                                                                                                       |                                                                                                                                                    
aaGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATggcg                                                                                                                                                       >  6:197701/1‑149 (MQ=255)
   tgtTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCCTATGgcggcg                                                                                                                                                   >  8:123625/1‑149 (MQ=255)
             tGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCt                                                                                                                                         >  6:197139/1‑149 (MQ=255)
              gggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCt                                                                                                                                         >  7:164534/1‑148 (MQ=255)
               ggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTa                                                                                                                                       <  1:325844/149‑1 (MQ=255)
               ggCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGACGCCGCCATACGGCGAGCGGCTTggttca                                                                                                                                         <  1:47744/147‑5 (MQ=38)
                 cagcatgagGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGGGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCTGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTAc                                                                                                                                      <  1:92235/149‑1 (MQ=255)
                 cagcagcagGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACg                                                                                                                                     <  2:109661/149‑1 (MQ=255)
                        agGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGAGTCCTGGGCGCGCCCTTCGTGCCGGCTGTCGTGAACCAAGCATCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCtac                                                                                                                               >  7:352752/1‑146 (MQ=255)
                            ttCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATcc                                                                                                                          >  1:25576/1‑149 (MQ=255)
                             tCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGGGCCTAAGGCGGCGTGATTGCCCTTACCGGCttgccta                                                                                                                           >  1:125600/1‑143 (MQ=255)
                                gACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTcgc                                                                                                                      <  6:80935/149‑1 (MQ=255)
                                  cGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGt                                                                                                                    <  8:164534/149‑1 (MQ=255)
                                   gATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGGGGTCGTGAACCAAGCCCCACGCCGTATGGCCGCGTCTTGGCCCTTACGGGCTTCCAACCTCGCGtc                                                                                                                   >  4:192151/1‑149 (MQ=255)
                                            cctcctCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTc                                                                                                          >  7:223186/1‑149 (MQ=255)
                                             ctcctcCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCACGCGTCTCCGCCTTcc                                                                                                         >  3:240639/1‑149 (MQ=255)
                                                 tccAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGc                                                                                                     >  4:298839/1‑149 (MQ=255)
                                                    aaaCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgc                                                                                                  <  4:240639/149‑1 (MQ=255)
                                                     aaCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGgcggcg                                                                                                                                                   <  6:240246/99‑1 (MQ=255)
                                                     aaCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGgcggcg                                                                                                                                                   >  5:240246/1‑99 (MQ=255)
                                                      aCGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAACAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCGCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgcgc                                                                                                <  5:197701/149‑1 (MQ=255)
                                                        ggggAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTgcgcgc                                                                                              >  5:225086/1‑149 (MQ=255)
                                                        ggggAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTAGGGCCGGCGGTCGTGAACCAATCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTAAGGGCTTGCATACTAGCCGCGCCGCCTTCCTGGCTtcgcgc                                                                                              >  3:214860/1‑149 (MQ=255)
                                                              ctctTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCt                                                                                        >  5:353859/1‑149 (MQ=255)
                                                                 ttctAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGGCTCCGCCTTCCTGGCTGCGCGCTCCGCGTgg                                                                                     >  4:43005/1‑148 (MQ=255)
                                                                           gCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAAc                                                                           >  6:380728/1‑149 (MQ=255)
                                                                              aGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGt                                                                         >  2:122885/1‑148 (MQ=255)
                                                                                  ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCaaa                                                                     <  3:64587/148‑1 (MQ=255)
                                                                                  ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCaaa                                                                     <  3:192151/148‑1 (MQ=255)
                                                                                  ctgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAAt                                                                    >  7:13230/1‑149 (MQ=255)
                                                                                   tgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGGTCTGCCCAACGTCAAAt                                                                    >  1:349009/1‑148 (MQ=255)
                                                                                   tgcAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAAt                                                                    >  6:350412/1‑148 (MQ=255)
                                                                                     cAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGTTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  <  1:297911/148‑1 (MQ=255)
                                                                                     cAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  <  1:170567/148‑1 (MQ=255)
                                                                                     cAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAg                                                                 <  8:223186/149‑1 (MQ=255)
                                                                                      aaaaGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGc                                                                >  6:169458/1‑149 (MQ=255)
                                                                                        atgCGTCCTGGGCGGGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGTGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGt                                                              <  8:352752/147‑1 (MQ=255)
                                                                                        aaGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCCTCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  >  8:329286/1‑145 (MQ=255)
                                                                                        aaGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGa                                                                  <  7:329286/145‑1 (MQ=255)
                                                                                        aaGCGTCCTGGGCGCGCCCATCGGGCCGGCTGTGGTGAAACAACCCGCACGCCGTATTGGGGCGTCGTTGCCCGTTAGGGCTTGCCTACTAGCGGGTCCGCCGTGCTGGCATCGCGCGCCGCGTGCAAGGGCCCACGTCAAATGAGCgg                                                              >  4:370784/1‑148 (MQ=255)
                                                                                         aGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGtc                                                             <  5:380728/149‑1 (MQ=255)
                                                                                          gCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGGCTCCGCCTGCCTGGCTGCGCGCTCCGCTTGCTCTGCCCACCGTCAAATGCGCGtct                                                            >  6:192398/1‑149 (MQ=255)
                                                                                               cTGGGCGCGCGCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTGTGCCCAACGTCAAATGAGCGtctttctt                                                       <  3:266226/149‑1 (MQ=255)
                                                                                               cTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGtctttctt                                                       <  7:123625/149‑1 (MQ=255)
                                                                                                  ggCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTg                                                      <  4:139119/147‑1 (MQ=255)
                                                                                                   gcgcgcCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGc                                                    >  8:133350/1‑148 (MQ=255)
                                                                                                   gcgcgcCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCa                                                   >  8:153118/1‑149 (MQ=255)
                                                                                                     gcgcCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATc                                                 <  3:43005/149‑1 (MQ=255)
                                                                                                      cgcCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCt                                                <  7:133350/149‑1 (MQ=255)
                                                                                                       gcCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTg                                               >  4:195759/1‑149 (MQ=255)
                                                                                                          cTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCAt                                            <  3:298839/149‑1 (MQ=255)
                                                                                                             cGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATg                                           <  2:25576/147‑1 (MQ=255)
                                                                                                                      cTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGtt                                 <  4:27605/149‑1 (MQ=255)
                                                                                                                        gTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTaa                               >  7:322814/1‑149 (MQ=255)
                                                                                                                        gTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTaa                               >  7:330225/1‑149 (MQ=255)
                                                                                                                         tCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAAc                              <  5:197139/149‑1 (MQ=255)
                                                                                                                            tGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCa                           >  7:270477/1‑149 (MQ=255)
                                                                                                                                ccAAGACGCCGCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGTCCAACGTCAAATGAGCGTCTTTATTGGCATCTGCATGCCACCCTGTTAAATCActt                        >  2:336138/1‑148 (MQ=255)
                                                                                                                                ccAAGACGCCGCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtc                       >  4:266954/1‑149 (MQ=255)
                                                                                                                                 cAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtct                      <  8:330225/149‑1 (MQ=255)
                                                                                                                                  aaGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtctt                     >  2:58445/1‑149 (MQ=255)
                                                                                                                                  aaGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCActtctt                     >  7:172753/1‑149 (MQ=255)
                                                                                                                                   aGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCttt                    <  6:225086/149‑1 (MQ=255)
                                                                                                                                    gCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTg                   <  1:58445/149‑1 (MQ=255)
                                                                                                                                    agacgccGCCATACGGCGAGCGGCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCttt                    <  3:266954/143‑1 (MQ=255)
                                                                                                                                            gCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTa            <  7:153118/149‑1 (MQ=255)
                                                                                                                                            gCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTa            <  3:370784/149‑1 (MQ=255)
                                                                                                                                               gTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGcc         >  2:69072/1‑149 (MQ=255)
                                                                                                                                                    gcgtcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtct     <  8:322814/148‑1 (MQ=255)
                                                                                                                                                    gcggcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtct     <  1:69072/148‑1 (MQ=255)
                                                                                                                                                     cggcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtctc    >  4:35082/1‑148 (MQ=255)
                                                                                                                                                      ggcgTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTtctctc  <  6:353859/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
AAGTGTTTGGCTCTGGGCAGCAGCAGGATTCAGACGATGAACCACCTCCTCCAAACGGGGAACTCTTCTAGGTAGGCCAGTGCTGCAAAAGCGTCCTGGGCGCGCCCTTCGGGCCGGCTGTCGTGAACCAAGCCGCTCGCCGTATGGCGGCGTCTTGGCCCTTACGGGCTTCCATCCTCGCGTCTCCGCCTTCCTGGCTGCGCGCTCCGCTTGCTCTGCCCAACGTCAAATGAGCGTCTTTCTTGGCATCTGCATGCCACCCTGTTAACTCACTTCTTTGCCGCCTAGCCTTCTCTC  >  NC_002947/5009570‑5009866

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: