Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,633,969 C→T 5.1% intergenic (+21/‑143) 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,9690CT 5.1% 509.6 / 10.7 156intergenic (+21/‑143)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.89e-01

Only 100 of 172 total aligned reads displayed.

CCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGT  >  NC_002947/1633820‑1634117
                                                                                                                                                     |                                                                                                                                                    
 ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                     >  1:3804068/1‑150 (MQ=255)
 ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                     >  1:507990/1‑150 (MQ=255)
 ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                     >  2:2948516/1‑150 (MQ=255)
  cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                    >  1:447305/1‑150 (MQ=255)
  cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                    >  1:814489/1‑150 (MQ=255)
  cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                    >  2:2178599/1‑150 (MQ=255)
  cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                    >  2:2270780/1‑150 (MQ=255)
   aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGa                                                                                                                                                   >  2:4237468/1‑150 (MQ=255)
         cTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAgcgc                                                                                                                                             <  1:3100205/150‑1 (MQ=255)
          tGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGAGCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                            >  1:2024903/1‑150 (MQ=255)
            cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGa                                                                                                                                          <  1:1654502/150‑1 (MQ=255)
            cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGa                                                                                                                                          <  2:1446425/150‑1 (MQ=255)
            cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGa                                                                                                                                          <  1:538435/150‑1 (MQ=255)
            cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGa                                                                                                                                          <  2:744299/150‑1 (MQ=255)
             aCGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGat                                                                                                                                         <  2:1013740/150‑2 (MQ=255)
             aCGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGaa                                                                                                                                         <  2:3580475/150‑1 (MQ=255)
              ccgCATGGTTACCGAACCGAACAATGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGatg                                                                                                                                        <  2:3544084/148‑3 (MQ=255)
               gacaatgTTCCGGTCCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATg                                                                                                                                       >  2:3277009/7‑150 (MQ=255)
                  aatgTTCCGGTCCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAt                                                                                                                                    >  2:513753/4‑150 (MQ=255)
                   tGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAtt                                                                                                                                   >  1:4282061/1‑150 (MQ=255)
                      ttCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCt                                                                                                                                <  2:2152944/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                                                                                                                        <  1:376906/150‑1 (MQ=255)
                              cggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAgg                                                                                                                        <  2:4462895/150‑1 (MQ=255)
                                gaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCt                                                                                                                      >  1:105914/1‑150 (MQ=255)
                                  aCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTGAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgt                                                                                                                    >  1:1641237/1‑150 (MQ=255)
                                   cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgtg                                                                                                                   <  1:2712478/150‑1 (MQ=255)
                                     ttGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGtt                                                                                                                 >  1:253908/1‑150 (MQ=255)
                                     ttGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGtt                                                                                                                 <  1:2964630/150‑1 (MQ=255)
                                         cGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCa                                                                                                             <  1:556340/150‑1 (MQ=255)
                                           gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                           <  2:4343441/150‑1 (MQ=255)
                                           gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                           >  1:2984670/1‑150 (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)
                                                   ttCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATGGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGga                                                                                                   <  1:4338432/150‑1 (MQ=255)
                                                      gcgGTGTGGAGGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGACTGCATTGCTGACATAGTCTGTGTTCTCAGGGATCGGgagaa                                                                                                >  2:4490123/1‑149 (MQ=255)
                                                       cgGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGAtt                                                                                               <  1:235450/150‑1 (MQ=255)
                                                            tgGATGAGCTGGCCAGAGCCCAAACTCATCCACTTGCCCAACCTGTCGCGGTTTGGCCTTATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCca                                                                                          <  1:1183497/150‑1 (MQ=255)
                                                              gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                            >  1:671725/1‑146 (MQ=255)
                                                              gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                            <  2:671725/146‑1 (MQ=255)
                                                              gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                            >  1:671775/1‑146 (MQ=255)
                                                                      ggCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATGGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATc                                                                                <  1:3226085/150‑1 (MQ=255)
                                                                      ggCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATc                                                                                <  2:799390/150‑1 (MQ=255)
                                                                            agCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGa                                                                          <  1:388223/150‑1 (MQ=255)
                                                                               ccAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACca                                                                       >  1:1740495/1‑150 (MQ=255)
                                                                                cAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcac                                                                      >  2:3673024/1‑150 (MQ=255)
                                                                                cAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcac                                                                      >  2:1323182/1‑150 (MQ=255)
                                                                                 aaaCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcaca                                                                     >  1:2408533/1‑150 (MQ=255)
                                                                                  aaCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAg                                                                    <  1:1469500/150‑1 (MQ=255)
                                                                                     tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCa                                                                 >  2:3768542/1‑150 (MQ=255)
                                                                                     tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCa                                                                 >  1:3220965/1‑150 (MQ=255)
                                                                                     tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCa                                                                 >  1:3035611/1‑150 (MQ=255)
                                                                                      cAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaa                                                                <  1:760971/150‑1 (MQ=255)
                                                                                       aGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                               >  1:694131/1‑150 (MQ=255)
                                                                                       aGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                               >  1:3293210/1‑150 (MQ=255)
                                                                                           aCTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGt                                                           >  2:3808121/1‑150 (MQ=255)
                                                                                                 ccAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCa                                                     <  1:403569/150‑1 (MQ=255)
                                                                                                 ccAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCa                                                     <  1:2073391/150‑1 (MQ=255)
                                                                                                  cAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGatgg                                                                                                                                       <  1:550916/67‑4 (MQ=18)
                                                                                                  cAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGatgg                                                                                                                                       >  2:550916/1‑64 (MQ=18)
                                                                                                   aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                   >  2:1829309/1‑150 (MQ=255)
                                                                                                   aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                   <  1:1931383/150‑1 (MQ=255)
                                                                                                   aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                   >  2:853264/1‑150 (MQ=255)
                                                                                                      cTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATg                                                >  1:109791/1‑150 (MQ=255)
                                                                                                        gTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGtt                                              <  2:4302993/150‑1 (MQ=255)
                                                                                                          cgcgGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGa                                            <  2:1940641/150‑1 (MQ=255)
                                                                                                           gcgGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAt                                           >  1:375601/1‑150 (MQ=255)
                                                                                                           gcgGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGGTGAt                                           >  1:3790039/1‑150 (MQ=255)
                                                                                                             gggTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAttt                                         <  1:2948516/150‑1 (MQ=255)
                                                                                                               gTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGa                                       >  2:1594808/1‑150 (MQ=255)
                                                                                                                ttGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaa                                      <  2:2998779/150‑1 (MQ=255)
                                                                                                                    ccTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaac                                  <  1:1585258/150‑1 (MQ=255)
                                                                                                                    ccTGATCCATTGATACCCATCGGCACTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaac                                  >  1:3882131/1‑150 (MQ=255)
                                                                                                                     cTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACCAcg                                 >  1:279591/1‑150 (MQ=255)
                                                                                                                           cATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGa                           <  1:2592979/150‑1 (MQ=255)
                                                                                                                           cATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGa                           <  1:3794900/150‑1 (MQ=255)
                                                                                                                            aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaa                          >  1:266266/1‑150 (MQ=255)
                                                                                                                            aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaa                          >  1:2792549/1‑150 (MQ=255)
                                                                                                                             ttGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaac                         <  1:1271751/150‑1 (MQ=255)
                                                                                                                                    ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                  <  1:2313004/150‑1 (MQ=255)
                                                                                                                                    ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                  <  1:339670/150‑1 (MQ=255)
                                                                                                                                    ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                  >  2:227268/1‑150 (MQ=255)
                                                                                                                                    ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                  >  1:1801674/1‑150 (MQ=255)
                                                                                                                                    ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                  >  1:84199/1‑150 (MQ=255)
                                                                                                                                     cATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCg                 >  2:122375/1‑150 (MQ=255)
                                                                                                                                     cATAGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTATGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGTACAACGATCg                 >  2:3534870/1‑150 (MQ=255)
                                                                                                                                      aTCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCgg                >  1:4233736/1‑150 (MQ=255)
                                                                                                                                         ggCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTg             >  2:2954586/1‑150 (MQ=255)
                                                                                                                                           cgcTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGc           >  2:335371/1‑150 (MQ=255)
                                                                                                                                              tGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGc        >  2:4096553/1‑150 (MQ=255)
                                                                                                                                                  cagcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGa    <  2:3944818/150‑1 (MQ=255)
                                                                                                                                                  cagcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGa    <  1:1912346/150‑1 (MQ=255)
                                                                                                                                                  cagcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGa    >  1:2222367/1‑150 (MQ=255)
                                                                                                                                                  cagcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAGTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGa    >  1:3352034/1‑150 (MQ=255)
                                                                                                                                                    gcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGt  <  2:814489/150‑1 (MQ=255)
                                                                                                                                                    gcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGt  <  1:2270780/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGT  >  NC_002947/1633820‑1634117

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: