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

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

GAAGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACG  >  NC_002947/4780823‑4781118
                                                                                                                                                    |                                                                                                                                                   
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                                                                                        <  4:317593/149‑1 (MQ=255)
                                                             ccGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgc                                                                                        <  7:296662/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                                                                                       <  5:286831/106‑1 (MQ=255)
                                                                                                         ggcggcGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGgcgcg                                                                                       >  6:286831/1‑106 (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)
                                                                                                                                               gCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGt      <  6:103007/149‑1 (MQ=255)
                                                                                                                                                cccAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTc     <  3:220815/149‑1 (MQ=255)
                                                                                                                                                   cggCAATGCTCAGCGCACCCCAGCACCCCCGCCTGGCTCAAGCCCAGCACAGCGGCCTGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACg  <  2:283216/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGTGAAGATCACGTAGGCCAGCAGCTTGGCGGCATTGATGGGCAGGTGCAGCATCTGTGCCGGCGCGGCGTGTGCCTGCACCTGCTCCCAGACCAACAAGGTCGGCGGCGCTTGGAAGCCTTCGAGCATCGCTTGGGCCAAGCCCAGGCAATGCTCATCGCACACCAGCACCCGCGCCTGGCTCAAGCGCAGCACATCGGCCAGGCGCGCCACCGGCAAGGACGGATCCAGCGAGAGATAACCGGCACCCGCCTGGAACACCCCGACCACCATGCCCAGCAACGACAGGTCACG  >  NC_002947/4780823‑4781118

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: