Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 614,495 A→T 6.6% L252Q (CTG→CAG)  ispA ← farnesyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947614,4950AT 6.6% 475.7 / 15.2 150L252Q (CTG→CAG) ispAfarnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (69/70);  new base T (7/3);  total (76/74)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

ATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATCCAG  >  NC_002947/614348‑614640
                                                                                                                                                   |                                                                                                                                                    
aTAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCAAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAgg                                                                                                                                                    <  1:1017687/150‑1 (MQ=255)
    gTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCg                                                                                                                                                >  1:1094629/1‑150 (MQ=255)
     tGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCgg                                                                                                                                               <  1:37414/150‑1 (MQ=255)
      gAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATGGGGTTTGGCGGCTTGCAGGCCCCGCTGGGCCgcg                                                                                                                                              >  2:1342561/1‑148 (MQ=255)
      gAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCggg                                                                                                                                              <  2:945748/150‑1 (MQ=255)
         gCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAg                                                                                                                                           >  2:2047270/1‑150 (MQ=255)
            gTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTc                                                                                                                                        <  1:1884232/150‑1 (MQ=255)
             tttAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCg                                                                                                                                       >  2:2324056/1‑150 (MQ=255)
                 gTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGttt                                                                                                                                  <  2:1970225/150‑1 (MQ=255)
                           tCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCt                                                                                                                        >  2:1374869/1‑150 (MQ=255)
                             taCAATATAACGCTCCAAGTCCCGCAGCGTCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCCAGTAGGTCGGTTTGTCACACGCTAt                                                                                                                      <  2:1108594/149‑1 (MQ=255)
                             gACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTAt                                                                                                                      >  2:596181/1‑150 (MQ=255)
                             gACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTAt                                                                                                                      >  1:1809042/1‑150 (MQ=255)
                             gACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGACTGTTCGCCGGACCCTTGCTGTGCGGCCAGCGGCAGGTCGCGCAGTTCGGGAACCTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTAGGTTTTTCACGCGCTAt                                                                                                                      >  1:1375695/1‑150 (MQ=255)
                                aaTATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGCAGGTCGGTTTGTCCCGCGCTATGTc                                                                                                                   >  2:1252315/1‑150 (MQ=255)
                                  tataACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTg                                                                                                                <  2:1809042/150‑1 (MQ=255)
                                     aaCGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCAc                                                                                                             >  1:169625/1‑150 (MQ=255)
                                      aCGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCAcc                                                                                                            >  1:1423225/1‑150 (MQ=255)
                                       cgcgCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCt                                                                                                           >  2:1962737/1‑150 (MQ=255)
                                        gcgcCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCCGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCtt                                                                                                          >  2:52152/1‑150 (MQ=255)
                                          gcCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGt                                                                                                                                         <  2:215727/119‑1 (MQ=255)
                                                gCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCATGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTg                                                                                                  >  2:965873/1‑150 (MQ=255)
                                                 cccGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCATGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGc                                                                                                 <  2:1177709/150‑1 (MQ=255)
                                                  ccGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGcc                                                                                                >  2:244788/1‑150 (MQ=255)
                                                  ccGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGcc                                                                                                >  1:980099/1‑150 (MQ=255)
                                                  ccGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGcc                                                                                                >  1:2346578/1‑150 (MQ=255)
                                                  ccGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGcc                                                                                                >  2:1587329/1‑150 (MQ=255)
                                                     cAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCaa                                                                                             <  2:1663989/150‑1 (MQ=255)
                                                      aGCGGCTCGGCCTTTTCGCCCAAACCTTGCAGTGCCGTCAGCGCCTGGTCGCTCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAg                                                                                            <  2:2200977/150‑1 (MQ=255)
                                                           cTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGc                                                                                       <  2:1094629/150‑1 (MQ=255)
                                                              ggCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCGTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGTCCCAGTAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGTCGTTTGCCCAAGGTAGCAgt                                                                                    <  1:1823236/150‑1 (MQ=255)
                                                              ggCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAgt                                                                                    <  1:596181/150‑1 (MQ=255)
                                                                ccTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAg                                                                                            >  1:167803/1‑140 (MQ=255)
                                                                       gCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCAcgc                                                                                                                           >  1:545729/1‑103 (MQ=255)
                                                                       gCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCAcgc                                                                                                                           <  2:545729/103‑1 (MQ=255)
                                                                       gCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAg                                                                                            >  1:1994382/1‑133 (MQ=255)
                                                                       gCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAg                                                                                            <  2:1994382/133‑1 (MQ=255)
                                                                        ccGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTa                                                                                                                                            <  1:1999980/86‑1 (MQ=255)
                                                                          gAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCa                                                                                      >  2:1659740/1‑136 (MQ=255)
                                                                          gAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCa                                                                                      <  1:1659740/136‑1 (MQ=255)
                                                                           aaCCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCCGGGTCCGGTAGGACGTGCTTGAACAGGTAGGTTTGGTACGTTTTCGGTTGCCCAAGGTAGCAGTGTCGCTTTCCACg                                                                       <  1:1252315/150‑1 (MQ=255)
                                                                            aCCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGt                                                                      <  2:2346578/150‑1 (MQ=255)
                                                                             cccTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTc                                                                     >  1:603733/1‑150 (MQ=255)
                                                                              ccTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTa                                                                                         >  1:805468/1‑129 (MQ=255)
                                                                              ccTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTa                                                                                         <  2:805468/129‑1 (MQ=255)
                                                                               cTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGa                                                                   >  1:1605051/1‑150 (MQ=255)
                                                                               cTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGa                                                                   <  2:439903/150‑1 (MQ=255)
                                                                                ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAg                                                                  <  1:1047443/150‑1 (MQ=255)
                                                                                  gCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGa                                                                >  1:1068649/1‑150 (MQ=255)
                                                                                         ggCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATgtcgtc                                                         <  1:1962737/150‑1 (MQ=255)
                                                                                           ccAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCt                                                       >  1:860482/1‑150 (MQ=255)
                                                                                           ccAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCt                                                       <  1:360742/150‑1 (MQ=255)
                                                                                            cAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTc                                                                           >  2:131774/1‑129 (MQ=255)
                                                                                               cgcCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATgtc                                                            >  1:286379/1‑141 (MQ=255)
                                                                                               cgcCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATgtc                                                            <  2:286379/141‑1 (MQ=255)
                                                                                                        gcgcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGc                                          >  1:1579105/1‑150 (MQ=255)
                                                                                                        gcgcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGCCGTTTGCCCATGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTTCTGCACCTGGAATGc                                          <  2:700749/150‑1 (MQ=255)
                                                                                                          gcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATTTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGGTGTCGTCCTGCACCTGGACTGCCa                                       >  2:928000/1‑150 (MQ=255)
                                                                                                          gcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCa                                       >  2:699876/1‑150 (MQ=255)
                                                                                                          gcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCa                                       <  1:52152/150‑1 (MQ=255)
                                                                                                          gcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCa                                       >  1:2102202/1‑150 (MQ=255)
                                                                                                          gcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTACGTCTGCACCTTGGCGTTTGCCCa                                                                                              >  1:2162193/1‑96 (MQ=255)
                                                                                                          gcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTACGTCTGCACCTTGGCGTTTGCCCa                                                                                              <  2:2162193/96‑1 (MQ=255)
                                                                                                           cAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAg                                      <  1:1342561/150‑1 (MQ=255)
                                                                                                              ttCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCACAGGc                                    >  2:1348820/1‑150 (MQ=255)
                                                                                                                 gATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAAt                                >  1:2323806/1‑150 (MQ=255)
                                                                                                                 gATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAAt                                >  1:1450938/1‑150 (MQ=255)
                                                                                                                 gATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCGGGAATGCCAGGCCAAt                                <  1:195250/150‑1 (MQ=255)
                                                                                                                      cATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCCCGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCt                           <  2:341433/150‑1 (MQ=255)
                                                                                                                       aTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCa                                       <  1:785742/137‑1 (MQ=255)
                                                                                                                       aTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCtg                          >  2:1452734/1‑150 (MQ=255)
                                                                                                                       aTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCtg                          <  2:1806257/150‑1 (MQ=255)
                                                                                                                         aGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCACAGGCCATTGGCCtgg                         <  2:603733/150‑2 (MQ=255)
                                                                                                                           gCTTTGGCGGCTTCCAGGCCCAGCTGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCa                      <  1:1393969/150‑1 (MQ=255)
                                                                                                                           gCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCACAGGCCATTGGCCtggca                       <  1:1348820/150‑4 (MQ=255)
                                                                                                                           gCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATgg                              >  2:1424573/1‑142 (MQ=255)
                                                                                                                             tttGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTcc                                                        <  2:1287408/115‑1 (MQ=255)
                                                                                                                             tttGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTcc                                                        >  1:1287408/1‑115 (MQ=255)
                                                                                                                                ggcggcTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCtg                          <  1:2310108/141‑1 (MQ=255)
                                                                                                                                   gtcTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGtt                                                                                                    >  2:2090911/3‑65 (MQ=255)
                                                                                                                                   gtcTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGtt                                                                                                    <  1:2090911/63‑1 (MQ=255)
                                                                                                                                   ggcTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGc                                    <  2:2251515/128‑1 (MQ=255)
                                                                                                                                     ctacccggCCCTGCTGGGCCTGGAAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATaa                   >  2:1052448/7‑143 (MQ=255)
                                                                                                                                     ctacccagCCCTGCTGGGCCTGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAAGCCAAAGGCCTGTGCATAAGTTTGCa            >  2:1705353/8‑150 (MQ=255)
                                                                                                                                     cTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa            <  2:169625/150‑1 (MQ=255)
                                                                                                                                     cTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa            <  2:1605051/150‑1 (MQ=255)
                                                                                                                                      ttCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGAt                                                               >  1:2351791/1‑99 (MQ=255)
                                                                                                                                      ttCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGc                             >  2:971321/1‑132 (MQ=255)
                                                                                                                                      ttCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGc                             <  1:971321/132‑1 (MQ=255)
                                                                                                                                      tacccggCCCTGCTGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAAt                                            <  1:878436/113‑1 (MQ=255)
                                                                                                                                      tacccggCCCTGCTGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAAt                                            >  2:878436/6‑118 (MQ=255)
                                                                                                                                         ccggCCCTGCTGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGc        >  1:310337/3‑150 (MQ=255)
                                                                                                                                         cAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGc                             <  1:1866217/129‑1 (MQ=255)
                                                                                                                                         cAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGc                             >  2:1866217/1‑129 (MQ=255)
                                                                                                                                          aGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAAt                                >  1:950521/1‑125 (MQ=255)
                                                                                                                                            gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCa                                                                                      <  1:2178175/70‑1 (MQ=255)
                                                                                                                                            gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa            >  1:967366/1‑143 (MQ=255)
                                                                                                                                            gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa            <  2:967366/143‑1 (MQ=255)
                                                                                                                                              ccAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAg           <  2:1368765/142‑1 (MQ=255)
                                                                                                                                              ccAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAg           >  1:1368765/1‑142 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATCCAG  >  NC_002947/614348‑614640

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