Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,850,713 G→C 5.6% intergenic (‑88/+11) PP_3394 ← / ← PP_3395 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,850,7130GC 5.6% 483.9 / 8.9 150intergenic (‑88/+11)PP_3394/PP_33953‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (84/56);  new base C (3/6);  total (88/62)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

TCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACC  >  NC_002947/3850565‑3850862
                                                                                                                                                    |                                                                                                                                                     
tCACTTGCAGGCTCTGCAAGCCATCGCGCCGGCCGACGTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATAGTCAGc                                                                                                                                                      >  1:541777/1‑150 (MQ=255)
tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                                                      >  1:3925230/1‑150 (MQ=255)
tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                                                      <  1:5044654/150‑1 (MQ=255)
    ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACg                                                                                                                                                  >  2:2614294/1‑150 (MQ=255)
    ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGCTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCAcc                                                                                                                                                  >  2:3822455/1‑149 (MQ=255)
     tGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                                                 <  1:1015212/150‑1 (MQ=255)
     tGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                                                 <  2:444469/150‑1 (MQ=255)
     tGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                                                 <  1:3185706/150‑1 (MQ=255)
     tGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                                                 <  2:1812833/150‑1 (MQ=255)
        aGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGAGGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATgg                                                                                                                                              <  1:2657648/150‑1 (MQ=255)
        aGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATgg                                                                                                                                              <  2:5402714/150‑1 (MQ=255)
        aGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATgg                                                                                                                                              <  1:2683614/150‑1 (MQ=255)
         ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                                             <  1:3545997/150‑1 (MQ=255)
         ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                                             <  1:3402058/150‑1 (MQ=255)
          gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATGGTTATTGTTGTCATCGCGGTCTCTCTGGCTGAAGCGTTGCGCCAGACTCAGAGCTGGATGACGGGCAGACAAGACGTGGTTCGCGGGCTGATCCAGCGTCCGCCCCGATGGag                                                                                                                                            >  2:5452424/1‑149 (MQ=255)
          gCTCTGCAAGCCATCGCGCAGGCCGACTTCATGGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                                                            >  2:2562839/1‑150 (MQ=255)
            tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTGCGCGGGCTGATCCATCGTGGCTGACGATGGATCa                                                                                                                                          <  2:774970/150‑1 (MQ=255)
               gCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAGGGATCAGcc                                                                                                                                       >  2:2532883/1‑150 (MQ=255)
                 aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCcg                                                                                                                                     >  2:92601/1‑150 (MQ=255)
                      aTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGccc                                                                                                                                >  1:933148/1‑150 (MQ=255)
                      aTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGccc                                                                                                                                >  2:681860/1‑150 (MQ=255)
                              ggCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                                                      <  2:1117811/120‑1 (MQ=255)
                              ggCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                                                      >  1:1117811/1‑120 (MQ=255)
                               gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACATGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGGCTGGcc                                                                                                                       >  2:4388990/1‑150 (MQ=255)
                                cctcctccATTTATTGTTATGGTTGTCACGGGGGTCTCTGTGTGGGACGGGTTGCGCCAGACTAATTGCTGTATGAGGGGCACAAAAGACGTGTTGCGCGGGCTCACCCAGGTACCCCAACGAGGCACCCGCCCGCGCCCTGCCAGGCCg                                                                                                                      <  2:1061358/143‑1 (MQ=255)
                                ccGACTACATTGATTGTTATGGTGGTCACCGCGGTCTCTCTGGCTGAAGCGTTGCGCCAGACTAAGCGCTGGAGGACGGGCTGAAAAGACGTGTTCCGCGGGCTGAGCCATCGTCACCCCCGATGGATCCGCCCGCGCCCTGCCAGGCCg                                                                                                                      <  2:4057572/150‑1 (MQ=255)
                                ccGACATCATTGATTGTTATTGTTGTCAACGCGGTCTCTCTGTCTCAAGCGTTGCGCCAGACTAAGAGCTGGATGCCGGGCCGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                                                      <  2:3253773/150‑1 (MQ=255)
                                  gACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                                                    >  1:1320050/1‑150 (MQ=255)
                                   aCTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg                                                                                                                   >  2:4898379/1‑150 (MQ=255)
                                   aCTTCATGGATTGTTATTGTTGTCCTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAACAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg                                                                                                                   <  2:4122674/150‑1 (MQ=255)
                                     ttCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCa                                                                                                                 >  1:4605293/1‑150 (MQ=255)
                                     ttCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCa                                                                                                                 >  2:4413019/1‑150 (MQ=255)
                                     ttCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCa                                                                                                                 >  1:1771263/1‑150 (MQ=255)
                                             ttgttatttttCTAATCGCGGTCTCTCTGTCTCAAGCGTTGCGCAACACTAACAGATGGATGACGGGCACAAAACACGTTTTTCGCGGCCTTATCCAACGAAAGCAACTATGTAACACCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                                                         <  2:449843/150‑1 (MQ=255)
                                              tgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAAc                                                                                                        >  1:5617891/1‑150 (MQ=255)
                                               gttcttgtgGCCATCGCGGACTCCCTGTCTTAAGGGTTGCCCCAGCCTAAGAGGGGGATGCGGGGCACAAAACCGGTGTGTCGCGGGCTGATCCATCTACACCCAGTAGGCATCCCCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACa                                                                                                       <  2:187554/150‑1 (MQ=255)
                                                         tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAg                                                                                             >  1:1647372/1‑150 (MQ=255)
                                                          cATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTTTTTCGCGGGCTGATCCATCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGc                                                                                            <  2:137073/150‑1 (MQ=255)
                                                           aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCa                                                                                           >  2:3503236/1‑150 (MQ=255)
                                                           aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCa                                                                                           >  1:2828439/1‑150 (MQ=255)
                                                             cgcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGGCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGGAGCAACAGGTGGTCAAGCAgg                                                                                         >  1:4909391/1‑150 (MQ=255)
                                                                   ctctctGTCTGACGCGTTGCGCCCGACTAAGAGCTGGTTGACGGGCAGAAAAGACGGGTTGCGCGGGCTTATCCATCGTCCGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccga                                                                                   <  2:3359056/150‑1 (MQ=255)
                                                                     ctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACCGGTGGTCAAGCAGGCgccgagc                                                                                 >  1:4946717/1‑150 (MQ=255)
                                                                      tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagcc                                                                                >  1:3315931/1‑150 (MQ=255)
                                                                      tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagcc                                                                                >  1:2065438/1‑150 (MQ=255)
                                                                      tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCAACCTCACCCAAGAAGGATAAGACCCCCCCCAGCCAGGCCGGTGCACGCAGCAACAGGTTGAAAAGCAGGCgcccagcc                                                                                >  1:3547048/1‑150 (MQ=255)
                                                                            tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                                                          >  1:1819185/1‑150 (MQ=255)
                                                                             gAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCg                                                                         >  2:3774677/1‑150 (MQ=255)
                                                                                  gTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACCCCGAGAGc                                                                    >  1:3358122/1‑150 (MQ=255)
                                                                                   ttGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCg                                                                   >  2:4031467/1‑150 (MQ=255)
                                                                                     gcgcCAGACTACGAGCTGTATCACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGCATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtg                                                                 <  2:4810806/150‑1 (MQ=255)
                                                                                         cAGACTAAGAGCTGTATGACGGGCATAAAAGACGTTTTCCGCGGGCTGAACCATCGTCAGCCACGATGCATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGcc                                                             <  2:1196758/150‑1 (MQ=255)
                                                                                           gaatAAGAGCTGTATGAGGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCAACTATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCgg                                                           <  2:3338362/147‑1 (MQ=255)
                                                                                             cTAAGTGCTGCATGAGGGGCATAAAAGCCGTGTTTCGCGGGCTGTCCCTTCGTCAGCCAGTTTGGACCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                                                         <  2:4714936/150‑1 (MQ=255)
                                                                                              tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccga                                                                                   >  1:4953109/1‑123 (MQ=255)
                                                                                              tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGTGTGTGCCGGTcg                                                        >  1:4961451/1‑150 (MQ=255)
                                                                                              tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGTGTGTGCCGGTcg                                                        >  1:4961459/1‑150 (MQ=255)
                                                                                                agagCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACCGCGAGAGCGTGTGCCGGTcgcg                                                      >  2:5009026/1‑150 (MQ=255)
                                                                                                     tGGCTGAGGGGCTCAAAAGAGGTGTGGCGGGGGCTGACCCATCGTCCGCCCCGCTGTATCCCCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGc                                                 <  1:4413019/150‑1 (MQ=255)
                                                                                                      ggATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCa                                                <  2:5151018/150‑1 (MQ=255)
                                                                                                       gATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCATCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAg                                               >  2:4006316/1‑150 (MQ=255)
                                                                                                         tGACGGGCAGAAAAGCCGTGTTTCGCGGGCTGATCCATCGTCAGCCATGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                                             <  1:2045377/150‑1 (MQ=255)
                                                                                                         ggcggggCATAAAAGACGTGTTTCGGGGGCTGATCCATGGTCCGCCACTTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                                             <  1:3138911/146‑1 (MQ=255)
                                                                                                             gggCAGAAAAGACGTGTTTCGCGGGCTGATACATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACg                                         >  1:49299/1‑150 (MQ=255)
                                                                                                               gCATAAATGAGGTGTTGGGCGGGCTTACCCATCGACCGCCACTATGTAACCCCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGaa                                       <  1:3309216/150‑1 (MQ=255)
                                                                                                               gCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGaa                                       >  2:5075355/1‑150 (MQ=255)
                                                                                                                 acaaaaCACGTGTTGCGCGGGCTCATCAATCTACAGCCACTTTGCAACACCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt                                     <  1:2613873/148‑1 (MQ=255)
                                                                                                                  gAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTc                                    >  2:1850164/1‑150 (MQ=255)
                                                                                                                     acgaggtgtTTCGGGGGCTGATCCATCGTCAGCCTCGTTGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAcc                                 <  2:1206911/145‑1 (MQ=255)
                                                                                                                       gACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCAGCCAGGCCGGTGCACCCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCaa                               >  1:1637175/1‑150 (MQ=255)
                                                                                                                        aCGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAg                              <  2:3925230/150‑1 (MQ=255)
                                                                                                                            gtTTCGGGGGCTTATCAATGGTCAGCAAGTTGGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGcc                          <  1:5014586/150‑1 (MQ=255)
                                                                                                                            gtTTCGCGGGCTGATCCTTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGcc                          <  2:1230087/150‑1 (MQ=255)
                                                                                                                             tttCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGccc                         >  2:3599327/1‑150 (MQ=255)
                                                                                                                                 gcgGGCTGATCCATCGACAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAgtc                     <  2:883240/150‑1 (MQ=255)
                                                                                                                                       tGATCCATCGTGGCTGACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtg                                                                                                   >  1:734163/1‑66 (MQ=18)
                                                                                                                                       tGATCCATCGTGGCTGACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtg                                                                                                   <  2:734163/66‑1 (MQ=18)
                                                                                                                                       tGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagt               >  1:2557061/1‑150 (MQ=255)
                                                                                                                                       tGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagt               >  1:2557177/1‑150 (MQ=255)
                                                                                                                                         aTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagtca             >  2:1750883/1‑150 (MQ=255)
                                                                                                                                         aTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagtca             >  1:674028/1‑150 (MQ=255)
                                                                                                                                         aTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagtca             >  1:5116503/1‑150 (MQ=255)
                                                                                                                                         aTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagtca             >  1:809119/1‑150 (MQ=255)
                                                                                                                                          tCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACGGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGtcagtcag            >  2:5137306/1‑150 (MQ=255)
                                                                                                                                          tCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGtcagtcag            >  1:4836709/1‑150 (MQ=255)
                                                                                                                                           ccATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgc           <  2:5517079/150‑1 (MQ=255)
                                                                                                                                            cATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcg          >  1:1243084/1‑150 (MQ=255)
                                                                                                                                            cATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcg          >  1:449305/1‑150 (MQ=255)
                                                                                                                                             aTCGTCAGCCCCGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         <  2:3320859/150‑1 (MQ=255)
                                                                                                                                             aTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         >  2:4710216/1‑150 (MQ=255)
                                                                                                                                             aTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         >  2:3067086/1‑150 (MQ=255)
                                                                                                                                             aTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         >  2:2639758/1‑150 (MQ=255)
                                                                                                                                             aTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         >  2:5562609/1‑150 (MQ=255)
                                                                                                                                             aTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         >  2:250607/1‑150 (MQ=255)
                                                                                                                                             aTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         >  2:1337625/1‑150 (MQ=255)
                                                                                                                                             aTCGCCCGCCACGTGGCAACACCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcgc         <  1:3822455/150‑1 (MQ=255)
                                                                                                                                              tCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCa        >  2:2600754/1‑150 (MQ=255)
                                                                                                                                                gTCAGCCACGATGGATCATCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAAt      <  1:954702/150‑1 (MQ=255)
                                                                                                                                                   aGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAc   <  2:4381876/150‑1 (MQ=255)
                                                                                                                                                    gCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAcc  <  2:2951789/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACC  >  NC_002947/3850565‑3850862

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: