Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,686,549 A→C 13.6% Y26S (TAC→TCC)  recJ → ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,686,5490AC13.6% 780.3 / 46.3 314Y26S (TAC→TCC) recJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base A (163/102);  new base C (2/47);  total (165/149)
Fisher's exact test for biased strand distribution p-value = 3.39e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

TCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCG  >  NC_002947/1686400‑1686695
                                                                                                                                                     |                                                                                                                                                  
tCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTa                                                                                                                                                    <  2:4324834/150‑1 (MQ=255)
      ctgacGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTAcgccgc                                                                                                                                              <  2:3063901/150‑1 (MQ=255)
      ctgacGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTAcgccgc                                                                                                                                              <  1:2047137/150‑1 (MQ=255)
       tgacGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCcgcg                                                                                                                                             >  2:5225594/1‑150 (MQ=255)
             ccGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCgcggc                                                                                                                                       <  2:4597797/150‑1 (MQ=255)
              cGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCgcggcg                                                                                                                                      >  1:6282552/1‑150 (MQ=255)
               gTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGt                                                                                                                                     <  2:4239214/150‑1 (MQ=255)
               gTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGt                                                                                                                                     <  1:104010/150‑1 (MQ=255)
                 ggCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGtaa                                                                                                                                   <  1:5843781/150‑3 (MQ=255)
                 ggCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGc                                                                                                                                   >  1:1421200/1‑150 (MQ=255)
                  gCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGtaaa                                                                                                                                  <  1:2843969/150‑4 (MQ=255)
                  gCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCAAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGtaaa                                                                                                                                  <  1:5101879/150‑4 (MQ=255)
                    cAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGtaaagg                                                                                                                                <  1:2711114/150‑6 (MQ=255)
                    cAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGtaaagg                                                                                                                                <  1:4900965/150‑6 (MQ=255)
                    cAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCaga                                                                                                                                >  1:2476592/1‑150 (MQ=255)
                       cATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCg                                                                                                                             <  2:6393689/150‑1 (MQ=255)
                       cATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTGTACGCCGCGCGCGGCGTGCAGAGCg                                                                                                                             >  2:4007019/1‑150 (MQ=255)
                         tGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAg                                                                                                                           <  2:3403569/150‑1 (MQ=255)
                            aGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGcc                                                                                                                        >  1:4197134/1‑150 (MQ=255)
                                 aaCTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCt                                                                                                                   <  2:6229936/150‑1 (MQ=255)
                                 aaCTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCt                                                                                                                   >  1:3542317/1‑150 (MQ=255)
                                  aCTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTg                                                                                                                  <  1:4627808/150‑1 (MQ=255)
                                    tggaGTGGCTGCTCGGTGAGCGTGTTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCTTTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGa                                                                                                                <  1:3469431/147‑1 (MQ=255)
                                    tGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGa                                                                                                                <  2:4103633/150‑1 (MQ=255)
                                    tGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGa                                                                                                                >  2:2605842/1‑150 (MQ=255)
                                    tGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGa                                                                                                                <  2:783730/150‑1 (MQ=255)
                                    tGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGa                                                                                                                >  1:4794341/1‑150 (MQ=255)
                                    tGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCGGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGa                                                                                                                <  1:5495569/150‑1 (MQ=255)
                                      ctgTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCAATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACa                                                                                                              <  1:3128626/148‑1 (MQ=255)
                                           gctgctCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGc                                                                                                         <  2:2282330/150‑1 (MQ=255)
                                              gctCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTg                                                                                                      <  1:5586531/150‑1 (MQ=255)
                                              gctCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTg                                                                                                      <  2:5610552/150‑1 (MQ=255)
                                               ctCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTgg                                                                                                     >  2:2465214/1‑150 (MQ=255)
                                                tCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGc                                                                                                    <  2:1698044/150‑1 (MQ=255)
                                                  ggTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCAc                                                                                                  >  2:1652941/1‑150 (MQ=255)
                                                      aGCTTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGTTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCct                                                                                              <  1:4954167/150‑1 (MQ=255)
                                                      aGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCAGGCACGCct                                                                                              >  1:4155655/1‑150 (MQ=255)
                                                           ggTTGAAGCAACACATGCGTATCGAACCTCGCCCGGTGCCGCCGACCCTGCCATTTCTGGGTAACCGGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCtgctgc                                                                                         <  1:5384564/150‑1 (MQ=255)
                                                            gTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCTCGcccagcagg                                                                                        >  2:5945318/1‑143 (MQ=255)
                                                            gTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGcc                                                                                        >  2:2557395/1‑150 (MQ=255)
                                                              tGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCt                                                                                      >  1:98561/1‑150 (MQ=255)
                                                                aaGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTAc                                                                                    <  2:3846427/150‑1 (MQ=255)
                                                                  gCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTgg                                                                                                     >  1:6090888/1‑131 (MQ=255)
                                                                   cAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTAcccg                                                                                 >  1:1802502/1‑148 (MQ=255)
                                                                   cAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACcag                                                                                 <  1:2336279/150‑1 (MQ=255)
                                                                     acacATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACcagca                                                                               <  1:3672452/150‑1 (MQ=255)
                                                                        caTGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTGGCTGACCCGCCTTTCCGCCGCGCGCGGCGTGCAGCGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCt                                                                            <  1:4370640/150‑1 (MQ=255)
                                                                         aTGGGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCGTGGCTGACCCGCCTTTCCGCCGCGCGCGGCGTGCAGAGCGCGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTc                                                                           <  1:2128332/150‑1 (MQ=255)
                                                                          tGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCa                                                                          >  1:2559039/1‑150 (MQ=255)
                                                                           gCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCTTTTCTGGTTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCaa                                                                         <  1:6704302/150‑1 (MQ=255)
                                                                           gCGTATCGAACCTCGCCCGCGGCCGCCGACCCTGCCTTTTCTGGGTACCCTGCCACCCTTGCTGACCCGCCTTTACGCGGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCaa                                                                         <  1:519027/150‑1 (MQ=255)
                                                                            cGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAg                                                                        <  2:1245088/150‑1 (MQ=255)
                                                                             gTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCTTTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAgg                                                                       <  1:5689180/150‑1 (MQ=255)
                                                                                tCGACCCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCa                                                                    <  1:869/150‑1 (MQ=255)
                                                                                 cGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCAt                                                                   >  2:5195413/1‑150 (MQ=255)
                                                                                   aaCCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGGAACCTGCCACCCGTGCTGACCCGCCTTTCCGCCGCGCGCGGCGGGCAGAGCGAGGCCGAGCTGGACCCGAGCCTGGCCCGCCTGCTGCCCGACCAGCAGCGCAAGGGCAGCg                                                                 >  2:1079187/1‑150 (MQ=255)
                                                                                   aaCCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGAAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCg                                                                 >  1:1378293/1‑150 (MQ=255)
                                                                                     ccTCGCCCGCGGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTGCCGCCGCGCGCGGCGTGCCGAGCGCGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAg                                                               <  2:1524769/150‑1 (MQ=255)
                                                                                      cTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAgg                                                              >  2:2379730/1‑150 (MQ=255)
                                                                                      cTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAgg                                                              >  2:1034166/1‑150 (MQ=255)
                                                                                         gCCCGCGGCCGCCGCCCCGGCCATGTCTGGGTAACCGGCCCCCCTTGCGGCCCCGCCTGGCCGCCGCGCGCGGCGTGCCGCGCGCGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCGGCCCGACCAGCAGCTCAAGGGCATCGAGGCCg                                                           <  2:5293126/150‑1 (MQ=255)
                                                                                         gCCCGCGGCCGCCGACCCGGCCTTGACTGGGTACCCTGCCACCCTTGCTGCCCCGCCCTTCCCCCGCGCGCGGCGTGCCGCGCGAGGCCGAGCTGGACACGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCg                                                           <  2:3881013/150‑1 (MQ=255)
                                                                                           ccGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGgt                                                         <  2:3919707/150‑1 (MQ=255)
                                                                                             gcggccgccCACCCGGCCATTTCTGGGTACCCTGCCACCCTTGCTGACCCGCCTTGCCGCCGCGCGCGGCGTGCAGCGCGCGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGgtgt                                                       <  1:6445160/147‑1 (MQ=255)
                                                                                               tGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTg                                                                                                      >  1:529392/1‑101 (MQ=255)
                                                                                                gccgccGACCCGGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGa                                                    <  2:886270/150‑1 (MQ=255)
                                                                                                 ccgccgACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGAc                                                   <  2:701799/150‑1 (MQ=255)
                                                                                                 ccgccgACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGACGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCCTCGAGGCCGGTGTGGAc                                                   >  1:291491/1‑150 (MQ=255)
                                                                                                   gccgACCCTGCCATGTCGGGGTACCCGGCCACCCGTGCGGACCCGCCTGGCCGCCGCGCGCGGCGTGCAGCGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACct                                                 <  1:728055/150‑1 (MQ=255)
                                                                                                    ccgACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTCCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACcagca                                                                               <  1:924076/119‑1 (MQ=255)
                                                                                                     cgACCCTGCCATTTCTGGGTAACCTGCCACCCGTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGCGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACctgc                                               <  2:3290573/150‑1 (MQ=255)
                                                                                                      caCCCTGCCATTTGTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGCGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACctgct                                              <  2:3009703/149‑1 (MQ=255)
                                                                                                       aCCCTGCCTTTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTCCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACctgctg                                             <  1:6565268/150‑1 (MQ=255)
                                                                                                         ccTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTggg                                           >  2:495980/1‑149 (MQ=255)
                                                                                                         ccTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGt                                           >  2:5486861/1‑150 (MQ=255)
                                                                                                         ccTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGt                                           >  1:1937824/1‑150 (MQ=255)
                                                                                                          cTGCCATTTCTGGTTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGGGGCGTGCCGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTc                                          <  1:1271905/150‑1 (MQ=255)
                                                                                                           tGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCg                                         >  2:2099248/1‑150 (MQ=255)
                                                                                                             ccATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTCACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGat                                       <  1:791129/150‑2 (MQ=255)
                                                                                                              cATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAg                                      >  2:6784308/1‑150 (MQ=255)
                                                                                                              cATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAg                                      >  2:6503748/1‑150 (MQ=255)
                                                                                                                     gggTACCCTGCCACCCTTGCCGCCCCGACTGCCCGCCGCGCGCGGCGTGCCGAGCGCGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGa                               <  2:5098534/150‑1 (MQ=255)
                                                                                                                      ggTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAc                              >  2:4109028/1‑150 (MQ=255)
                                                                                                                         aaCCTGCCCCCCTGGCTGACCCGCCTTTCCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccag                           <  1:2157313/150‑1 (MQ=255)
                                                                                                                           ccTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcg                         >  1:6003138/1‑150 (MQ=255)
                                                                                                                             tGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcc                       >  2:6728738/1‑150 (MQ=255)
                                                                                                                             tGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcc                       >  1:2216272/1‑150 (MQ=255)
                                                                                                                             ggCCACCCTTGCGGCCCCGCCTTCACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcc                       <  1:1164612/149‑1 (MQ=255)
                                                                                                                                 aCCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgccagcg                   >  1:6370734/1‑150 (MQ=255)
                                                                                                                                   ccTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCa                 >  2:615261/1‑150 (MQ=255)
                                                                                                                                     ttGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATc               <  2:2185936/150‑1 (MQ=255)
                                                                                                                                              cGCCTTTCCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgtc      <  2:5828008/150‑1 (MQ=255)
                                                                                                                                              cGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcca                      >  2:642855/1‑134 (MQ=255)
                                                                                                                                                tctgcaCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg    >  2:602187/6‑150 (MQ=255)
                                                                                                                                                tctgcaCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg    >  1:4015542/6‑150 (MQ=255)
                                                                                                                                                 cTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGc   >  1:3585752/1‑150 (MQ=255)
                                                                                                                                                  tttACGCCGCGCGTGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCg  >  2:3345229/1‑150 (MQ=255)
                                                                                                                                                  tttACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCg  >  1:2374989/1‑150 (MQ=255)
                                                                                                                                                  tttACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCg  >  1:3298668/1‑150 (MQ=255)
                                                                                                                                                  tgcaCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCg  >  2:4159338/4‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCGTGACTGACGACCGTGGCCAGCATGAAGTGCAACTGCAGTGGCTGCTCGGTGAGCGTGGTTGAAGCAACACATGCGTATCGAACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCG  >  NC_002947/1686400‑1686695

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: