Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,633,967 A→G 5.1% intergenic (+19/‑145) 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,9670AG 5.1% 519.0 / 10.4 156intergenic (+19/‑145)lepB/rncsignal peptidase I/ribonuclease 3
Reads supporting (aligned to +/- strand):  ref base A (71/77);  new base G (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 173 total aligned reads displayed.

TTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAG  >  NC_002947/1633818‑1634116
                                                                                                                                                     |                                                                                                                                                     
ttCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCca                                                                                                                                                       >  2:198758/1‑150 (MQ=255)
   ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                    >  1:3804068/1‑150 (MQ=255)
   ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                    >  1:507990/1‑150 (MQ=255)
   ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                    >  2:2948516/1‑150 (MQ=255)
    cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagcag                                                                                                                                                   >  2:2270780/1‑150 (MQ=255)
     aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGa                                                                                                                                                  <  2:2171328/150‑1 (MQ=255)
       ggAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACa                                                                                                                                                >  1:2950455/1‑150 (MQ=255)
          aCTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAgcg                                                                                                                                             <  1:2912534/150‑1 (MQ=255)
            tGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGAGCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                           >  1:2024903/1‑150 (MQ=255)
              cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGa                                                                                                                                         <  1:1654502/150‑1 (MQ=255)
              cACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGa                                                                                                                                         <  1:538435/150‑1 (MQ=255)
               aCGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCTGCTGGACAGCGCCGat                                                                                                                                        <  2:1013740/150‑2 (MQ=255)
               aCGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGaa                                                                                                                                        <  2:2931346/150‑1 (MQ=255)
               aCGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGaa                                                                                                                                        <  2:3580475/150‑1 (MQ=255)
                 gacaatgTTCCGGTCCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATg                                                                                                                                      >  2:3277009/7‑150 (MQ=255)
                    aatgTTCCGGTCCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAt                                                                                                                                   >  2:513753/4‑150 (MQ=255)
                     tGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAtt                                                                                                                                  >  1:4282061/1‑150 (MQ=255)
                      ggTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAttt                                                                                                                                 <  2:3905003/150‑1 (MQ=255)
                      ggTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAttt                                                                                                                                 >  1:4134315/1‑150 (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)
                                cggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAgg                                                                                                                       <  1:376906/150‑1 (MQ=255)
                                  gaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCt                                                                                                                     >  1:105914/1‑150 (MQ=255)
                                   aaCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtg                                                                                                                    >  2:3578151/1‑150 (MQ=255)
                                    aCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTGAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgt                                                                                                                   >  1:1641237/1‑150 (MQ=255)
                                    aCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgt                                                                                                                   <  2:2181528/150‑1 (MQ=255)
                                    aCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgt                                                                                                                   <  1:4515325/150‑1 (MQ=255)
                                       ttGTCGGCAAGGCCTTCGCGGTGTGGATGTGCTGGCCAGAGCCCAAACTCAGCCTCTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGtt                                                                                                                <  1:2474470/150‑1 (MQ=255)
                                       ttGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGtt                                                                                                                >  2:3537381/1‑150 (MQ=255)
                                          tCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTtctc                                                                                                             <  2:4515486/150‑1 (MQ=255)
                                           cGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCa                                                                                                            <  1:556340/150‑1 (MQ=255)
                                             gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                          <  2:1303477/150‑1 (MQ=255)
                                             gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                          <  2:4343441/150‑1 (MQ=255)
                                             gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                          >  2:2223165/1‑150 (MQ=255)
                                             gCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAgg                                                                                                          >  1:2984670/1‑150 (MQ=255)
                                                aGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGAt                                                                                                       <  2:1618861/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                                                                                        <  2:3770080/150‑1 (MQ=255)
                                                               ggATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCcac                                                                                        <  1:652273/150‑1 (MQ=255)
                                                                gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                           <  2:671775/146‑1 (MQ=255)
                                                                gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                           >  1:671775/1‑146 (MQ=255)
                                                                      cTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCaa                                                                                 >  1:2788152/1‑150 (MQ=255)
                                                                       tGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAAt                                                                                >  1:751736/1‑150 (MQ=255)
                                                                        ggCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATGGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATc                                                                               <  1:3226085/150‑1 (MQ=255)
                                                                        ggCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATc                                                                               <  2:799390/150‑1 (MQ=255)
                                                                           cAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATcggc                                                                            <  2:4282061/150‑1 (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)
                                                                                  cAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcac                                                                     >  2:1323182/1‑150 (MQ=255)
                                                                                  cAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcac                                                                     >  2:3673024/1‑150 (MQ=255)
                                                                                   aaaCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcaca                                                                    >  1:2408368/1‑150 (MQ=255)
                                                                                    aaCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAg                                                                   <  1:1469500/150‑1 (MQ=255)
                                                                                     aCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGc                                                                  <  2:102994/150‑1 (MQ=255)
                                                                                       tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCa                                                                >  1:3035611/1‑150 (MQ=255)
                                                                                         aGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                              >  1:694131/1‑150 (MQ=255)
                                                                                         aGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                              >  1:3293210/1‑150 (MQ=255)
                                                                                         aGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                              <  1:3316600/150‑1 (MQ=255)
                                                                                           ccACTTGCCCTAGCTGTCGGGGGTTGGGCTGATCCATTGTTACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCTTTTCTGACATTGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                            <  1:3948430/150‑1 (MQ=255)
                                                                                           ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                            >  1:1874160/1‑150 (MQ=255)
                                                                                             aCTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGt                                                          <  2:2599290/150‑1 (MQ=255)
                                                                                                 gCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCttt                                                      >  2:4233908/1‑150 (MQ=255)
                                                                                                   ccAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCa                                                    <  1:2073391/150‑1 (MQ=255)
                                                                                                   ccAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCa                                                    <  1:3037914/150‑1 (MQ=255)
                                                                                                     aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                  >  2:853264/1‑150 (MQ=255)
                                                                                                     aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                  <  1:1931383/150‑1 (MQ=255)
                                                                                                     aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                  >  2:1829342/1‑150 (MQ=255)
                                                                                                     aaCCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAgg                                                  >  2:1829309/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)
                                                                                                                  ttGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaa                                     <  2:576757/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)
                                                                                                                             cATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGa                          <  1:3794900/150‑1 (MQ=255)
                                                                                                                             cATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGa                          <  1:2592979/150‑1 (MQ=255)
                                                                                                                              aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaa                         >  1:266266/1‑150 (MQ=255)
                                                                                                                               ttGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaac                        <  1:1271751/150‑1 (MQ=255)
                                                                                                                                  aTACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaacaac                     >  1:3207569/1‑150 (MQ=255)
                                                                                                                                      ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                 >  1:84199/1‑150 (MQ=255)
                                                                                                                                      ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                 >  1:762853/1‑150 (MQ=255)
                                                                                                                                      ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                 <  1:2313004/150‑1 (MQ=255)
                                                                                                                                      ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATc                 <  1:339670/150‑1 (MQ=255)
                                                                                                                                        ttCGGCGCTGTCTGCTGGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCgg               >  1:551725/2‑150 (MQ=255)
                                                                                                                                        aTCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCgg               <  1:2665980/150‑1 (MQ=255)
                                                                                                                                        aTCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCgg               >  1:4233736/1‑150 (MQ=255)
                                                                                                                                          cGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGt             <  1:3788501/150‑1 (MQ=255)
                                                                                                                                           ggCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTg            >  2:2954586/1‑150 (MQ=255)
                                                                                                                                             cgcTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGc          >  2:335371/1‑150 (MQ=255)
                                                                                                                                                 gTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCa      >  2:758938/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)
                                                                                                                                                     |                                                                                                                                                     
TTCCCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAG  >  NC_002947/1633818‑1634116

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: