Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,779,408 A→T 10.4% intergenic (‑142/+19) pvdD ← / ← pvdJ non‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,779,4080AT10.4% 170.6 / 9.0 65intergenic (‑142/+19)pvdD/pvdJnon‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (4/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAG  >  NC_002947/4779260‑4779556
                                                                                                                                                    |                                                                                                                                                    
cTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCa                                                                                                                                                      <  8:198899/149‑1 (MQ=255)
      cTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCgg                                                                                                                                                >  4:267008/1‑149 (MQ=255)
         gAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGc                                                                                                                                              <  4:262620/148‑1 (MQ=255)
          aaCATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTATGACCGGGccc                                                                                                                                            <  7:127842/149‑1 (MQ=255)
            cATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCg                                                                                                                                           >  5:64324/1‑148 (MQ=255)
             aTCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCa                                                                                                                                         <  2:320985/149‑1 (MQ=255)
                 ccTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAgc                                                                                                                                     >  6:186772/1‑149 (MQ=255)
                  cTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCagag                                                                                                                                      <  3:320066/147‑1 (MQ=255)
                  cTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCACaga                                                                                                                                       <  7:287248/146‑1 (MQ=255)
                       aGGGCACCGGTCTTTGCCTGCGGCACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCGGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACt                                                                                                                               <  3:253630/149‑1 (MQ=255)
                       aGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACt                                                                                                                               >  6:252893/1‑149 (MQ=255)
                        gggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTg                                                                                                                              <  5:286348/149‑1 (MQ=255)
                        gggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTg                                                                                                                              <  5:186772/149‑1 (MQ=255)
                         ggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTgg                                                                                                                             <  7:308035/149‑1 (MQ=255)
                           cACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATAcc                                                                                                                                                  >  3:70361/1‑126 (MQ=255)
                           cACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGCACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATAcc                                                                                                                                                  <  4:70361/126‑1 (MQ=255)
                            aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTc                                                                                                                          <  7:191782/149‑1 (MQ=255)
                            aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTATGACCGGGcc                                                                                                                                             <  7:177411/130‑1 (MQ=255)
                            aCCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTATGACCGGGcc                                                                                                                                             >  8:177411/1‑130 (MQ=255)
                               ggTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAAc                                                                                                                       <  3:5066/149‑1 (MQ=255)
                                  cTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACtcg                                                                                                                    >  3:222113/1‑149 (MQ=255)
                                      gCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAgcgc                                                                                                                                   >  6:193618/1‑130 (MQ=255)
                                      gCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAgcgc                                                                                                                                   <  5:193618/130‑1 (MQ=255)
                                            ggAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATctcct                                                                                                          <  8:128744/149‑1 (MQ=255)
                                                    ccAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGt                                                                                                   <  4:222113/148‑1 (MQ=255)
                                                                  gCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggc                                                                                    >  8:234126/1‑149 (MQ=255)
                                                                         atgCTGCCGAACATATGAAAGAGCCGGACGGCTGGCGGCTCTGGGTGTTCTGAAGTTGACGTGCCGGGCCCTTTCATACCGTGCCCGCAGAGCGCTACTGGGTCAACTCGTCTATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCaa                                                                              <  5:105264/148‑1 (MQ=255)
                                                                         atgCTGCCGAACAAATGAACGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCaa                                                                              <  3:267008/148‑1 (MQ=255)
                                                                         atgCTGCCGAACAAATGAAAGAGCCGGACTGCTGCCGGCTCTGGGTGTAATGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTAAACTAGTCGATCTCCTCTGCAGTAAGACGTTTCAAGgccgcc                                                                                >  6:105264/1‑146 (MQ=255)
                                                                         atgCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCaa                                                                              <  5:252893/148‑1 (MQ=255)
                                                                           gctgcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGa                                                                           >  4:79344/1‑149 (MQ=255)
                                                                             tgcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAtt                                                                         >  3:246350/1‑149 (MQ=255)
                                                                               ccGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAtt                                                                         <  3:79344/147‑1 (MQ=255)
                                                                               ccGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGGCATACCCGGCCCCCACCGCCGTACTGGGTCAACACGACGATCTCCTGAGCACTAAGAAGTTTCAAGGCGGCCAACGATTTa                                                                       >  6:231632/1‑149 (MQ=255)
                                                                                  aaCAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGc                                                                     >  8:13452/1‑148 (MQ=255)
                                                                                   aCAAATTAAAGAGCCGGACGGCTGCGGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCGCGGTCATACCGGTCCCGCAGAGCGATACTGGGTCAACTCGTGGATCTCCTCTGCATTAAGACGTTTCAAGGCGGCCATGGATTTAGcc                                                                    <  2:330470/148‑1 (MQ=255)
                                                                                     aaaTGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAAc                                                                 >  3:63530/1‑149 (MQ=255)
                                                                                         gAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGAAGTTTCAAGGCGGCCAAGGATTTAGCCAACTCg                                                              >  7:230644/1‑148 (MQ=255)
                                                                                         gAAAGAGCCGGACGGCTGACGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGTGTAAACTAGTCGATCTCCTCTGCAGTAAGACTTTTCAATGCGGCCAAGCATTTAGACAACTCg                                                              >  2:273877/1‑148 (MQ=255)
                                                                                           aaGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTcc                                                           <  4:63530/149‑1 (MQ=255)
                                                                                            agagCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCt                                                          <  8:193089/149‑1 (MQ=255)
                                                                                                   gACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTGATACCGGGCCCGCAGAGCGCTACTGGGTCAACACGTCGATCTCCTCTGGAGTTAGACGTTTCAAGGCGGCCAAGGATTTAGGCAACTCGTCCTGAATctc                                                   >  6:213758/1‑147 (MQ=255)
                                                                                                      ggCTGCCGGCTCTGGGTGTAATGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTGGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAGTTAGCCAACTCGTcc                                                           >  3:48373/1‑138 (MQ=255)
                                                                                                        cTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc                                              >  3:272214/1‑149 (MQ=255)
                                                                                                         tGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgcc                                             >  2:246705/1‑149 (MQ=255)
                                                                                                           ccGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccct                                           <  8:70406/149‑1 (MQ=255)
                                                                                                               ctctGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgccc                                       >  7:185379/1‑149 (MQ=255)
                                                                                                                   tggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccct                                           <  7:89732/140‑1 (MQ=255)
                                                                                                                   tggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccct                                           >  8:89732/2‑141 (MQ=255)
                                                                                                                    ggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCtt                                     <  7:234126/146‑1 (MQ=255)
                                                                                                                    ggTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGc                                   >  5:30135/1‑148 (MQ=255)
                                                                                                                     gtgtACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCaa                                 >  4:83407/1‑149 (MQ=255)
                                                                                                                          cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGCTCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTg                            >  5:65620/1‑149 (MQ=255)
                                                                                                                          cTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACt                             <  6:64324/148‑1 (MQ=255)
                                                                                                                           tGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGc                           <  1:273877/149‑1 (MQ=255)
                                                                                                                             aaGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGa                         <  6:30135/149‑1 (MQ=255)
                                                                                                                               gTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACg                       <  5:231632/149‑1 (MQ=255)
                                                                                                                                ttGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACgg                      >  5:83155/1‑149 (MQ=255)
                                                                                                                                    ctctgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTg                            >  5:95817/6‑139 (MQ=255)
                                                                                                                                    ctctgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTg                            <  6:95817/134‑1 (MQ=255)
                                                                                                                                    cGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTc                  <  4:272214/149‑1 (MQ=255)
                                                                                                                                       tgcGGGCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCt               >  8:143918/3‑149 (MQ=255)
                                                                                                                                       gCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGc                <  8:230644/148‑1 (MQ=255)
                                                                                                                                          gggCCCGGTATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGaa            >  5:294135/1‑149 (MQ=255)
                                                                                                                                              ccGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCt        >  4:211024/1‑149 (MQ=255)
                                                                                                                                                ggTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGGTTCACGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGCCGCGCCCTGCCCTTGCAACCCTGCGACGGCCTCGCTGAACGCTcc      >  3:240384/1‑149 (MQ=255)
                                                                                                                                                    aTACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAg  <  4:246350/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTGCACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCATACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGCTCCCAAG  >  NC_002947/4779260‑4779556

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