Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,538,604 A→T 7.6% 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,6040AT 7.6% 646.3 / 7.6 234intergenic (+53/‑121)PP_4020/cpooxidoreductase/non‑heme chloroperoxidase
Reads supporting (aligned to +/- strand):  ref base A (82/125);  new base T (3/20);  total (89/145)
Fisher's exact test for biased strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

AGGCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGC  >  NC_002947/4538455‑4538751
                                                                                                                                                     |                                                                                                                                                   
aGGCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTcaca                                                                                                                                                     <  1:2014365/150‑1 (MQ=255)
aGGCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTcaca                                                                                                                                                     <  2:4098876/150‑1 (MQ=255)
  gCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAgt                                                                                                                                                   <  1:2367488/150‑1 (MQ=255)
  gCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAgt                                                                                                                                                   <  2:6316962/150‑1 (MQ=255)
   cATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCAGAgtg                                                                                                                                                  <  1:1002376/150‑1 (MQ=255)
     tCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTa                                                                                                                                                <  1:8605014/150‑1 (MQ=255)
      cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                               <  2:2638135/150‑1 (MQ=255)
      cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                               <  1:2524071/150‑1 (MQ=255)
       gACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGGCGGTGCGGGCCGCGTAGCGGGTGGAAGGGGCCGCACAGTGTAgg                                                                                                                                              >  2:7426238/1‑150 (MQ=255)
             ggCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCggg                                                                                                                                        <  1:6633670/150‑1 (MQ=255)
             ggCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGGGAAAGGGCCCGCACAGTGTAGGAGCggg                                                                                                                                        >  2:5746113/1‑150 (MQ=255)
                cGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCGCGTCGCGGGGGGAAGGGGCCTCACAGGGTAGGAGCGGGGTc                                                                                                                                     >  2:3725915/1‑150 (MQ=255)
                 gCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGCCGCGTCGCGGGGGGAAGGGCCCGCACCGGGGTGGAGCGGGGTCa                                                                                                                                    >  2:6092374/1‑150 (MQ=255)
                  cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGGGGGCCGCTTAGCGGGGGAAAGGGACCACAAAATGTTGGAGCGGGTTCAc                                                                                                                                   >  1:7188315/1‑150 (MQ=255)
                  cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCAc                                                                                                                                   <  2:4151926/150‑1 (MQ=255)
                  cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCAc                                                                                                                                   <  1:4623575/150‑1 (MQ=255)
                  cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGGGTAAGGGCCCTCACACGGTAGGTGCGGGTTCAc                                                                                                                                   >  2:8554540/1‑150 (MQ=255)
                  cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGTGCGGGCCACATCGCGGGGGAAAAGGACCACACACTGTAGGAGCGGGTTGAc                                                                                                                                   >  1:7939560/1‑150 (MQ=255)
                    tCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCt                                                                                                                                 <  2:2979965/150‑1 (MQ=255)
                      atagctcatggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGc                                                                                                                               >  2:5611910‑M1/85‑150 (MQ=255)
                         ggcggcGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGaa                                                                                                                            <  2:2926185/150‑1 (MQ=255)
                           tcatggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAaga                                                                                                                          >  1:4014117‑M1/80‑150 (MQ=255)
                           tcatggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAaga                                                                                                                          <  2:3949477‑M1/71‑1 (MQ=255)
                           cggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAaga                                                                                                                          <  2:6774789/150‑1 (MQ=255)
                            ggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                         <  1:6277897/150‑1 (MQ=255)
                               gACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAgtg                                                                                                                                                  >  1:1991996/1‑122 (MQ=255)
                                 ccTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCca                                                                                                                    <  1:6866618/150‑1 (MQ=255)
                                     ctggctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCATCGCGGGTGAAAGGGACCGCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcaca                                                                                                                <  2:8710538‑M2/150‑5 (MQ=255)
                                      tctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACTGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAg                                                                                                               >  1:3786242‑M1/69‑150 (MQ=255)
                                       ggctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacagg                                                                                                              >  1:1171411‑M2/1‑144 (MQ=255)
                                       ggctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacagg                                                                                                              >  1:8397558‑M2/1‑144 (MQ=255)
                                       cgcctgggtgcaggacgcgaagggagctattgagaccagcagggattgtggcgggtcacaggaagttGTGGTGCCTGGGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAgg                                                                                                              >  1:6471849‑M1/68‑150 (MQ=255)
                                        gctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggc                                                                                                             <  1:7708403‑M2/150‑8 (MQ=255)
                                            gggtgcaggaagcgaagggagctattgagaataggagggattttggcgggcatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTCGGTGCGGTTTCACCTGCGAAGAGGCCCACACAGGCAaca                                                                                                         >  2:443701‑M1/63‑150 (MQ=255)
                                            cgcCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaaca                                                                                                         >  2:1426648‑M2/1‑139 (MQ=255)
                                             gcctacggAAGCAGCGCCAAAGGCCCGGTACCCCAAGCCCGACCGTGTCACAACGGCATCAGAGTTGCCTGTGCGGGCCACTTCGCGGGAGAAACGGACCTCACAGTCTAGGAGCCCGTTCACCTGCGAAGAGGCCCACACAGGCAacac                                                                                                        <  2:5083270/144‑1 (MQ=255)
                                               cTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCTGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACaa                                                                                                      <  2:6549486/150‑1 (MQ=255)
                                               cTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGTTTCACCTGCGAAGAGGCCCACACAGGCAACACaa                                                                                                      <  2:3467586/150‑1 (MQ=255)
                                                  aggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGCGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAgg                                                                                                                        >  1:2732896‑M1/57‑129 (MQ=255)
                                                  aggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGCGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAgg                                                                                                                        <  2:2732896‑M1/73‑1 (MQ=255)
                                                    gAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTcc                                                                                                 <  2:6553683/150‑1 (MQ=255)
                                                     aaGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCATAGCGGGTGAAAGGGACCTCACACGGTAGGAGCGGGTTGACCAGCGAAGAGGCCCACACAGGCCACACAACttaca                                                                                                >  2:4998332/1‑147 (MQ=255)
                                                       gcggCGGCAACGGCCGGGTACCGCAGGCGCGACGGTGTCTGACGGGCATCAGTGTTGCCTGTGGGGGCCTCTTCGGGGGTGAAAGGTACCTCACTGTGTGGGCGCGGGTTCACGTGCGATGTGGCCCACACTGGCAACACAACTTCCTGa                                                                                              <  2:2327298/147‑1 (MQ=255)
                                                               aaaGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGCCACATCCCCGGGGGAAGGGCCCCCACCGGGGCGGGGCGGGGTCACCAGCGCAGAGGCCCCCACCCGCAACACCACAACCAGAAAcccccc                                                                                      >  2:2664286/1‑150 (MQ=255)
                                                                   gCCTGGTACCGCATGCGCGACCGTTTCTGACGGGCACCTGTGTTGCCTGTGCGGGCCTCTGCGGGGGTGAAAGGGACCTCACTGTGTTGGAGGGGTTTAACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCaaaa                                                                                  <  2:468126/150‑1 (MQ=255)
                                                                      tGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGGGAAAGGGACCTCACAGTGTAGGTGCGGGTTGACCTGCGAAGAGGCCCACACAGGCCACACAACATCAGGATACCCGCCAAAATcc                                                                               >  2:8666875/1‑150 (MQ=255)
                                                                         tACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCcctc                                                                            <  2:7930625/150‑1 (MQ=255)
                                                                              cttGCGCGACCGTGTCTGACGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAATGGGACCTCACTGTGTTGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTAtt                                                                       <  1:5837112/148‑1 (MQ=255)
                                                                              cATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTAGCGGGTGAAAGGGGCCTCACAGTGTAGGATCGGGGTCACCTGCGAAGAGGCCCACACACGCCACACAACATCCTGATAACCGCCAAAATACCACCTAtt                                                                       >  2:6480095/1‑150 (MQ=255)
                                                                                 gcgcGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATtctc                                                                    <  2:2805105/150‑1 (MQ=255)
                                                                                     gACCGTGTCTTATGGGAATCTGTTTTGCGTGTGCGGGCCTCTTCGGGGTTTAACGGGACCTCACTGTTTTGGAGCGGGTTCACCTGCTAAGTGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATa                                                                <  2:3303984/150‑1 (MQ=255)
                                                                                         gtgtCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCGTCGCGGGTGAAAGGGACCTCACAGGGTAGGGGCGGGTTGACCGGCGAAGAGGCCCCCAAAGGCCACACAACATCCTGGTACCCGCCCAAATACCACCGATTGTGAATAGcgc                                                            >  2:3258261/1‑148 (MQ=255)
                                                                                              tGAAGGGCATCAGTGTTGCCGGGGCGGGCCGCTTCGCGGGTGTAAGGGACCACACACTGTTGGGGGCGGGTGGCCGGCGCACAGGCCCACACACGCCACACAACATGCCGAGAACCCCCAAAATAACACCGATTGTCAATAGCTcccagc                                                       >  2:1487025/1‑147 (MQ=255)
                                                                                                aaGGGCATCTGTGTTGCCTGTGCGGGCCTCTTCGCGGTTGAATGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGc                                                     <  2:8236870/150‑1 (MQ=255)
                                                                                                aaGGGCATCTGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGc                                                     <  2:4367014/150‑1 (MQ=255)
                                                                                                   ggCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTc                                                  <  2:1970539/150‑1 (MQ=255)
                                                                                                    gCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAACGGTACCTCACAGTGTTGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTcc                                                 <  2:5348460/150‑1 (MQ=255)
                                                                                                     catctttgtTGCGTGTGCGGGCCTTTTCGGGGTTTAATGGGACCTCACTGTGTAGGTGGGGTTTCCCCTGCGACGTGGCCCACACAGGAAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCt                                                <  1:5285829/144‑1 (MQ=255)
                                                                                                           tgtTGCCGGTGCGGGCCACATCGCGGGGGAAAAGGACCACACAGGGGAGGAGCGGGGTCACCAGCGCAGAGGCCCACACAGGCAACACAACATCCAGATACCCCCCACAATCCCACCGATTCTCAATAGCTCCCATAGCTTACAGCcccc                                          >  2:363437/1‑150 (MQ=255)
                                                                                                             TTGCCTGTGCGGGCCTCTTCGCGGGTGAACGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAg                                        >  2:6279288‑M1/1‑150 (MQ=255)
                                                                                                              TGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTTTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAgg                                       >  2:5355008‑M1/1‑150 (MQ=255)
                                                                                                               gcgtgtgCGGGCCTCTTCGCGGGTTAATGGGACCTCACTGTGTAGGAGCGGGTTCACCTGCTAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGga                                      <  2:1074749/147‑1 (MQ=255)
                                                                                                               GCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTTGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGga                                      >  2:1339208‑M1/1‑150 (MQ=255)
                                                                                                                CCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGgag                                     <  2:3738105‑M1/150‑1 (MQ=255)
                                                                                                                  tgtgCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccag                                   <  2:7099973‑M2/150‑82 (MQ=255)
                                                                                                                   GTGGGGGCCTCTTGGGGGGTTAAAGGGACCACACTGTTTAGGCGCGGTTTCACGTGCTAAGCGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGaca                                  >  1:6972655‑M1/1‑150 (MQ=255)
                                                                                                                    tgggggCCTCTTCGGGGGTGAACGGGACCTCACTGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGacac                                 <  1:602098/147‑1 (MQ=255)
                                                                                                                    tgggggCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTTGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGacac                                 <  1:4978833/147‑1 (MQ=255)
                                                                                                                    tgCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagca                                 <  2:2971553‑M2/150‑84 (MQ=255)
                                                                                                                    TGCGGGCCTCTTCGGGGGTGAATGGGACCTCACCGTGTAGGAGGGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGacac                                 >  2:3879916‑M1/1‑150 (MQ=255)
                                                                                                                         gCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCa                            <  2:5195077/150‑1 (MQ=255)
                                                                                                                            tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                         <  2:5241968/150‑1 (MQ=255)
                                                                                                                             ctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgc                        <  2:4797109‑M2/150‑93 (MQ=255)
                                                                                                                             ctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgc                        <  2:6930156‑M2/150‑93 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgcc                       <  1:1279994‑M2/150‑94 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                       >  2:5982396/1‑150 (MQ=255)
                                                                                                                              ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                       <  2:7832277/150‑1 (MQ=255)
                                                                                                                                 gcgGTTTAAAGGCACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGCCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                    <  2:6129977/150‑1 (MQ=255)
                                                                                                                                 gcgGGTTAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                    <  2:1222621/150‑1 (MQ=255)
                                                                                                                                 gcgGGTGACAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                    <  2:7205008/150‑1 (MQ=255)
                                                                                                                                 gcgGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgccggg                    <  2:769163‑M2/150‑97 (MQ=255)
                                                                                                                                 gcgGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                    <  2:2050040/150‑1 (MQ=255)
                                                                                                                                   gggTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGt                  <  2:4039737/150‑1 (MQ=255)
                                                                                                                                    ggTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgccggggac                 <  1:3027111‑M2/150‑100 (MQ=255)
                                                                                                                                     gTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTaa                >  1:8204181/1‑150 (MQ=255)
                                                                                                                                      tGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAAc               >  1:4598755/1‑150 (MQ=255)
                                                                                                                                         aaGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCAc            >  1:2459581/1‑150 (MQ=255)
                                                                                                                                          aGGGACCTCACAGTGTAGGCGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACg           >  1:2228189/1‑150 (MQ=255)
                                                                                                                                          aGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACg           >  2:4859585/1‑150 (MQ=255)
                                                                                                                                           gggACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATccc                                                                              >  2:2180717/1‑82 (MQ=255)
                                                                                                                                           gggACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGa          <  2:8417541/150‑1 (MQ=255)
                                                                                                                                            ggACCTCACAGTGTAGGAGGGGGTTCACCTGCGAAGTGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGaa         <  2:7381929/150‑1 (MQ=255)
                                                                                                                                            ggACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGaa         >  1:2258822/1‑150 (MQ=255)
                                                                                                                                                cTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTa                 <  2:2675965/138‑1 (MQ=255)
                                                                                                                                                 tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACg    >  1:4510186/1‑150 (MQ=255)
                                                                                                                                                 tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACg    >  1:964253/1‑150 (MQ=255)
                                                                                                                                                 tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACg    >  1:8621097/1‑150 (MQ=255)
                                                                                                                                                   acaGTGTTGGTGCGGGTTCACCTGCGACGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGc  <  2:1591622/150‑1 (MQ=255)
                                                                                                                                                   acaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGc  >  2:7554803/1‑150 (MQ=255)
                                                                                                                                                   acaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGc  >  1:3712106/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGGCATCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGC  >  NC_002947/4538455‑4538751

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