Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,538,604 A→T 10.4% intergenic (+53/‑121) PP_4020 → / → cpo oxidoreductase/non‑heme chloroperoxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,538,6040AT10.4% 529.9 / 11.6 194intergenic (+53/‑121)PP_4020/cpooxidoreductase/non‑heme chloroperoxidase
Reads supporting (aligned to +/- strand):  ref base A (86/80);  new base T (2/22);  total (91/103)
Fisher's exact test for biased strand distribution p-value = 3.94e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

AGGCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGCGT  >  NC_002947/4538455‑4538753
                                                                                                                                                     |                                                                                                                                                     
    aTCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGGGCGGGGCTCTTAGCGGGGGGAAGGGACCGcacacagg                                                                                                                                                   >  2:1789935/1‑146 (MQ=255)
      cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                                 <  1:4608076/150‑1 (MQ=255)
      cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                                 <  1:1897065/150‑1 (MQ=255)
      cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                                 <  2:4263022/150‑1 (MQ=255)
          aaCGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGc                                                                                                                                             <  1:1516089/150‑1 (MQ=255)
          aaCGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGc                                                                                                                                             <  2:1028526/150‑1 (MQ=255)
            cGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCgg                                                                                                                                           <  1:1707565/150‑1 (MQ=255)
            cGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCgg                                                                                                                                           <  2:3702561/150‑1 (MQ=255)
             ggCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCggg                                                                                                                                          <  2:6819224/150‑1 (MQ=255)
             ggCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCGCTTCGCGGGGGGAAGGGGCCTCACAGGGTAGGTGCggg                                                                                                                                          >  2:4149695/1‑150 (MQ=255)
                cGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGGCGCATCGCGGGTGAAAGGGGCCCCCCCGGGGAGGAGCGGGGTc                                                                                                                                       >  2:1387242/1‑150 (MQ=255)
                   ttCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCTTGCGCGACCGTGTCTGAAGGGCATCAGTGTTGGCGGTGCGGGGCGCGTCGCGGGGGGGAGGGGCCGAACACTGTAGGGGGGGGGTcgac                                                                                                                                    >  1:5485136/1‑147 (MQ=255)
                   ttCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGGCGGTGGGGGCCTCGTAGCGGGGGGACGGGACCGCACACAGTAGGAGCGGGGTCAcc                                                                                                                                    >  1:5691794/1‑150 (MQ=255)
                    tCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCt                                                                                                                                   <  2:1309330/150‑1 (MQ=255)
                    tCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCt                                                                                                                                   <  1:6126097/150‑1 (MQ=255)
                      cccGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgc                                                                                                                                 >  1:4559985‑M2/1‑148 (MQ=255)
                       ccGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGTGCGGGGCGCATCGCGGGGGGAAGGGACCACACAGTGTAGGAGCGGGTTAACCAgcg                                                                                                                                <  2:260801‑M2/150‑4 (MQ=255)
                        cggcggCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGGCGGTGGGGGCCACGTCGCGGGGGTCAGGGGCCGCACAGAGTAGGAGCGGGTTAACCAgcga                                                                                                                               <  1:526183‑M2/150‑5 (MQ=255)
                        cggcggCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcga                                                                                                                               >  2:716869‑M2/1‑146 (MQ=255)
                         ggcggcGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaa                                                                                                                              >  1:1617338‑M2/1‑145 (MQ=255)
                         ggcggcGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaa                                                                                                                              >  1:5640474‑M2/1‑145 (MQ=255)
                          gcggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCggg                                                                                                                                          >  2:4667420/1‑137 (MQ=255)
                          gcggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaag                                                                                                                             >  1:2953564‑M2/1‑144 (MQ=255)
                            ggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCACAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagag                                                                                                                           >  1:4519752‑M2/1‑142 (MQ=255)
                            ggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagag                                                                                                                           >  2:4210096‑M2/1‑142 (MQ=255)
                            catggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                           >  1:3334493‑M1/79‑150 (MQ=255)
                             gcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagagg                                                                                                                          >  2:6693961‑M2/1‑141 (MQ=255)
                             atggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGTGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAgg                                                                                                                          >  2:5331137‑M1/78‑150 (MQ=255)
                                  ctgctgGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGGCGCGTCGCGGGTGGAGGGGGCCCCACAGGGGAGGGGCGGGGTGACCGGCGAAGAGGccccc                                                                                                                     >  2:2289414/1‑148 (MQ=255)
                                     gtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCcacaca                                                                                                                  >  2:4351460‑M1/70‑150 (MQ=255)
                                      tggctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAg                                                                                                                 <  2:5956840/150‑1 (MQ=255)
                                      tggccgCCCCTTTGGAAGCTGCGGCAAACGCCGGGCACCGCACGCGCGACCGGGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAg                                                                                                                 <  2:6537953/150‑1 (MQ=255)
                                        gctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGc                                                                                                               <  1:1672838/150‑1 (MQ=255)
                                          tgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCaa                                                                                                             <  1:2069970/150‑1 (MQ=255)
                                          tgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGCCTCTTCGCGGGGGGAAGGGCCCGCCCCGTGGAGGTGCGGGGTGACCAGCGACGAGGCCCACACAGGCac                                                                                                             >  2:3182334/1‑149 (MQ=255)
                                            cgcCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCTGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAACGGGACCTCACTGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAaca                                                                                                           <  2:4118444/150‑1 (MQ=255)
                                             gcCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAacac                                                                                                          <  2:4652591/150‑1 (MQ=255)
                                               cTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGTGCGGGCCGCGTCGCGGGGGGAACGGACCGCACAGGGTGGGGGCGGGGGGACCAGCGAAGACGCCCCCACACCCcacacac                                                                                                        >  2:3758518/1‑149 (MQ=255)
                                                tttGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAAc                                                                                                       <  2:443494/150‑1 (MQ=255)
                                                tttGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGTGGGGGGCGCGTCGCGGGTGGAAGGGCCCCCACACTGTGGGGGCGGGGTGGACGGCGAAGGGGCCCCCACCCGCCACCCAAc                                                                                                       >  2:1842787/1‑150 (MQ=255)
                                                   ggAACCTGCGGCCCCGGCCGGGCACCGCACGCGCGCCCGGGTCTGACGGGCATCAGTGTTGCCTGTGCTTGCCTGTTCGCGGGTGAACGGTAAACAACAGTGTTTGGGCGGGTTCACCGGCGAAGAGGCCCACACAGGCAACACAACTTc                                                                                                    <  2:3408863/150‑1 (MQ=255)
                                                       gCTGCGGCAACGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTGGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGa                                                                                                <  2:6614847/150‑1 (MQ=255)
                                                             gCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGCCACGTCGCGGGTGGGAGGGGCCGCCCACTGTATGAGCCGGGTGAACAGCGAACAGGCCCACCCACGCCACACAACATACAGATAcccc                                                                                          >  2:4761226/1‑149 (MQ=255)
                                                                 cggCCTGGTACCGCTTGCGCGACCGTGTCTGAAGGGCATCTGTGTTGCGTGTGGGGGCCTCTTCGGGGGTGAAAGGGACCTCACTGTTTAGGAGGGGGTTCACCTGCGATGAGGCCCACACTGGCAACACAACTTCCTGATACCCGCCaa                                                                                      <  2:2370958/149‑1 (MQ=255)
                                                                     cTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGa                                                                                                                                               <  1:4905742/89‑1 (MQ=255)
                                                                     cTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGa                                                                                                                                               >  2:4905742/1‑89 (MQ=255)
                                                                       ggTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGTGCGGGCCGCGTCGCGGGTGAAAGGGACCGCACAGTGTAGGAGCGGGTTCACCGGCGAAGAGGCCCACACAGGCAACACCACTTCCTGGTACCCGCCAAAATacc                                                                                >  2:3976690/1‑150 (MQ=255)
                                                                         tcccGCAGGCGCGCCCGGGTCTGACGGGCATCCGTGTGGCCTGTGCGGGCCTCTGCGGGGTTGAAAGGTACCTCACAGTGTTGGCGCGGTTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCcctc                                                                              <  2:2486006/148‑1 (MQ=255)
                                                                           ccGCATGCGCGCCCGTGTCTGAAGGGCATCAGTGGTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCcctcct                                                                            <  2:1200469/150‑1 (MQ=255)
                                                                                tGCGCGACGTTGTCTTATGGGATTCTGTGTTGCGTGTGGGGGCCTGTTCGGGGGTGAATGGGACCTCACTGTTTTGGTGGGGTTTCACCTGCGATGAGGCCAAAACTGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATtct                                                                       <  1:4441860/150‑1 (MQ=255)
                                                                                tGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCAGTGCGGGCCTCTTCGCGGGTGAAAGGGACCGCACAGTGTAGGAGCGGGGTCACCAGCGAAGAGGCCCACACAGGCAACCCAACATCCAGATACCCGCCAAAATCCCACCGAttga                                                                       >  2:2806678/1‑148 (MQ=255)
                                                                                          tgtCTTACGGGCATCTGTGTTGCCTGTGGGGGCCTCTTGGCGGTTTAATGGGACCTCAATTTGTTGGAGGGGGTTCACCTGCTAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTcc                                                             <  2:1929453/150‑1 (MQ=255)
                                                                                           gtgtGACGGGCATCTGTTTTGCGTGTGGGGGCCTGTTCGGGGGTTAATGGGACCTCACTTTTTTGGCGGGGTTTCACCTGCTAAGTGGCCCACACTGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTccc                                                            <  1:3358905/147‑1 (MQ=255)
                                                                                             tttacgggAATCTGTGTTGCGTGTGGGGGCCTCTTGGGGGTTGAATGGGACCTAACTGTTTTGGCGGGGGTTCACCTGCGAAGAGGCCCACACTGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCtt                                                          <  1:2705160/145‑1 (MQ=255)
                                                                                                   ggcttatgtgtTGCCTTTGGGGGCTTTTGCGGGGGTTAATGGCCCCTCACTGTGTTGGGGGGGTTTCCCGTGCCCAGCGCCCCACACAGCCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTc                                                    <  1:1248077/143‑1 (MQ=255)
                                                                                                             ttGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagaggcccacacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgca                                          <  1:4861037‑M2/150‑90 (MQ=255)
                                                                                                              TGCCTGTGCGGGCCTCTTCGGGGGTGAATGGTACCTAACTGTGTAGGAGCGGGTTCACCTGCTAAGAGGCCAACACAGGAAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAgg                                         >  1:5627633‑M1/1‑150 (MQ=255)
                                                                                                               GCCTGTGCGGGCCTCTTCGCGGGTTAAAGGGACCTCACAGTGTTGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGga                                        >  2:679997‑M1/1‑150 (MQ=255)
                                                                                                               GCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTTGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGga                                        >  2:4090845‑M1/1‑150 (MQ=255)
                                                                                                                    tgCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagag                                                                                                                           <  2:5799787‑M2/62‑9 (MQ=255)
                                                                                                                    TGCGGGCCTCTTCGCGGGTGAAAGGGCCCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGacac                                   >  2:747362‑M1/1‑150 (MQ=255)
                                                                                                                      cGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGa                                 >  2:3007898/1‑150 (MQ=255)
                                                                                                                       gggCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagaggcccgcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggt                                <  1:5189965‑M2/150‑100 (MQ=255)
                                                                                                                         gCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCa                              >  2:5625256/1‑150 (MQ=255)
                                                                                                                         gCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCa                              >  1:4795869/1‑150 (MQ=255)
                                                                                                                          ccTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagaggcccacacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgc                             <  2:3622630‑M2/150‑103 (MQ=255)
                                                                                                                           ctctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATg                            >  2:6443744/1‑150 (MQ=255)
                                                                                                                            tctTCGCGGGTTAAAGGTACCTCACTGTGTAGGAGCGGTTTCACCGGCTAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                           <  2:2460061/150‑1 (MQ=255)
                                                                                                                            tctTCGCGGGTGAAAGGTACCTCACTGTTTTGGTGGGGGTTCACCTGCTAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                           <  2:1902751/150‑1 (MQ=255)
                                                                                                                            tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagaggcccgcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccg                           <  1:1940156‑M2/150‑105 (MQ=255)
                                                                                                                            tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagaggcccacacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccg                           <  2:1109402‑M2/150‑105 (MQ=255)
                                                                                                                            tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGAcc                               >  1:5913383/1‑146 (MQ=255)
                                                                                                                            tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                           >  2:5551841/1‑150 (MQ=255)
                                                                                                                            tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                           >  1:417754/1‑150 (MQ=255)
                                                                                                                              ttCGCGGGTTAATGGGACCTCACAGTTTAGGAGCGGGTTCACCTGCTAAGTGGCCCACACTGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                         <  1:2790386/150‑1 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTgcgaagaggcccgcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgcc                         <  1:5342498‑M2/150‑107 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCCAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                         >  1:13618/1‑150 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                         >  1:620723/1‑150 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                         >  1:907615/1‑150 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                         <  2:6198002/150‑1 (MQ=255)
                                                                                                                                 gggggTTAAGGGGACCACACTGTTTTGGGGGGGGTACCCCTGCTACGGGGCCCACACCGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                      <  1:6555334/148‑1 (MQ=255)
                                                                                                                                 gcgGGTGAAAGGGACCTCACAGTGTTGGAGGGGGTTCACCTGCGAAGTGGCCAACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                      <  1:2145253/150‑1 (MQ=255)
                                                                                                                                      tGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCAGGAGCTATGTAAc                 >  1:5775439/1‑150 (MQ=255)
                                                                                                                                        aatgggACCTCACTGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCa               <  2:5485136/147‑1 (MQ=255)
                                                                                                                                        aaaGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCa               <  1:5551841/150‑1 (MQ=255)
                                                                                                                                          aGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACg             >  1:6426536/1‑150 (MQ=255)
                                                                                                                                           gggACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGa            <  2:3865483/150‑1 (MQ=255)
                                                                                                                                             gACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAg          >  2:1509451/1‑150 (MQ=255)
                                                                                                                                              aCCTCACAGTGTAGGAGCGGGTTCACCTgcgaagaggcccgcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgccggggaagcggccgttg         >  1:3129479‑M2/1‑28 (MQ=255)
                                                                                                                                              aCCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAgg         >  1:1130576/1‑150 (MQ=255)
                                                                                                                                                 tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                         <  1:3753256/131‑1 (MQ=255)
                                                                                                                                                 tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACg      >  1:241055/1‑150 (MQ=255)
                                                                                                                                                  cacaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACgg     <  2:6389057/150‑1 (MQ=255)
                                                                                                                                                  cacaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACgg     >  2:3256301/1‑150 (MQ=255)
                                                                                                                                                  cacaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACgg     >  1:729500/1‑150 (MQ=255)
                                                                                                                                                  cacaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACgg     >  2:2721522/1‑150 (MQ=255)
                                                                                                                                                   acaGTGTAGGAGCGGGTTCACCTgcgaagaggcccacacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgccggggaagcggccgttgtcgat    >  2:5528072‑M2/1‑23 (MQ=255)
                                                                                                                                                   acaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGc    <  2:2181954/150‑1 (MQ=255)
                                                                                                                                                     aGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGCGt  >  1:305634/1‑150 (MQ=255)
                                                                                                                                                     aGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGCGt  >  1:6755551/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGCGT  >  NC_002947/4538455‑4538753

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: