Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,785,336 G→A 5.3% H660Y (CAC→TAC)  pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,785,3360GA 5.3% 230.1 / 2.1 75H660Y (CAC→TAC) pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (33/38);  new base A (2/2);  total (35/40)
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.13e-01

AGGTGCAGCAGACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCGGGCACGCG  >  NC_002947/4785188‑4785481
                                                                                                                                                    |                                                                                                                                                 
aGGTGCAGCAGACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTg                                                                                                                                                   >  5:151917/1‑149 (MQ=255)
 ggTGCAGCAGACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTgg                                                                                                                                                  <  6:271295/149‑1 (MQ=255)
  gTGCAGCAGACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTggg                                                                                                                                                 >  1:348599/1‑149 (MQ=255)
  gTGCAGCAGACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTggg                                                                                                                                                 >  4:350554/1‑149 (MQ=255)
     cagcagACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCaa                                                                                                                                              <  2:174743/149‑1 (MQ=255)
        cagACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGa                                                                                                                                           >  5:57765/1‑149 (MQ=255)
           aCGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCa                                                                                                                                        <  4:182709/149‑1 (MQ=255)
             ggCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACg                                                                                                                                      >  4:236913/1‑149 (MQ=255)
               ccGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGcc                                                                                                                                    >  4:144537/1‑149 (MQ=255)
                  gATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCtt                                                                                                                                 <  4:222978/149‑1 (MQ=255)
                      ccAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGc                                                                                                                             <  8:343106/149‑1 (MQ=255)
                       cAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCt                                                                                                                            >  6:137548/1‑149 (MQ=255)
                       cAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCt                                                                                                                            <  1:266883/149‑1 (MQ=255)
                       cAGCGCGCCGAACAGCACGGTGGGCCCCAGGTCCGCCGCCACGGGCGAGCCCATCGCCCGGCTGGCGGCGGCTGGCCCCTGGCGGCGCGGCCGGAGGCCCTGCCCGTAGGTGTCCAGCAGGCGGTGGGCCACGACCACGCACTTgggag                                                                                                                            >  4:232273/1‑147 (MQ=255)
                        aGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGc                                                                                                                                     <  5:367373/139‑1 (MQ=255)
                        aGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGc                                                                                                                                     >  6:367373/1‑139 (MQ=255)
                         gcgcgcCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTgg                                                                                                                          <  8:76718/149‑1 (MQ=255)
                             gcCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCgg                                                                                                                      >  1:371780/1‑149 (MQ=255)
                             gcCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCgg                                                                                                                      <  4:216198/149‑1 (MQ=255)
                               cGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTg                                                                                                                    <  4:347466/149‑1 (MQ=255)
                               cGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTg                                                                                                                    >  7:86601/1‑149 (MQ=255)
                               cGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTg                                                                                                                    <  2:371780/149‑1 (MQ=255)
                                gAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTg                                                                                                                    >  6:306249/1‑148 (MQ=255)
                                 aaCAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGGAACGACCACGCCCTTGGGCTGGCCGGTGGa                                                                                                                  <  1:30995/149‑1 (MQ=255)
                                 aaCAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGa                                                                                                                  <  6:151917/149‑1 (MQ=255)
                                 aaCAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGa                                                                                                                  <  6:57765/149‑1 (MQ=255)
                                         ggTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAg                                                                                                          <  2:348599/149‑1 (MQ=255)
                                          gtgtGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGTCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAgt                                                                                                         <  3:136449/149‑1 (MQ=255)
                                             tgCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAgtgt                                                                                                       <  8:53517/148‑1 (MQ=255)
                                             tgCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTa                                                                                                      <  7:233910/149‑1 (MQ=255)
                                               ccccAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAgtg                                                                                                        >  4:178060/1‑145 (MQ=255)
                                               ccccAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAg                                                                                                     <  3:236913/148‑1 (MQ=255)
                                                cccAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGAt                                                                                                   >  4:148412/1‑149 (MQ=255)
                                                cccAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGAt                                                                                                   <  3:253789/149‑1 (MQ=255)
                                                 ccAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATg                                                                                                  >  6:339777/1‑149 (MQ=255)
                                                      tCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATa                                                                                             <  5:258695/149‑1 (MQ=255)
                                                           cgccACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCg                                                                                        <  3:350554/149‑1 (MQ=255)
                                                               aCGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTg                                                                                    >  3:209051/1‑149 (MQ=255)
                                                               aCGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTg                                                                                    <  6:123401/149‑1 (MQ=255)
                                                                tggTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTgg                                                                                   <  7:247435/148‑1 (MQ=255)
                                                                  gTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGggc                                                                                 <  3:144537/149‑1 (MQ=255)
                                                                  gTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGggc                                                                                 <  1:29568/149‑1 (MQ=255)
                                                                  gTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTCGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTgt                                                                                   >  8:247435/1‑146 (MQ=255)
                                                                     gagaCCATCGCCATGCTGGCGTCGGCTGGCAACTGTAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTAGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGggcggc                                                                              >  5:303298/1‑149 (MQ=255)
                                                                     gagaCCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTAGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGggcggc                                                                              >  5:242091/1‑149 (MQ=255)
                                                                          cATCGCCATGTTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACTACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGa                                                                         >  3:270570/1‑149 (MQ=255)
                                                                           aTCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGAc                                                                        <  5:137548/149‑1 (MQ=255)
                                                                                cATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTgg                                                                                                                          <  4:311712/94‑1 (MQ=255)
                                                                                cATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTgg                                                                                                                          >  3:311712/1‑94 (MQ=255)
                                                                                cATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGa                                                                   <  3:148412/149‑1 (MQ=255)
                                                                                     tGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGc                                                               >  4:218973/1‑148 (MQ=255)
                                                                                       gCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACTACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCa                                                             <  4:270570/148‑1 (MQ=255)
                                                                                       gCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAg                                                            <  6:5065/149‑1 (MQ=255)
                                                                                       gCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGAAGACCTCGATGGCCAg                                                            >  7:231440/1‑149 (MQ=255)
                                                                                           cGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAgcggc                                                        <  3:218973/149‑1 (MQ=255)
                                                                                               tGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTc                                                                     <  4:32794/132‑1 (MQ=255)
                                                                                               tGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTc                                                                     >  3:32794/1‑132 (MQ=255)
                                                                                                 gCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCgtcgtc                                                  <  5:6840/149‑1 (MQ=255)
                                                                                                    aCTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCggg                                               <  4:209051/149‑1 (MQ=255)
                                                                                                     cTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCggg                                               >  8:325561/1‑148 (MQ=255)
                                                                                                        cAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACt                                           >  5:110298/1‑149 (MQ=255)
                                                                                                        cAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACt                                           >  7:79760/1‑149 (MQ=255)
                                                                                                                 cAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAgg                                  >  5:380308/1‑149 (MQ=255)
                                                                                                                   ggATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAggcg                                >  7:111435/1‑149 (MQ=255)
                                                                                                                      gccccTGCACGTAGGTGGCCAGCAGCCCGTGGGCAACCACCACGCCCTGTGGCTGGCAGGTGGAGCCGGAAGTGCAGATGACATACGCCGGGTGGGCGGCGCAGACCTCGATGGCCAGCGGCGGCGTCGGGCACTGTCGCCAGGCGgca                             <  6:324127/147‑1 (MQ=255)
                                                                                                                         ccTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACAACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGgcagca                          >  5:280741/1‑149 (MQ=255)
                                                                                                                             cACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTGGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCg                       <  7:161551/148‑1 (MQ=255)
                                                                                                                                 tAGTTGGGCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCTCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAgg                  <  2:123786/149‑1 (MQ=255)
                                                                                                                                  aGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGGCAGGCGGCAGCATCGAAGGGc                 >  4:288753/1‑149 (MQ=255)
                                                                                                                                      ggCCAGCATGCCGTAGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAAt              <  6:303298/148‑1 (MQ=255)
                                                                                                                                      ggCCAGCATGCCGTAGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAAt              <  6:242091/148‑1 (MQ=255)
                                                                                                                                         cAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCg          >  1:203713/1‑149 (MQ=255)
                                                                                                                                          aGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCgg         >  6:148803/1‑149 (MQ=255)
                                                                                                                                             aTGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCGGGc       >  5:314416/1‑148 (MQ=255)
                                                                                                                                               gCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCGGGCAcg    >  8:182345/1‑149 (MQ=255)
                                                                                                                                                ccGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCGGGCAcgc   >  1:3703/1‑149 (MQ=255)
                                                                                                                                                 cGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGc                                                                                          >  2:211464/1‑61 (MQ=255)
                                                                                                                                                 cGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGAAGATGACATAGGACGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCGGGCAcgcg  >  2:23199/1‑149 (MQ=255)
                                                                                                                                                 cGTGGGCAACGACCACGCCCTTGGGCTGGCCGCTGGAGCCGGAAGTGTAGATGACATAGGc                                                                                          <  1:211464/61‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGGTGCAGCAGACGGCCGGATGCCAGCGCGCCGAACAGCACGGTGTGCCCCAGGTCCGCCGCCACGGTCGAGACCATCGCCATGCTGGCGTCGGCTGGCAACTGCAGGCGCTGCAGGATGCCCTGCACGTAGTTGGCCAGCATGCCGTGGGCAACGACCACGCCCTTGGGCTGGCCGGTGGAGCCGGAAGTGTAGATGACATAGGCCGGTTGGGCGGCGCAGACCTCGATGGCCAGCGGCGTCGTCGGGCACTGTCGCCAGGCGGCAGCATCGAAAGGCAATACCGGGCACGCG  >  NC_002947/4785188‑4785481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: