Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,780,965 A→G 6.9% L2117L (TTG→CTG) ‡ pvdJ ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,780,9650AG 6.9% 165.2 / 2.5 58L2117L (TTG→CTG) ‡pvdJnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base G (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

GAACCAGAAGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAG  >  NC_002947/4780817‑4781112
                                                                                                                                                    |                                                                                                                                                   
gAACCAGAAGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCaa                                                                                                                                                     <  7:233130/149‑1 (MQ=255)
     aGAAGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCa                                                                                                                                                <  8:234564/149‑1 (MQ=255)
      gAAGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAg                                                                                                                                               <  1:31986/149‑1 (MQ=255)
       aaGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAgg                                                                                                                                              <  1:253558/149‑1 (MQ=255)
        aGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAgg                                                                                                                                              <  3:19588/148‑1 (MQ=255)
        aGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGGCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGc                                                                                                                                             <  7:198444/149‑1 (MQ=255)
         gTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGc                                                                                                                                             >  7:225786/1‑148 (MQ=255)
            aaGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAAt                                                                                                                                          <  2:347112/148‑1 (MQ=255)
            aaGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATg                                                                                                                                         >  8:296662/1‑149 (MQ=255)
             aGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGc                                                                                                                                        <  2:33320/149‑1 (MQ=255)
               aTCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATTGCCTGGGCTTGGCCCAAGCGATGCTc                                                                                                                                      <  4:96010/149‑1 (MQ=255)
                tCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGTCGCGGCGTGTGCCTGCAGCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATTGCCTGGGCTTGGCCCAAGCGATGCTcg                                                                                                                                     <  8:164508/149‑2 (MQ=255)
                     tAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGca                                                                                                                                >  5:383386/1‑149 (MQ=255)
                          cagcagCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCacac                                                                                                                             >  6:118689/1‑147 (MQ=255)
                           agcagcTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAg                                                                                                                          <  5:151728/149‑1 (MQ=255)
                           agcagcTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAg                                                                                                                          >  5:216651/1‑149 (MQ=255)
                                 ttGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCcg                                                                                                                    >  7:41332/1‑149 (MQ=255)
                                     cggcATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGcac                                                                                                                               <  2:24037/134‑1 (MQ=255)
                                     cggcATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGcac                                                                                                                               >  1:24037/1‑134 (MQ=255)
                                         aTTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGAACCTGCTCCCAGACCAACAAGGTCGGCGGAGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACAAGCACCCGCGCCTGGc                                                                                                            >  7:47172/1‑149 (MQ=255)
                                          ttGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCt                                                                                                           <  7:295595/149‑1 (MQ=255)
                                              tGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCaa                                                                                                        <  2:23366/148‑1 (MQ=255)
                                                 gCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAgcgc                                                                                                    <  8:225786/149‑1 (MQ=255)
                                                 gCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAgcg                                                                                                     >  1:272070/1‑148 (MQ=255)
                                                   aGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCgcag                                                                                                  >  3:317593/1‑149 (MQ=255)
                                                                 tgCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAg                                                                                      <  6:362465/147‑1 (MQ=255)
                                                                   ccGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgc                                                                                  <  7:296662/149‑1 (MQ=255)
                                                                   ccGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgc                                                                                  <  4:317593/149‑1 (MQ=255)
                                                                             gtgtGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGc                                                                         >  2:322070/1‑148 (MQ=255)
                                                                                tgCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAgg                                                                     >  4:104724/1‑149 (MQ=255)
                                                                                      cACCTGCTCCCAGACCAACACGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGa                                                                <  6:357020/148‑1 (MQ=255)
                                                                                      cACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGAt                                                               <  5:118689/149‑1 (MQ=255)
                                                                                      cACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGAt                                                               >  4:65924/1‑149 (MQ=255)
                                                                                         cTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGCAGCCTTCGAGCATCGCTTGGGCCAAGACCAGGCAATGCTCATCGCCCACCAGCACCCGCGCCTGGCGCAAGCGCAGCACATCGGCCAGGCGCGCCCCCGGCAAGGACGGATc                                                              >  3:359002/1‑147 (MQ=255)
                                                                                         cTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAg                                                                                                       <  5:304011/106‑1 (MQ=255)
                                                                                         cTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAg                                                                                                       >  6:304011/1‑106 (MQ=255)
                                                                                            cTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTGCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGAACACCAGCACCCGCGCGTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAg                                                           <  7:246433/147‑1 (MQ=255)
                                                                                               ccAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCgaga                                                      <  7:213489/149‑1 (MQ=255)
                                                                                                cAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCga                                                        >  8:213489/1‑146 (MQ=255)
                                                                                                 aGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCga                                                        >  2:221483/1‑145 (MQ=255)
                                                                                                 aGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCga                                                        >  2:279014/1‑145 (MQ=255)
                                                                                                    ccAACACGGGGGGGGGCGCTGGTAAGCCTGCGCGCATCGCTTGGGCCAACCCCAGGAAATGCTCATCCCACCCCCGCCCCCGCGCCGGGCTCAAGCGCAGCACATCGGCCTGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATaa                                                 <  8:41332/149‑1 (MQ=255)
                                                                                                               ggcggcGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgcg                                                                                 >  6:286831/1‑106 (MQ=255)
                                                                                                               ggcggcGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgcg                                                                                 <  5:286831/106‑1 (MQ=255)
                                                                                                                 cggcgCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCt                                    <  6:383386/149‑1 (MQ=255)
                                                                                                                     gcTTGGAACCGTTCGTGCATCGTTTGGGCAACGCCAAGTAATTGCTCATCGCACACAACCACCCGCGCCTGGCTCACGCCCAGCACAGCGGCCCGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGaa                                <  8:47172/149‑1 (MQ=255)
                                                                                                                       ttGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAaca                              >  7:86849/1‑149 (MQ=255)
                                                                                                                        tGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAacac                             >  4:213167/1‑149 (MQ=255)
                                                                                                                           aaGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACAcccc                          <  2:272070/149‑1 (MQ=255)
                                                                                                                                ttCGAGCATCGCTTGGGCCAAGCCCAGGCAATGATCATCGCACACCAGCACACCCGCCTGGCTCAAGCGCAGCACAACGGCCAGCAGCCCCACCCGCAAGGACGGATCCAGCGACAGATAACCCGCACCAGCCAGGAACACCCCGacca                     >  5:55221/1‑149 (MQ=255)
                                                                                                                                  tgAGCATTGCCTGGGCTTGGCCCAAGCGATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGaccac                    >  2:291514/2‑148 (MQ=255)
                                                                                                                                  cGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGacca                     <  1:279014/147‑1 (MQ=255)
                                                                                                                                   gAGCATTGCCTGGGCTTGGCCCAAGCGATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCCGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGAccacca                  >  3:97636/1‑149 (MQ=255)
                                                                                                                                      cATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGc               <  3:213167/149‑1 (MQ=255)
                                                                                                                                            ttGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGc          >  5:103007/1‑148 (MQ=255)
                                                                                                                                                gCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGa     >  4:220815/1‑149 (MQ=255)
                                                                                                                                                 ccAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGa     >  4:184474/1‑148 (MQ=255)
                                                                                                                                                  cAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACa   <  1:322070/149‑1 (MQ=255)
                                                                                                                                                  cAAGCCCAGGCAATGCTCATCGCACACCACCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACa   <  8:279384/149‑1 (MQ=255)
                                                                                                                                                  aaaGCCCAGGAAATGCTATGCCCACCCCCCCCCCGGCCCCTGGCTCACGCGCAGAACATCGCCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACa   <  6:216651/148‑1 (MQ=255)
                                                                                                                                                   aaGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAg  >  1:283216/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAACCAGAAGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAG  >  NC_002947/4780817‑4781112

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