Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,538,609 A→T 8.8% intergenic (+58/‑116) 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,6090AT 8.8% 650.1 / 10.7 233intergenic (+58/‑116)PP_4020/cpooxidoreductase/non‑heme chloroperoxidase
Reads supporting (aligned to +/- strand):  ref base A (85/121);  new base T (4/20);  total (91/142)
Fisher's exact test for biased strand distribution p-value = 2.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

TCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGCGTACAGA  >  NC_002947/4538460‑4538758
                                                                                                                                                     |                                                                                                                                                     
tCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTa                                                                                                                                                       <  1:8605014/150‑1 (MQ=255)
 cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                                      <  1:2524071/150‑1 (MQ=255)
 cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                                      <  2:2638135/150‑1 (MQ=255)
 cGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCATTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAg                                                                                                                                                      <  1:2003246/150‑1 (MQ=255)
  gACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGGCGGTGCGGGCCGCGTAGCGGGTGGAAGGGGCCGCACAGTGTAgg                                                                                                                                                     >  2:7426238/1‑150 (MQ=255)
       cGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCgg                                                                                                                                                <  1:4770564/150‑1 (MQ=255)
       cGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGCGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCgg                                                                                                                                                <  2:1476091/150‑1 (MQ=255)
        ggCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGGGAAAGGGCCCGCACAGTGTAGGAGCggg                                                                                                                                               >  2:5746113/1‑150 (MQ=255)
           cGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTc                                                                                                                                            <  2:5103281/150‑1 (MQ=255)
            gCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGCCGCGTCGCGGGGGGAAGGGCCCGCACCGGGGTGGAGCGGGGTCa                                                                                                                                           >  2:6092374/1‑150 (MQ=255)
             cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGGGGGCCGCTTAGCGGGGGAAAGGGACCACAAAATGTTGGAGCGGGTTCAc                                                                                                                                          >  1:7188315/1‑150 (MQ=255)
             cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCAc                                                                                                                                          <  1:8378287/150‑1 (MQ=255)
             cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCAc                                                                                                                                          <  2:4151926/150‑1 (MQ=255)
             cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCAc                                                                                                                                          <  1:4623285/150‑1 (MQ=255)
             cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCATCGCGGGGGGAAGGGACCGCACCCTGTAGGAGCGGGTTGAc                                                                                                                                          >  2:8554610/1‑150 (MQ=255)
             cTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGTGCGGGCCACATCGCGGGGGAAAAGGACCACACACTGTAGGAGCGGGTTGAc                                                                                                                                          >  1:7939560/1‑150 (MQ=255)
               tCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCt                                                                                                                                        <  2:2979932/150‑1 (MQ=255)
                 atagctcatggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGc                                                                                                                                      >  2:5611910‑M1/85‑150 (MQ=255)
                    ggcggcGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGaa                                                                                                                                   <  2:2926185/150‑1 (MQ=255)
                      tcatggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAaga                                                                                                                                 <  2:3949790‑M1/71‑1 (MQ=255)
                      cggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGTGAaga                                                                                                                                 <  1:4556290/150‑1 (MQ=255)
                      cggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAaga                                                                                                                                 <  1:3754604/150‑1 (MQ=255)
                       ggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                                <  1:396407/150‑1 (MQ=255)
                       ggcgACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                                <  2:7865687/150‑1 (MQ=255)
                       catggtcgtgtctcctgggtgcaggaagcgatgggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                                >  2:7217425‑M1/79‑150 (MQ=255)
                       catggtcgtgtctcctgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                                >  1:6760541‑M1/79‑150 (MQ=255)
                                ctggctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcaca                                                                                                                       >  2:3778371‑M2/1‑146 (MQ=255)
                                cgggctgcgCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcaca                                                                                                                       >  1:5940326‑M2/3‑146 (MQ=255)
                                  ggctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacagg                                                                                                                     >  1:1171411‑M2/1‑144 (MQ=255)
                                   gctgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggc                                                                                                                    <  1:7708403‑M2/150‑8 (MQ=255)
                                     tgCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaa                                                                                                                  >  2:873793‑M2/1‑141 (MQ=255)
                                      tgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                                >  1:1016977‑M1/64‑135 (MQ=255)
                                      tgggtgcaggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAagag                                                                                                                                <  2:1016977‑M1/72‑1 (MQ=255)
                                      gcgcCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaac                                                                                                                 >  1:2830809‑M2/1‑140 (MQ=255)
                                       gggtgcaggaagcgaagggagctattgagaataggagggattttggcgggcatcaggaagttGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTCGGTGCGGTTTCACCTGCGAAGAGGCCCACACAGGCAaca                                                                                                                >  2:443701‑M1/63‑150 (MQ=255)
                                       cgcCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaaca                                                                                                                >  2:1426648‑M2/1‑139 (MQ=255)
                                          cTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCATTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACaa                                                                                                             >  1:3293364/1‑150 (MQ=255)
                                          cTTTGGAAGCGGCGCCAAAGGCCAGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTCTTGCCTGTGCGGGCCTCTTCGCGGGTGAACGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACaa                                                                                                             <  2:6549522/150‑1 (MQ=255)
                                             aggaagcgaagggagctattgagaataggagggattttggcgggtatcaggaagttGCGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAgg                                                                                                                               >  1:2732896‑M1/57‑129 (MQ=255)
                                              ggACGCTGCGCAAACGGCCTGGTACCGATGGCGCGACGGTTTCTGACGGGAATCTGTGTTGCCTGGGGGGGCCTCTTCGGGGGTGAACGGGACCTCACTGTTTTGGCGGGGGTTCACCTGCGAAGCGGCCCACACAGGCAACACAACTTc                                                                                                         <  2:3299552/150‑1 (MQ=255)
                                              ggAAGCTGCGGCAAAGGCCGGTCACCGCTTGCGCGACCGTTTCTGATGGGCTTCTGTGTTGCCTGTGCGGGCCTCTTCGGGGGTGAAAGGCACATCACTGTGTTGTAGCGGTTACACCTGCGAAGAGGCCCACACAGGCAACACAACTTc                                                                                                         <  2:5471974/150‑1 (MQ=255)
                                               gacgCTGCGGCAACGCCCTGGTCCCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAACGGGACCTAACAGTGTTGGCGCGGGTTCACCTGCGAAGAGGCCCCCACAGGCAACACAACTTcc                                                                                                        <  2:3198638/147‑1 (MQ=255)
                                                  gcggCGGCAACGGCCGGGTACCGCAGGCGCGACGGTGTCTGACGGGCATCAGTGTTGCCTGTGGGGGCCTCTTCGGGGGTGAAAGGTACCTCACTGTGTGGGCGCGGGTTCACGTGCGATGTGGCCCACACTGGCAACACAACTTCCTGa                                                                                                     <  2:2327298/147‑1 (MQ=255)
                                                         cAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTTCGGGCCTCTTAGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTGACCAGCGAAGAGGCCCACACAGGCAACACAAATTCCAGaaaacccc                                                                                              >  2:4934799/1‑143 (MQ=255)
                                                          aaaGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCGGGGCGGGCCACATCCCCGGGGGAAGGGCCCCCACCGGGGCGGGGCGGGGTCACCAGCGCAGAGGCCCCCACCCGCAACACCACAACCAGAAAcccccc                                                                                             >  2:2664286/1‑150 (MQ=255)
                                                               ccTGGGACCGCAGGCGCGCCCGGGGCGGCCGGGCCGCCGGGTGGCCGGTGGGGGCCTGTTCGGGGGTGAACGGTCCGTCACAGTGTTGGTGCGGTTTCACCTGCGAAGAGGCCCACACTGGCAACACAACTTCCTGATACCCGCCAAAAt                                                                                        <  2:353815/150‑1 (MQ=255)
                                                                 tGGTACCGCTTGCGCGACCGTGTCTGATGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGGGGGTGAAAGGGACCTCACAGTGTTGGAGCGGTTTCACCTGCGAAGTGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATcc                                                                                      <  2:4831131/150‑1 (MQ=255)
                                                                   gTACCGCATGCGCGACCGTGTCTGACGGGCATCTGTGTTGCCTGTGCGGGCCTCTTCGGGGGTGAATGGGACCTCACTGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCcct                                                                                    <  2:8308957/150‑1 (MQ=255)
                                                                   gTACCGAATGCGCGACGGTGTCTGATGGGCATCTGTGTTGCCTGGGCGGGCATCTTGGCGGGTGAACGGGACCTCACTGTGTAGGCGCGGGTTAACCTGCGAAGTGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCcct                                                                                    <  1:6508019/150‑1 (MQ=255)
                                                                    tACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCcctc                                                                                   <  2:6824726/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)
                                                                                 aCCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAg                                                                      <  2:603818/150‑1 (MQ=255)
                                                                                         tGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCATAGCGGGTGAAAGGGCCCGCCCAGTGTAGGAGCGGGGTCGCCGGCGCAGAGGCCCCCACAGGCCACACCACCTCCAGAAACCCGCCCAAATCCCACCGAATCTCAATATCGCCCGTc                                                              >  2:8475595/1‑150 (MQ=255)
                                                                                           aaGGGCATCTGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGc                                                            <  2:4367014/150‑1 (MQ=255)
                                                                                              gtcATCTGTGTTGCCTGTGGGGGCCTCTTGGCGGGTTAAAGGGACCTCACTGTTTTGCTGCGGTTCCACCTGCTATGAGGCCAACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTc                                                         <  2:3884468/148‑1 (MQ=255)
                                                                                               gCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAACGGTACCTCACAGTGTTGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTcc                                                        <  2:5348460/150‑1 (MQ=255)
                                                                                                cATCTGTGTTGCGTGGGCGGGCCTCTGCGGGGGTTAATGGGACCTCACTGTGTTGGTGCGGGTTCACGGGCGATGCGCCCCAAACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCt                                                       <  1:6357715/150‑1 (MQ=255)
                                                                                                        tgGCCTGTGCGGGCCTCTTCGGGGGTGAATGGGACCCCACTGTGTAGGTGGGGGTTCACCTGCTAAGTGGCCCACACAGGAAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAg                                               >  1:8567045‑M1/3‑150 (MQ=255)
                                                                                                        ggGCCTGTGGGGGCCTGTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACGTGCGACGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAg                                               >  2:5284682‑M1/3‑150 (MQ=255)
                                                                                                        TTGCCTGTGCGGGCCTCTTCGCGGGTGAACGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAg                                               >  2:6279288‑M1/1‑150 (MQ=255)
                                                                                                        TTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAg                                               >  2:6643417‑M1/1‑150 (MQ=255)
                                                                                                        TTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAg                                               >  2:5032457‑M1/1‑150 (MQ=255)
                                                                                                          GCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAAGCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGga                                             <  1:959904‑M1/150‑1 (MQ=255)
                                                                                                           CCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGgag                                            <  2:3738004‑M1/150‑1 (MQ=255)
                                                                                                             tgtgCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccag                                          <  2:7099973‑M2/150‑82 (MQ=255)
                                                                                                              gtgGGGGCCTTTTCGGGGTTTAATGGTACAACAATGTGTTGGCGGGGTTTAACCTGCTAAGAGGCCAACACTGGAAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGaca                                         <  1:6972998/150‑1 (MQ=255)
                                                                                                               tgCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagca                                        <  2:2712448‑M2/150‑84 (MQ=255)
                                                                                                               gggggggCTCTTGGGGGGGGAATGGGCCCTCACTGTGTTGGAGGGGGTTCCCCTGCGATGCGGCCCACACCGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGacac                                        <  1:4978922/147‑1 (MQ=255)
                                                                                                                    gCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcg                                   <  2:6941116‑M2/150‑89 (MQ=255)
                                                                                                                    gCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCa                                   <  2:5195077/150‑1 (MQ=255)
                                                                                                                     ccTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCAt                                  >  2:4715760/1‑150 (MQ=255)
                                                                                                                       tctTCGGGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                                <  2:8154704/150‑1 (MQ=255)
                                                                                                                       tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCGcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccg                                <  2:6331446‑M2/150‑92 (MQ=255)
                                                                                                                       tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                                >  2:8708321/1‑150 (MQ=255)
                                                                                                                       tctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGa                                <  2:1324509/150‑1 (MQ=255)
                                                                                                                        ctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAg                               <  2:3641213/150‑1 (MQ=255)
                                                                                                                        ctTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAg                               >  2:407578/1‑150 (MQ=255)
                                                                                                                         ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgcc                              <  1:1279994‑M2/150‑94 (MQ=255)
                                                                                                                         ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                              >  1:2367021/1‑150 (MQ=255)
                                                                                                                         ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGc                              <  2:7832277/150‑1 (MQ=255)
                                                                                                                         ttCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTAACCTGCGAAGAGGCCCAcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgaagacagcaggtcgccgcc                              <  1:6922617‑M2/150‑94 (MQ=255)
                                                                                                                            gcgGTTTAAAGGCACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGCCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                           <  2:6129977/150‑1 (MQ=255)
                                                                                                                            gcgGGTTAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                           <  2:1222621/150‑1 (MQ=255)
                                                                                                                            gcgGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTAt                           <  2:2050040/150‑1 (MQ=255)
                                                                                                                              gggTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCGcacaggcaacactgatgcccttcagacacggtcgcgcatgcggtaccaggcctttgccgcagcttccaaaggcgcagccagcaggtcgccgccgggga                         <  2:3331016‑M2/150‑99 (MQ=255)
                                                                                                                                    aaGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCAc                   >  1:2459581/1‑150 (MQ=255)
                                                                                                                                      gggACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATccc                                                                                     >  2:2180717/1‑82 (MQ=255)
                                                                                                                                       ggACCTCACAGTGTAGGAGGGGGTTCACCTGCGAAGTGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGaa                <  2:7381929/150‑1 (MQ=255)
                                                                                                                                       ggACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGaa                >  1:2258657/1‑150 (MQ=255)
                                                                                                                                       ggACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGaa                <  2:1831111/150‑1 (MQ=255)
                                                                                                                                           cTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTa                        <  2:2675965/138‑1 (MQ=255)
                                                                                                                                            tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACg           >  1:2221530/1‑150 (MQ=255)
                                                                                                                                            tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACg           >  1:6122969/1‑150 (MQ=255)
                                                                                                                                            tCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACg           >  1:4510186/1‑150 (MQ=255)
                                                                                                                                              acaGTGTAGGAGCGGGTTCACCTGCTAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGc         <  2:1591523/150‑1 (MQ=255)
                                                                                                                                              acaGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGc         >  2:7554803/1‑150 (MQ=255)
                                                                                                                                                 gtgtAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATGCTCACTAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGCGTa      <  2:7607334/150‑1 (MQ=255)
                                                                                                                                                    tAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGAAACCACGAAGCACGGCGTACAg   >  1:8707802/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGACAACGGCCGCTTCCCCGGCGGCGACCTGCTGGCTGCGCCTTTGGAAGCTGCGGCAAAGGCCTGGTACCGCATGCGCGACCGTGTCTGAAGGGCATCAGTGTTGCCTGTGCGGGCCTCTTCGCGGGTGAAAGGGACCTCACAGTGTAGGAGCGGGTTCACCTGCGAAGAGGCCCACACAGGCAACACAACTTCCTGATACCCGCCAAAATCCCTCCTATTCTCAATAGCTCCCTTCGCTTCCTGCACCCAGGAGACACGACCATGAGCTATGTAACCACGAAGGACGGCGTACAGA  >  NC_002947/4538460‑4538758

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: