Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,770,958 Δ1 bp 5.1% coding (10276/10413 nt) pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,770,9580G. 5.1% 589.1 / 25.1 175coding (10276/10413 nt)pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (67/99);  new base . (4/5);  total (71/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 175 total aligned reads displayed.

CCCCGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGTTCATCGACG  >  NC_002947/4770811‑4771106
                                                                                                                                                     |                                                                                                                                                    
 cccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTg                                                                                                                                                   >  1:1156432/1‑150 (MQ=255)
 cccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTg                                                                                                                                                   >  2:4243091/1‑150 (MQ=255)
 cccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTg                                                                                                                                                   >  1:4355405/1‑150 (MQ=255)
  ccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTgc                                                                                                                                                  <  1:4781856/150‑1 (MQ=255)
    gAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTgcgc                                                                                                                                                >  1:3923639/1‑150 (MQ=255)
          taaTGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:3997675/148‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:2201522/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          >  1:4171117/1‑150 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:1192248/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          >  1:83915/1‑150 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          >  1:2307902/1‑150 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:1313151/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:3022426/150‑1 (MQ=255)
           cATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCGGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctgct                                                                                                                                         <  2:3109042/150‑5 (MQ=37)
           cATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCGGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctgct                                                                                                                                         <  2:2982587/150‑5 (MQ=37)
           cATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCACATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctgct                                                                                                                                         <  1:2535804/150‑5 (MQ=37)
               aaaGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCtt                                                                                                                                     <  1:1049188/150‑1 (MQ=255)
                  gTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTcgc                                                                                                                                  <  1:1357387/150‑1 (MQ=255)
                   ttCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTcgcg                                                                                                                                 >  2:2814270/1‑150 (MQ=255)
                   ttCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTcgcg                                                                                                                                 >  2:10297/1‑150 (MQ=255)
                    tCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGATTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCTCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCgcgc                                                                                                                                <  2:4014653/150‑1 (MQ=255)
                      aaTGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTc                                                                                                                              >  2:2355000/1‑150 (MQ=255)
                      aaTGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTc                                                                                                                              >  2:5162248/1‑150 (MQ=255)
                         gAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATc                                                                                                                           <  1:5204592/150‑1 (MQ=255)
                               ggCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACg                                                                                                                     <  1:3688972/150‑1 (MQ=255)
                                gcgtgcgtGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACgg                                                                                                                    >  1:45866/1‑150 (MQ=255)
                                        gAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTAc                                                                                                            <  2:1050445/150‑1 (MQ=255)
                                           gTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACcca                                                                                                         <  1:633039/150‑1 (MQ=255)
                                            ttGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACccac                                                                                                        >  2:1491151/1‑150 (MQ=255)
                                            ttGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACccac                                                                                                        <  2:753583/150‑1 (MQ=255)
                                            ttGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACccac                                                                                                        <  2:754451/150‑1 (MQ=255)
                                               ggTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCaa                                                                                                     <  1:2734389/150‑1 (MQ=255)
                                                gTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctg                                                                                                                                           >  2:2542836/1‑109 (MQ=25)
                                                gTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctg                                                                                                                                           <  1:2542836/111‑3 (MQ=25)
                                                 ttCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACt                                                                                                   >  2:4367904/1‑150 (MQ=255)
                                                     agcagcGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGc                                                                                               <  2:496732/150‑1 (MQ=255)
                                                        agcGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGt                                                                                            >  2:1329006/1‑150 (MQ=255)
                                                         gcgcTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTa                                                                                           >  2:2893468/1‑150 (MQ=255)
                                                         gcgcTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTa                                                                                           >  2:2893613/1‑150 (MQ=255)
                                                         gcgcTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTa                                                                                           <  1:4143869/150‑1 (MQ=255)
                                                          cgcTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTac                                                                                          <  1:3155511/150‑1 (MQ=255)
                                                             ttGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTaccac                                                                                       >  2:3724356/1‑150 (MQ=255)
                                                              tGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTaccacc                                                                                      <  2:340164/150‑1 (MQ=255)
                                                                 aCGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGc                                                                                   <  1:4367904/150‑1 (MQ=255)
                                                                 aCGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGc                                                                                   >  1:1230134/1‑150 (MQ=255)
                                                                      gtcgAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACggg                                                                              <  1:2179039/150‑1 (MQ=255)
                                                                      gtcgAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACggg                                                                              >  2:1993285/1‑150 (MQ=255)
                                                                       tgacccAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑GGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGc                                                                             >  2:4591665/6‑150 (MQ=255)
                                                                       tcgAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGc                                                                             >  1:4701070/1‑150 (MQ=255)
                                                                         gAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTg                                                                           <  2:4810180/150‑1 (MQ=255)
                                                                           gCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGc                                                                         <  1:2602544/150‑1 (MQ=255)
                                                                           gCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGc                                                                         >  1:1992780/1‑150 (MQ=255)
                                                                                  ggCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGtgtg                                                                  <  1:4599654/150‑1 (MQ=255)
                                                                                  gcccTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGtgtg                                                                  <  2:597958/148‑1 (MQ=255)
                                                                                   gCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGc                                                                 >  1:783148/1‑150 (MQ=255)
                                                                                    ccTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCa                                                                <  2:783148/150‑1 (MQ=255)
                                                                                     cTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAg                                                               >  1:4629254/1‑150 (MQ=255)
                                                                                     cTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAg                                                               >  1:4629418/1‑150 (MQ=255)
                                                                                      tCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGc                                                              <  1:5023827/150‑1 (MQ=255)
                                                                                        aTGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTc                                                            >  1:3484858/1‑150 (MQ=255)
                                                                                          gCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCAc                                                          <  2:2283500/150‑1 (MQ=255)
                                                                                            ggCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCAcgc                                                        <  1:5449740/150‑1 (MQ=255)
                                                                                             gCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCAcgcg                                                       >  2:1903084/1‑150 (MQ=255)
                                                                                                 gtccaggTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAg                                                   <  2:1401753/144‑1 (MQ=255)
                                                                                                     aaGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGgcgc                                               <  1:3373706/150‑1 (MQ=255)
                                                                                                        tGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTg                                            <  2:1710699/150‑1 (MQ=255)
                                                                                                         gCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACGAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGa                                           <  2:536827/150‑1 (MQ=255)
                                                                                                         gCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGa                                           <  2:4103378/150‑1 (MQ=255)
                                                                                                         gCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGa                                           <  2:537562/150‑1 (MQ=255)
                                                                                                           ccACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCtt                                             >  1:4555571/1‑146 (MQ=255)
                                                                                                             acaGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCt                                       >  1:1867560/1‑150 (MQ=255)
                                                                                                             acaGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCt                                       <  2:4355405/150‑1 (MQ=255)
                                                                                                             acaGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCt                                       >  2:2091019/1‑150 (MQ=255)
                                                                                                             acaGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCt                                       <  1:677991/150‑1 (MQ=255)
                                                                                                               aGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgc                                     <  1:4847509/150‑1 (MQ=255)
                                                                                                               aGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgc                                     <  1:5094641/150‑1 (MQ=255)
                                                                                                                ggCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgcg                                    <  2:321684/150‑1 (MQ=255)
                                                                                                                ggCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgcg                                    <  2:2475488/150‑1 (MQ=255)
                                                                                                                 gCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgcgc                                   <  1:2948056/150‑1 (MQ=255)
                                                                                                                  ccAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgcgcg                                  <  2:3058579/150‑1 (MQ=255)
                                                                                                                     ttggcCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGCTAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAg                               >  1:1557299/4‑150 (MQ=37)
                                                                                                                      cagcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGt                              <  1:10297/150‑1 (MQ=255)
                                                                                                                      cagcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGt                              <  2:3103776/150‑1 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTgg                                                                          >  2:248702/1‑105 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTgg                                                                          <  1:248702/105‑1 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtg                             <  2:45866/150‑1 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCAGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtg                             >  1:4236232/1‑150 (MQ=255)
                                                                                                                          cgcgCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtgatg                          <  2:5339425/150‑1 (MQ=255)
                                                                                                                          cgcgCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtga                            <  1:245327/148‑1 (MQ=255)
                                                                                                                           gcgcAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtgatga                         <  2:1096839/150‑1 (MQ=255)
                                                                                                                                   tgtCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGTTGAACACCt                   <  2:3583346/150‑1 (MQ=255)
                                                                                                                                   tgtCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCt                   <  1:3498280/150‑1 (MQ=255)
                                                                                                                                     tCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACTTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTgt                 <  1:747038/150‑2 (MQ=255)
                                                                                                                                     tCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTgg                 >  1:4812953/1‑150 (MQ=255)
                                                                                                                                     tCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTgg                 >  1:2271353/1‑150 (MQ=255)
                                                                                                                                     tCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCCCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTgg                 >  2:4365822/1‑150 (MQ=255)
                                                                                                                                       cAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGt               >  2:4619988/1‑150 (MQ=255)
                                                                                                                                           gCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGt           >  1:4880708/1‑150 (MQ=255)
                                                                                                                                                gCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGTTCATc      <  2:2731708/150‑1 (MQ=255)
                                                                                                                                                  cAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGTTCATCGa    >  2:3408631/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCCGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGTTCATCGACG  >  NC_002947/4770811‑4771106

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: