Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,633,966 C→T 5.1% intergenic (+18/‑146) lepB → / → rnc signal peptidase I/ribonuclease 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,633,9660CT 5.1% 509.4 / 11.1 156intergenic (+18/‑146)lepB/rncsignal peptidase I/ribonuclease 3
Reads supporting (aligned to +/- strand):  ref base C (71/77);  new base T (3/5);  total (74/82)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 175 total aligned reads displayed.

ATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGA  >  NC_002947/1633817‑1634115
                                                                                                                                                     |                                                                                                                                                     
aTTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTcc                                                                                                                                                       >  1:3028546/1‑150 (MQ=255)
aTTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATAGTCCGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTAAGCCACTAGCCCAACCTGTCGCGGGTAGGCCTGATCCAATGATACCCATCGGCGCTGTcc                                                                                                                                                       >  2:4360587/1‑150 (MQ=255)
 ttCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCca                                                                                                                                                      >  2:198758/1‑150 (MQ=255)
    ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                   >  1:507990/1‑150 (MQ=255)
    ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                   >  1:3804068/1‑150 (MQ=255)
     cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                  >  2:2178599/1‑150 (MQ=255)
     cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                  >  1:814489/1‑150 (MQ=255)
     cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                  >  2:2270780/1‑150 (MQ=255)
      aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGa                                                                                                                                                 >  2:4237468/1‑150 (MQ=255)
      aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGa                                                                                                                                                 <  2:2171328/150‑1 (MQ=255)
          aaCTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAgc                                                                                                                                             <  1:3237302/150‑1 (MQ=255)
           aCTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAgcg                                                                                                                                            <  1:952008/150‑1 (MQ=255)
             tGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                          >  1:2262951/1‑150 (MQ=255)
             tGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGAGCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                          >  1:2024903/1‑150 (MQ=255)
               cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGa                                                                                                                                        <  1:1654502/150‑1 (MQ=255)
                aCGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGaa                                                                                                                                       <  2:2931346/150‑1 (MQ=255)
                 ccgCATGGTTACCGAACCGAACAATGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGatg                                                                                                                                      <  2:3544084/148‑3 (MQ=255)
                   gCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGc                                                                                                                                    <  2:1785617/150‑1 (MQ=255)
                     aatgTTCCGGTCCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAt                                                                                                                                  >  2:513753/4‑150 (MQ=255)
                      tGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAtt                                                                                                                                 >  1:4282061/1‑150 (MQ=255)
                       ggTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAttt                                                                                                                                >  1:4134315/1‑150 (MQ=255)
                       ggTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAttt                                                                                                                                <  2:3905003/150‑1 (MQ=255)
                           ccggaccggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGa                                                                                                                            <  1:1639848/150‑1 (MQ=255)
                            cggaccggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGAc                                                                                                                           <  1:1117295/150‑1 (MQ=255)
                                ccggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAg                                                                                                                       >  2:916936/1‑150 (MQ=255)
                                 cggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAgg                                                                                                                      <  2:4463050/150‑1 (MQ=255)
                                   gaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCt                                                                                                                    >  1:105914/1‑150 (MQ=255)
                                    aaCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtg                                                                                                                   >  2:3578151/1‑150 (MQ=255)
                                     aCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgt                                                                                                                  <  1:4515325/150‑1 (MQ=255)
                                      cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgtg                                                                                                                 <  1:2712478/150‑1 (MQ=255)
                                        ttGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGtt                                                                                                               <  1:2964630/150‑1 (MQ=255)
                                        ttGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGtt                                                                                                               >  1:1293966/1‑150 (MQ=255)
                                        ttGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGtt                                                                                                               >  1:253908/1‑150 (MQ=255)
                                           tCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTtctc                                                                                                            <  2:4515486/150‑1 (MQ=255)
                                            cGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCa                                                                                                           <  1:556340/150‑1 (MQ=255)
                                              gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                         >  2:2223165/1‑150 (MQ=255)
                                              gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                         <  2:1303477/150‑1 (MQ=255)
                                                 aGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGAt                                                                                                      <  2:1618861/150‑1 (MQ=255)
                                                    ccTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGCTCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCgg                                                                                                   <  2:2801063/150‑1 (MQ=255)
                                                         gcgGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGAt                                                                                              >  2:1537137/1‑150 (MQ=255)
                                                          cgGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGAtt                                                                                             <  1:235450/150‑1 (MQ=255)
                                                               tgGATGAGCTGGCCAGAGCCCAAACTCATCCACTTGCCCAACCTGTCGCGGTTTGGCCTTATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCca                                                                                        <  1:1183497/150‑1 (MQ=255)
                                                                ggATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCcac                                                                                       <  1:652273/150‑1 (MQ=255)
                                                                 gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                          >  1:671725/1‑146 (MQ=255)
                                                                 gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                          <  2:671775/146‑1 (MQ=255)
                                                                        tGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAAt                                                                               >  1:751736/1‑150 (MQ=255)
                                                                         ggCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATGGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATc                                                                              <  1:3226085/150‑1 (MQ=255)
                                                                         ggCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATc                                                                              <  2:799390/150‑1 (MQ=255)
                                                                          gCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATcg                                                                             >  2:330894/1‑150 (MQ=255)
                                                                            aagagCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTTTTCTCATGGATCGGGAGATTCGCCACACATCAATcggc                                                                           >  1:381550/2‑150 (MQ=255)
                                                                              gaacccAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATcggcgg                                                                         <  2:1329450/147‑1 (MQ=255)
                                                                               agCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGa                                                                        <  1:388223/150‑1 (MQ=255)
                                                                                  ccAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACca                                                                     >  1:1740495/1‑150 (MQ=255)
                                                                                  ccAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGAGCGGGAGATTCGCCACACATCAATCGGCGGACca                                                                     >  2:3063839/1‑150 (MQ=255)
                                                                                    aaaCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcaca                                                                   >  1:2408533/1‑150 (MQ=255)
                                                                                     aaCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAg                                                                  <  1:1469500/150‑1 (MQ=255)
                                                                                      aCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGc                                                                 <  2:102994/150‑1 (MQ=255)
                                                                                        tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCa                                                               >  2:3768542/1‑150 (MQ=255)
                                                                                        tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCa                                                               >  1:3035602/1‑150 (MQ=255)
                                                                                          aGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                             <  1:3316600/150‑1 (MQ=255)
                                                                                            ccACTTGCCCTAGCTGTCGGGGGTTGGGCTGATCCATTGTTACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCTTTTCTGACATTGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                           <  1:3948430/150‑1 (MQ=255)
                                                                                            ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                           <  1:1652876/150‑1 (MQ=255)
                                                                                            ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                           >  1:1874160/1‑150 (MQ=255)
                                                                                              aCTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGt                                                         <  2:2599290/150‑1 (MQ=255)
                                                                                              aCTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGt                                                         >  2:3808121/1‑150 (MQ=255)
                                                                                                   cccAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTc                                                    <  1:690389/150‑1 (MQ=255)
                                                                                                    ccAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCa                                                   <  1:3037914/150‑1 (MQ=255)
                                                                                                     cAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGatgg                                                                                                                                     <  1:550916/67‑4 (MQ=18)
                                                                                                      aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                 >  2:853264/1‑150 (MQ=255)
                                                                                                      aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                 <  1:1931383/150‑1 (MQ=255)
                                                                                                      aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                 >  2:1829342/1‑150 (MQ=255)
                                                                                                         cTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATg                                              >  1:109791/1‑150 (MQ=255)
                                                                                                           gTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGtt                                            <  2:4302993/150‑1 (MQ=255)
                                                                                                              gcgGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAt                                         >  1:375601/1‑150 (MQ=255)
                                                                                                                gggTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAttt                                       <  1:2948516/150‑1 (MQ=255)
                                                                                                                   ttGGCCTGATCCATTGATACCCCTCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaa                                    <  2:4303114/150‑1 (MQ=255)
                                                                                                                   ttGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaa                                    <  2:576757/150‑1 (MQ=255)
                                                                                                                   ttGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaa                                    <  2:2998779/150‑1 (MQ=255)
                                                                                                                      gCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaa                                 <  2:2950455/150‑1 (MQ=255)
                                                                                                                       ccTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaac                                <  1:4421210/150‑1 (MQ=255)
                                                                                                                       ccTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaac                                <  1:1585258/150‑1 (MQ=255)
                                                                                                                        cTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACCAcg                               >  1:279591/1‑150 (MQ=255)
                                                                                                                        cTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACAGAGGCTGTGTGCTCAGGGATCGGGAGATTCGCCACACCTCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGCGGTTGATTTGAACAAcg                               >  2:3629460/1‑150 (MQ=255)
                                                                                                                              cATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGa                         <  1:3794900/150‑1 (MQ=255)
                                                                                                                              cATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGa                         <  1:2592979/150‑1 (MQ=255)
                                                                                                                               aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaa                        >  1:2792549/1‑150 (MQ=255)
                                                                                                                               aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaa                        >  1:266266/1‑150 (MQ=255)
                                                                                                                               aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaa                        <  1:3768542/150‑1 (MQ=255)
                                                                                                                                   aTACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaacaac                    >  1:3207569/1‑150 (MQ=255)
                                                                                                                                       ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                <  1:2313004/150‑1 (MQ=255)
                                                                                                                                       ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                >  1:762853/1‑150 (MQ=255)
                                                                                                                                       ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                <  1:339670/150‑1 (MQ=255)
                                                                                                                                        cATAGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTATGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGTACAACGATCg               >  2:3534870/1‑150 (MQ=255)
                                                                                                                                         ttCGGCGCTGTCTGCTGGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCgg              >  2:1315407/2‑150 (MQ=255)
                                                                                                                                         aTCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCgg              >  1:4233736/1‑150 (MQ=255)
                                                                                                                                           cGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGt            <  1:2178599/150‑1 (MQ=255)
                                                                                                                                           cGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGt            <  1:3788501/150‑1 (MQ=255)
                                                                                                                                                 tGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGc      >  2:4096553/1‑150 (MQ=255)
                                                                                                                                                  gTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCa     >  2:758938/1‑150 (MQ=255)
                                                                                                                                                     cagcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGa  <  1:1912433/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGA  >  NC_002947/1633817‑1634115

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