Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,770,957 Δ1 bp 5.2% coding (10277/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,9570A. 5.2% 591.9 / 24.0 174coding (10277/10413 nt)pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (66/99);  new base . (4/5);  total (70/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 = 7.58e-01

Only 100 of 177 total aligned reads displayed.

CCCCGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGTTCATCGACG  >  NC_002947/4770811‑4771106
                                                                                                                                                    |                                                                                                                                                     
ccccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCt                                                                                                                                                    <  1:644044/150‑1 (MQ=255)
 cccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTg                                                                                                                                                   >  1:1156432/1‑150 (MQ=255)
 cccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTg                                                                                                                                                   <  2:2163011/150‑1 (MQ=255)
 cccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTg                                                                                                                                                   >  1:2593993/1‑150 (MQ=255)
  ccGAAAGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTgc                                                                                                                                                  <  1:3613040/150‑1 (MQ=255)
     aaaGCTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGcgcg                                                                                                                                               >  1:509333/1‑150 (MQ=255)
         cTCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGTGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          >  2:4304939/1‑150 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCGGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctgc                                                                                                                                          >  1:2982587/1‑147 (MQ=37)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:1313151/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:1192248/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          >  1:83915/1‑150 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:2201522/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:27183/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:3022389/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          <  2:3022490/150‑1 (MQ=255)
          tCATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCaa                                                                                                                                          >  1:4171117/1‑150 (MQ=255)
           cATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCACATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctgct                                                                                                                                         <  1:2535804/150‑5 (MQ=37)
           cATGAAAGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCACATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctgct                                                                                                                                         <  1:2535848/150‑5 (MQ=37)
               aaaGTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCtt                                                                                                                                     <  1:1049188/150‑1 (MQ=255)
                  gTTCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTcgc                                                                                                                                  <  1:1357387/150‑1 (MQ=255)
                   ttCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTcgcg                                                                                                                                 >  1:3317407/1‑150 (MQ=255)
                    tCAATGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGATTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCTCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCgcgc                                                                                                                                <  2:4014653/150‑1 (MQ=255)
                        tGAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCAt                                                                                                                            >  2:4153179/1‑150 (MQ=255)
                         gAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCaa                                                                                                                                                       <  2:1453696/122‑1 (MQ=255)
                         gAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCaa                                                                                                                                                       >  1:1453696/1‑122 (MQ=255)
                         gAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATc                                                                                                                           <  1:4243091/150‑1 (MQ=255)
                         gAAGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATc                                                                                                                           <  1:5204592/150‑1 (MQ=255)
                          aaGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCGGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCg                                                                                                                          >  1:480202/1‑150 (MQ=255)
                           aGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTcg                                                                                                                                   <  1:2294935/140‑1 (MQ=255)
                           aGCAGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTcg                                                                                                                                   >  2:2294935/1‑140 (MQ=255)
                              aGGCGTGCGTGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCAc                                                                                                                      >  2:4589523/1‑150 (MQ=255)
                                gcgtgcgtGAGGTTGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACgg                                                                                                                    <  1:2352208/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:753583/150‑1 (MQ=255)
                                            ttGGGTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACccac                                                                                                        <  2:754451/150‑1 (MQ=255)
                                                gTTCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAGCTTGGCCAGGCCCTGGACA‑‑CTGCGCGGctg                                                                                                                                           >  2:2542836/1‑109 (MQ=25)
                                                 ttCGAGCAGCGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACt                                                                                                   >  2:4367904/1‑150 (MQ=255)
                                                    gAGCAGCGCTTGGACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACtgtg                                                                                                <  2:276012/150‑1 (MQ=255)
                                                     agcagcGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGc                                                                                               <  2:496732/150‑1 (MQ=255)
                                                      gcagcGCTTGTACGTCGTCGAGCAGGGTGGCCACATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCg                                                                                              >  2:3959554/1‑150 (MQ=255)
                                                        agcGCTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGt                                                                                            >  2:1329006/1‑150 (MQ=255)
                                                         gcgcTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTa                                                                                           >  2:2893613/1‑150 (MQ=255)
                                                            cTTGTACGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACtgtg                                                                                                <  1:4366349/142‑1 (MQ=255)
                                                                 aCGTCGTCGAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGc                                                                                   >  1:1230134/1‑150 (MQ=255)
                                                                      gtcgAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACggg                                                                              <  1:2179039/150‑1 (MQ=255)
                                                                       tcgAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGc                                                                             >  1:4701070/1‑150 (MQ=255)
                                                                         gAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTg                                                                           <  2:4810180/150‑1 (MQ=255)
                                                                         gAGCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTg                                                                           <  1:1993285/150‑1 (MQ=255)
                                                                           gCAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGc                                                                         <  1:2602616/150‑1 (MQ=255)
                                                                            cAGGGTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCa                                                                        <  2:2056140/150‑1 (MQ=255)
                                                                             cgggTGGCCTCATGCGGCACCTGGCAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGGACCACCTGCACGGGCTGGCAc                                                                       <  2:5049583/149‑1 (MQ=255)
                                                                               ggTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACgt                                                                     <  2:1000832/150‑1 (MQ=255)
                                                                               ggTGGCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACgt                                                                     <  2:2894750/150‑1 (MQ=255)
                                                                                   gCCTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGc                                                                 >  1:783148/1‑150 (MQ=255)
                                                                                    ccTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCa                                                                <  2:783148/150‑1 (MQ=255)
                                                                                    ccTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCa                                                                <  2:3107654/150‑1 (MQ=255)
                                                                                     cTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAg                                                               <  2:3163923/150‑1 (MQ=255)
                                                                                     cTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAg                                                               >  1:4629181/1‑150 (MQ=255)
                                                                                     cTCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAg                                                               >  1:4629418/1‑150 (MQ=255)
                                                                                      tCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGc                                                              >  2:585031/1‑150 (MQ=255)
                                                                                      tCATGCGGCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGc                                                              <  1:5023827/150‑1 (MQ=255)
                                                                                            ggCACCTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCAcgc                                                        <  1:5449740/150‑1 (MQ=255)
                                                                                                 gtccaggTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAg                                                   <  2:1401753/144‑1 (MQ=255)
                                                                                                 cTGGAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTCCCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAg                                                   >  2:2839819/1‑150 (MQ=255)
                                                                                                    gAAGTGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGgcg                                                <  1:4078885/150‑1 (MQ=255)
                                                                                                        tGCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTg                                            <  2:1710699/150‑1 (MQ=255)
                                                                                                         gCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACGAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGa                                           <  2:536827/150‑1 (MQ=255)
                                                                                                         gCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGa                                           <  2:537562/150‑1 (MQ=255)
                                                                                                         gCCCACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGa                                           <  2:4103378/150‑1 (MQ=255)
                                                                                                           ccACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCtt                                             <  2:4555571/146‑1 (MQ=255)
                                                                                                           ccACAGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCtt                                             >  1:4555571/1‑146 (MQ=255)
                                                                                                             acaGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGccc                                                                                                                 <  1:1923908/76‑1 (MQ=255)
                                                                                                             acaGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGccc                                                                                                                 >  2:1923908/1‑76 (MQ=255)
                                                                                                               aGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgc                                     <  1:4847509/150‑1 (MQ=255)
                                                                                                               aGGCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgc                                     >  1:164620/1‑150 (MQ=255)
                                                                                                                 gCCAGCAGCCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTgcgc                                   <  1:2948056/150‑1 (MQ=255)
                                                                                                                      cagcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGt                              <  1:10297/150‑1 (MQ=255)
                                                                                                                      cagcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGt                              <  2:3103776/150‑1 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTgg                                                                          >  2:248680/1‑105 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTgg                                                                          <  1:248680/105‑1 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtg                             >  1:3892546/1‑150 (MQ=255)
                                                                                                                       agcCGCGCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCAGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtg                             >  1:4236232/1‑150 (MQ=255)
                                                                                                                          cgcgCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtgatg                          <  1:4751971/150‑1 (MQ=255)
                                                                                                                          cgcgCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtgatg                          <  2:5339425/150‑1 (MQ=255)
                                                                                                                          cgcgCAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtga                            >  2:245327/1‑148 (MQ=255)
                                                                                                                           gcgcAG‑‑TGTCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGtgatga                         <  2:1096839/150‑1 (MQ=255)
                                                                                                                                   tgtCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGTTGAACACCt                   <  2:3583346/150‑1 (MQ=255)
                                                                                                                                   tgtCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCt                   <  1:3498280/150‑1 (MQ=255)
                                                                                                                                    gtCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTg                  <  2:1597877/150‑1 (MQ=255)
                                                                                                                                     tCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACTTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTgt                 <  1:747038/150‑2 (MQ=255)
                                                                                                                                     tCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTgg                 <  2:3693419/150‑1 (MQ=255)
                                                                                                                                     tCCAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCCCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTgg                 >  2:4365822/1‑150 (MQ=255)
                                                                                                                                       cAGGGCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGt               >  2:4619988/1‑150 (MQ=255)
                                                                                                                                         gggCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCa             <  1:4591665/150‑1 (MQ=255)
                                                                                                                                           gCCTGGCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGt           >  1:4880708/1‑150 (MQ=255)
                                                                                                                                                gCCAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGTTCATc      <  2:3317407/150‑1 (MQ=255)
                                                                                                                                                  cAAGCTGCGCGGCCAAGGCTTCGCGCTCATCGTCACGGCCCGCTACCCACCAACTGTGCGGTACCACCTGCACGGGCTGGCACGTGTGCAGCTCACGCGACAGGCGCTTGAGGCTGCGCGCAGTGATGAACACCTGGGTCAGTTCATCGa    >  1:3374091/1‑150 (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: