Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,828,109 A→T 5.9% S14C (AGC→TGC)  yggW → coproporphyrinogen/heterocyclic compound oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,828,1090AT 5.9% 636.3 / 4.2 198S14C (AGC→TGC) yggWcoproporphyrinogen/heterocyclic compound oxidase
Reads supporting (aligned to +/- strand):  ref base A (31/152);  new base T (10/2);  total (44/154)
Fisher's exact test for biased strand distribution p-value = 3.07e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

TGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAGCTG  >  NC_002947/5827960‑5828258
                                                                                                                                                     |                                                                                                                                                       
tGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTa                                                                                                                                                         <  2:5502056/150‑1 (MQ=255)
tGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTa                                                                                                                                                         <  2:8505708/150‑1 (MQ=255)
tGGGGACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGACATGCTGTCGCCCCCCGGCGCGGCAGGCTTTa                                                                                                                                                         <  2:2842087/150‑1 (MQ=255)
 gggTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTag                                                                                                                                                        <  1:4440187/150‑1 (MQ=255)
 gggTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTag                                                                                                                                                        <  1:4124656/150‑1 (MQ=255)
 gggTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTag                                                                                                                                                        <  1:7783919/150‑1 (MQ=255)
   gTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTagca                                                                                                                                                      <  1:7245601/150‑1 (MQ=255)
       cGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCt                                                                                                                                                  <  2:1836177/150‑1 (MQ=255)
        ggAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTg                                                                                                                                                 <  1:6385873/150‑1 (MQ=255)
         gAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGCGTTTGGAGccagg                                                                                                                                                >  2:3199941/1‑147 (MQ=255)
          aGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGccgcgc                                                                                                                                               <  2:2349275/150‑5 (MQ=255)
          aGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTgcc                                                                                                                                               <  2:1536747/150‑1 (MQ=255)
          aGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTgcc                                                                                                                                               <  2:6071798/150‑1 (MQ=255)
              cAACTGCTCCTGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTGCTAAAGCCTgccgcgc                                                                                                                                           <  1:2713641/150‑3 (MQ=255)
                 cTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTg                                                                                                                                        <  2:3472174/150‑1 (MQ=255)
                  tGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTgg                                                                                                                                       <  1:5171167/150‑1 (MQ=255)
                     tCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGCGGTTGTAGACGGCCGCCCCCGgcgc                                                                                                                                    >  1:2312516/1‑150 (MQ=255)
                     tCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGCTTTATCAACCAGCCCCCCCCGgcgc                                                                                                                                    >  1:7195564/1‑150 (MQ=255)
                     tCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGgcgc                                                                                                                                    <  1:259191/150‑1 (MQ=255)
                      ccAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCt                                                                                                                                   <  1:1843918/150‑1 (MQ=255)
                        aGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGCGGTTTCAGCCAGCCCCCCCCGGCGCGGt                                                                                                                                 >  1:6490456/1‑150 (MQ=255)
                        aGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGt                                                                                                                                 <  2:8608679/150‑1 (MQ=255)
                           gcCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACa                                                                                                                              <  1:5991408/150‑1 (MQ=255)
                             cGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATc                                                                                                                            <  1:4775983/150‑1 (MQ=255)
                                      cccAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCg                                                                                                                  <  1:3577180/150‑1 (MQ=255)
                                            agACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGc                                                                                                            <  2:8214190/150‑1 (MQ=255)
                                              aCAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGGTTTTGGAGCCAGCCGCCGCGGGCGCGGTACAACAACAGTACGTGGTGCGt                                                                                                          >  2:6773248/1‑150 (MQ=255)
                                               cAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTa                                                                                                        <  2:7228085/150‑1 (MQ=255)
                                                aaGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTAc                                                                                                       <  2:3436387/150‑1 (MQ=255)
                                                  gAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGGGGCGGTAGGCGCCGGCCGCCCCGGGCGCGGGACAACCACAATCCCGGGTGCggaagc                                                                                                      >  2:1600567/1‑145 (MQ=255)
                                                     ccAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGCGTTTGCAGCCTGCCGCCGCGGGCGCGGTACATACACAATGCGGTGTGCGTAAGCaaa                                                                                                  >  2:2820207/1‑150 (MQ=255)
                                                     ccAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCaaa                                                                                                  <  2:6509383/150‑1 (MQ=255)
                                                      cAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAAt                                                                                                 <  1:1560165/150‑1 (MQ=255)
                                                      cAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAAt                                                                                                 <  1:1560251/150‑1 (MQ=255)
                                                        gCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCGTTATCGGCCCGCCGCCGCGGGCGCGGTTCATACACAATCACGGGTGCGGACGCaaaatc                                                                                               >  2:4225066/1‑147 (MQ=255)
                                                         cTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGcc                                                                                              <  1:4843827/150‑1 (MQ=255)
                                                         cTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGcc                                                                                              <  1:3058882/150‑1 (MQ=255)
                                                              ccACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTAt                                                                                         <  2:7602885/150‑1 (MQ=255)
                                                               cACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGCGTTTGCAGCCAGCCCCCCCGGGCGCGGTAAAAACACAATACGGGGTGCGTACCCCAAAGCCCATaat                                                                                        >  2:2130754/1‑148 (MQ=255)
                                                                aCCGCGCCCGTTCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTg                                                                                       <  2:2350530/150‑1 (MQ=255)
                                                                 ccGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGGGGCGGTGGCAGCCCGCCCCCGCGGGCGCGGGACAACCACAATCCGGGGGGCGGACGCCAAAGccccgaatgg                                                                                      >  2:3743746/1‑143 (MQ=255)
                                                                   gcgcCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCCGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGa                                                                                    <  1:3507084/150‑1 (MQ=255)
                                                                    cgcCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGAc                                                                                   <  1:7511164/150‑1 (MQ=255)
                                                                      cccGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACtt                                                                                 <  2:3592294/150‑1 (MQ=255)
                                                                      cccGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACtt                                                                                 <  2:3866681/150‑1 (MQ=255)
                                                                        cGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCATTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCa                                                                               <  2:2039230/150‑1 (MQ=255)
                                                                          tGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAAc                                                                             <  1:507081/150‑1 (MQ=255)
                                                                           gCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGGCGGGGTTAGCAGCCGGCCCCCGCGGGCGCGGTACAACCACCATACCGGGGGCGGGCGCAAATGCCCCTAATGCGACATCAAca                                                                            >  2:6405060/1‑149 (MQ=255)
                                                                             cATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTcc                                                                          <  1:1245681/150‑1 (MQ=255)
                                                                               tGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCa                                                                        <  1:4213789/150‑1 (MQ=255)
                                                                                ggCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTccccc                                                                       <  2:8703212/150‑3 (MQ=255)
                                                                                ggCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCAc                                                                       <  2:4100157/150‑1 (MQ=255)
                                                                                ggCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCAc                                                                       <  2:864002/150‑1 (MQ=255)
                                                                                ggCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCAc                                                                       <  1:7615777/150‑1 (MQ=255)
                                                                                ggCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCAc                                                                       <  1:6155350/150‑1 (MQ=255)
                                                                                ggCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCAc                                                                       <  2:2683714/150‑1 (MQ=255)
                                                                                     tgctgcGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCGTTAGAAGCCAGCCCCCGCAGGCGCGGTTCAAACACAATACGGTGTGCGTAAGCACATTCCCGAATTTCGACATCAACACCCCCCCAGc                                                                  >  2:7922718/1‑150 (MQ=255)
                                                                                      gctgcGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGcc                                                                 <  1:297054/150‑1 (MQ=255)
                                                                                            cAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCt                                                           <  1:3191526/150‑1 (MQ=255)
                                                                                             aaaCGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCGGGCGTTAGCAGCCTGCCCCCCCGGGCGCTGTACAACCACAAACCCAGGTGCGTACGCCAAAGCCCCGATTGCGACCTCAACGCCCCCCccgccgcggcgg                                                          >  2:6294083/1‑144 (MQ=255)
                                                                                              aaCGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGa                                                         <  1:2650146/150‑1 (MQ=255)
                                                                                                 gTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTATCACCCAGCCCCCCCAGGCGCAGTAAAACCACAATACGTGGTTCCTAACCCAAATCCCATAATTCGACATCAAAAACCACCCCCCCGGGCCAGAAca                                                      >  2:5837225/1‑149 (MQ=255)
                                                                                                 gTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACt                                                      <  1:549905/150‑1 (MQ=255)
                                                                                                  tCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTa                                                     <  1:3977013/150‑1 (MQ=255)
                                                                                                       gCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGa                                                <  1:7074608/150‑1 (MQ=255)
                                                                                                       gCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGa                                                <  2:3507902/150‑1 (MQ=255)
                                                                                                        ccTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGaa                                               <  1:7060413/150‑1 (MQ=255)
                                                                                                          tGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGa                                             <  1:5642379/150‑1 (MQ=255)
                                                                                                           ggCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCATGCAGCCGGGCCTGAACTACCGGAAGAt                                            <  1:2251738/150‑1 (MQ=255)
                                                                                                             cATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGc                                          <  1:7811860/150‑1 (MQ=255)
                                                                                                              aTGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGcc                                         <  1:5289051/150‑1 (MQ=255)
                                                                                                                      aaTGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGAAGCCTGCCGCCGCTGGCGCTGTTCATCCGGATTCCGTGGTGCGTTCGCACATTCCCTTTTTGCGACTTCAACTCCCACGCTGCCGGGCCTGAACTACCGGAAGTTGCCTACGTcac                                 <  2:3524013/150‑1 (MQ=255)
                                                                                                                       aTGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTcaca                                <  1:2803660/150‑1 (MQ=255)
                                                                                                                         gCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCGGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGc                              >  2:493151/1‑150 (MQ=255)
                                                                                                                           tGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCCCGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGccc                            <  1:3327044/150‑1 (MQ=255)
                                                                                                                             ttgCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTg                          <  2:4470565/148‑1 (MQ=255)
                                                                                                                                ccccccGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTc                       <  2:5649042/150‑1 (MQ=255)
                                                                                                                                   cccGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTcgac                    <  2:3657956/150‑1 (MQ=255)
                                                                                                                                     cGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTcgacga                  <  2:7010150/150‑1 (MQ=255)
                                                                                                                                     cGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTcgacga                  <  1:8057779/150‑1 (MQ=255)
                                                                                                                                       gcgcGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGAcc                <  2:106553/150‑1 (MQ=255)
                                                                                                                                           ggCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGa            >  1:8604194/1‑150 (MQ=255)
                                                                                                                                             cAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGAcc          >  2:7481030/1‑150 (MQ=255)
                                                                                                                                             cAGGCTTTAGCAGCCTGCCCCCGCTGGCGCTGTACATCCACATTACGTGGTGCGTACGCAAAAGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCAGCTCGACGACCTCGAcc          >  2:6095727/1‑150 (MQ=255)
                                                                                                                                              aGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCa         <  1:4661403/150‑1 (MQ=255)
                                                                                                                                              aGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCa         >  2:5351683/1‑150 (MQ=255)
                                                                                                                                              aGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCa         <  2:4584598/150‑1 (MQ=255)
                                                                                                                                              aGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCa         >  1:5960893/1‑150 (MQ=255)
                                                                                                                                                gCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAgg       >  1:402242/1‑150 (MQ=255)
                                                                                                                                                gCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAgg       >  2:3837057/1‑150 (MQ=255)
                                                                                                                                                 cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGa      <  1:4049364/150‑1 (MQ=255)
                                                                                                                                                 cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGa      <  2:2678114/150‑1 (MQ=255)
                                                                                                                                                 cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGa      >  2:754962/1‑150 (MQ=255)
                                                                                                                                                  tttAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAg     >  1:8291925/1‑150 (MQ=255)
                                                                                                                                                   ttAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgc    <  2:4007682/150‑1 (MQ=255)
                                                                                                                                                   tgagcagcCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgc    <  1:537084/148‑1 (MQ=255)
                                                                                                                                                    tAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgct   >  1:6650128/1‑150 (MQ=255)
                                                                                                                                                    tAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgct   <  1:1700757/150‑1 (MQ=255)
                                                                                                                                                     agcagcCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAACTGCCCTTATTGCGACTTCAACTCCCACGCCGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgctg  <  2:3053451/150‑1 (MQ=255)
                                                                                                                                                     agcagcCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgctg  >  1:6780918/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGGGTACCGGAGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAGCTG  >  NC_002947/5827960‑5828258

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: