Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,510,710 G→C 11.2% L115L (CTG→CTC urtA → urea ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,510,7100GC11.2% 864.1 / 79.6 277L115L (CTG→CTCurtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (118/128);  new base C (12/19);  total (130/147)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 279 total aligned reads displayed.

GGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAG  >  NC_002947/5510561‑5510859
                                                                                                                                                     |                                                                                                                                                     
ggCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTg                                                                                                                                                       <  1:8107802/150‑1 (MQ=255)
     gAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGgt                                                                                                                                                  <  2:488822/150‑1 (MQ=255)
     gAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGgt                                                                                                                                                  >  1:260875/1‑150 (MQ=255)
         ggCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTc                                                                                                                                              >  2:4765491/1‑150 (MQ=255)
           cAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGcc                                                                                                                                                     <  2:8840135/141‑1 (MQ=255)
            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                           <  1:1400698/150‑1 (MQ=255)
            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                           >  2:5849319/1‑150 (MQ=255)
            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                           >  2:4240227/1‑150 (MQ=255)
            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                           <  2:3372884/150‑1 (MQ=255)
            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                           <  2:2369260/150‑1 (MQ=255)
            aaGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCgaa                                                                                                                                           <  1:7122292/150‑1 (MQ=255)
                     gAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAAc                                                                                                                                  >  1:4817446/1‑150 (MQ=255)
                         cGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGcc                                                                                                                              <  2:6427939/150‑1 (MQ=255)
                          ggtggtggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTTAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCct                                                                                                                             <  1:3980972/150‑1 (MQ=255)
                           gtggtggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCctg                                                                                                                            >  1:3917769/1‑150 (MQ=255)
                             ggtggtCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACACAGTGGCCGTGGTGTTCGGCTGCTGGACATCGGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCctgga                                                                                                                          <  2:7024649/150‑3 (MQ=255)
                                     aCCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTAcc                                                                                                                  >  1:4217760/1‑150 (MQ=255)
                                         ggCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGt                                                                                                              <  1:2992559/150‑1 (MQ=255)
                                         ggCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGt                                                                                                              <  1:5653833/150‑1 (MQ=255)
                                          gCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTg                                                                                                             <  1:325352/150‑1 (MQ=255)
                                            gTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCa                                                                                                           <  1:5184658/150‑1 (MQ=255)
                                              ccgattGGCCGCTCTTCGCTGAAAAGGCGCGCCAGCTGATTTCTCAGGATAAGGTCGCAGCCGTGTTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGt                                                                                                         <  1:8042894/145‑1 (MQ=255)
                                                aaCTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTAt                                                                                                       >  1:7933091/1‑150 (MQ=255)
                                                aaCTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTAt                                                                                                       <  2:3745156/150‑1 (MQ=255)
                                                  cTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGa                                                                                                     <  2:2955395/150‑1 (MQ=255)
                                                   tGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAg                                                                                                    <  1:2827508/150‑1 (MQ=255)
                                                   tGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAg                                                                                                    <  1:5251285/150‑1 (MQ=255)
                                                    ggCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAgg                                                                                                   <  1:8867098/150‑1 (MQ=255)
                                                         cTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGaa                                                                                              <  1:6608531/150‑1 (MQ=255)
                                                         cTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGaa                                                                                              <  1:6733824/150‑1 (MQ=255)
                                                         cTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGaa                                                                                              <  2:5228460/150‑1 (MQ=255)
                                                          tGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAag                                                                                             >  2:8656440/1‑150 (MQ=255)
                                                           gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                                                            <  1:5883345/150‑1 (MQ=255)
                                                           gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                                                            <  2:2452890/150‑1 (MQ=255)
                                                           gTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAaga                                                                                            <  2:3091134/150‑1 (MQ=255)
                                                              cGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGAt                                                                                         >  2:2781394/1‑150 (MQ=255)
                                                               gCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATg                                                                                        >  2:1192772/1‑150 (MQ=255)
                                                                cTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGt                                                                                       >  2:3292220/1‑150 (MQ=255)
                                                                cTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGt                                                                                       <  1:4066880/150‑1 (MQ=255)
                                                                cTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGt                                                                                       >  1:9031709/1‑150 (MQ=255)
                                                                 tGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTc                                                                                      >  1:4983241/1‑150 (MQ=255)
                                                                 tGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTc                                                                                      >  1:3736860/1‑150 (MQ=255)
                                                                 tGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTc                                                                                      >  2:8456185/1‑150 (MQ=255)
                                                                   aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGc                                                                                    <  2:1973812/150‑1 (MQ=255)
                                                                   aaaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGc                                                                                    <  1:3803313/150‑1 (MQ=255)
                                                                     aaGAGCCGTCAACTGCTGACCCAGGACCAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCg                                                                                  >  2:4498434/1‑150 (MQ=255)
                                                                       gagCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGac                                                                                >  2:7904814/1‑149 (MQ=255)
                                                                        agCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAAc                                                                               >  1:3825401/1‑150 (MQ=255)
                                                                          ccGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACgt                                                                             >  2:1622652/1‑150 (MQ=255)
                                                                            gTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACgtgt                                                                           >  2:4056532/1‑150 (MQ=255)
                                                                               aaCTGCTGTCCCTGGTCAAAGTGGCCGTGGTGTTCGGCTGGTGGACATCCGTGTCGCGCAAGTCTGTTCTGGCGGTGTTCGAAGAACTCAACGGCCTGCTTTTCTACCCGGTGCAGTTTGAGGGCGAAGAGCTGTCGCCGAACGTGTTCt                                                                        <  2:1533082/150‑1 (MQ=255)
                                                                                  tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTaca                                                                     >  1:2431019/1‑150 (MQ=255)
                                                                                       aCCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgc                                                                <  2:6215037/150‑1 (MQ=255)
                                                                                         ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                                                              >  1:6788706/1‑150 (MQ=255)
                                                                                         ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                                                              >  2:6396603/1‑150 (MQ=255)
                                                                                          cAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGcgcc                                                             >  2:1040046/1‑150 (MQ=255)
                                                                                          cAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGcgcc                                                             >  1:8890816/1‑150 (MQ=255)
                                                                                          cAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGcgcc                                                             >  1:5777569/1‑150 (MQ=255)
                                                                                           aGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGcgccg                                                            >  2:7480813/1‑150 (MQ=255)
                                                                                           aGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGcgccg                                                            >  1:7195317/1‑150 (MQ=255)
                                                                                              aCAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGccc                                                         <  1:6076969/150‑1 (MQ=255)
                                                                                                aaaGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCg                                                       >  1:2797766/1‑150 (MQ=255)
                                                                                                   gTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGATGAACTCAACGGCGTGCTGTTCTACCCTGTGCAGTATGAGGGCGATGAGATGTCGGCGAACGTGTTCTACACAGGCGCCGCCCCGAAc                                                    <  2:4302549/150‑1 (MQ=255)
                                                                                                   gTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAc                                                    <  1:1736233/150‑1 (MQ=255)
                                                                                                    tGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCAAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAcc                                                   >  2:7956552/1‑150 (MQ=255)
                                                                                                     ggCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCa                                                  <  2:6874821/150‑1 (MQ=255)
                                                                                                     ggCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCa                                                  <  1:848016/150‑1 (MQ=255)
                                                                                                       ccGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAATGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGc                                                <  2:5556755/150‑1 (MQ=255)
                                                                                                         gtggtgTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTATACAGGCGCCGCCCCGAACCAGCaa                                              <  1:4508808/150‑1 (MQ=255)
                                                                                                         gtggtgTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCaa                                              >  1:1137402/1‑150 (MQ=255)
                                                                                                         gccgtgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCag                                              <  1:4787977/147‑2 (MQ=255)
                                                                                                          tggtgTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAg                                             >  2:7640215/1‑150 (MQ=255)
                                                                                                             tgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGcaa                                          >  1:3318181/1‑148 (MQ=255)
                                                                                                             tgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGcaa                                          >  2:1349927/1‑148 (MQ=255)
                                                                                                             tgtTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCa                                          >  1:5622904/1‑150 (MQ=255)
                                                                                                                       ctgGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGt                                <  1:5553595/150‑1 (MQ=255)
                                                                                                                       ctgGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCGACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGt                                <  2:8721148/150‑1 (MQ=255)
                                                                                                                        tgGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGtt                               <  1:2090360/150‑2 (MQ=255)
                                                                                                                        tgGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGtt                               >  2:8403409/1‑149 (MQ=255)
                                                                                                                        tgGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTg                               >  2:5995975/1‑150 (MQ=255)
                                                                                                                         ggACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTgg                              >  1:8629759/1‑150 (MQ=255)
                                                                                                                          gACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                             >  2:8924606/1‑150 (MQ=255)
                                                                                                                             aTCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATa                          >  1:2161645/1‑150 (MQ=255)
                                                                                                                             aTCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATa                          >  2:1466066/1‑150 (MQ=255)
                                                                                                                               ccGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGTGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATAcc                        <  1:1132499/150‑1 (MQ=255)
                                                                                                                               ccGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGACTAcc                        >  1:969568/1‑150 (MQ=255)
                                                                                                                                cGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCt                       <  2:3141162/150‑1 (MQ=255)
                                                                                                                                  tgtCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtga                     <  2:976849/150‑1 (MQ=255)
                                                                                                                                   gtCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgat                    <  2:3842044/150‑1 (MQ=255)
                                                                                                                                   gtCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgat                    <  2:8092018/150‑1 (MQ=255)
                                                                                                                                   gtCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGTTGACTATCtgat                    <  1:4859433/150‑1 (MQ=255)
                                                                                                                                    tCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGTTGACTATCtgatg                   <  2:6669703/150‑1 (MQ=255)
                                                                                                                                       cgcAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGc                <  1:8767688/150‑1 (MQ=255)
                                                                                                                                           aGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAgg            <  2:5077540/150‑1 (MQ=255)
                                                                                                                                             tCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGAc          >  2:7044679/1‑150 (MQ=255)
                                                                                                                                               tGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGAcgg        >  1:694867/1‑150 (MQ=255)
                                                                                                                                                 ttCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                             >  2:7775467/1‑127 (MQ=255)
                                                                                                                                                 ttCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                             <  1:7775467/127‑1 (MQ=255)
                                                                                                                                                 ttCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGAcggcg      >  2:4677618/1‑150 (MQ=255)
                                                                                                                                                     gCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAg  <  1:6635772/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGTGAACGGCAAGATGCTGGAACCGGTGGTGGTCGACCCGGCGTCCAACTGGCCGCTGTTCGCTGAAAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAG  >  NC_002947/5510561‑5510859

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: