Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,826,283 A→C 5.9% I1636S (ATC→AGC) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,826,2830AC 5.9% 215.8 / 4.4 68I1636S (ATC→AGC) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base C (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATG  >  NC_002947/4826136‑4826429
                                                                                                                                                    |                                                                                                                                                      
gCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGa                                                                                                                                                        <  2:320844/148‑1 (MQ=255)
     gAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCcg                                                                                                                                                  >  5:282960/1‑149 (MQ=255)
      aaCAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGAt                                                                                                                                                       <  8:269088/143‑1 (MQ=255)
      aaCAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGAt                                                                                                                                                       >  7:269088/1‑143 (MQ=255)
           ggCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCt                                                                                                                                            >  6:235375/1‑149 (MQ=255)
            gCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCtg                                                                                                                                           <  1:288215/149‑1 (MQ=255)
                  aGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGCCCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGttgttg                                                                                                                                     <  8:219496/149‑1 (MQ=255)
                      gCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGc                                                                                                                                 >  7:105504/1‑149 (MQ=255)
                         gCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCg                                                                                                                              <  5:235375/149‑1 (MQ=255)
                          cTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCgg                                                                                                                             >  4:174042/1‑149 (MQ=255)
                          cTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCgg                                                                                                                             >  3:151666/1‑149 (MQ=255)
                             gACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGTCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCTGTTc                                                                                                                          >  6:20492/1‑149 (MQ=255)
                             gACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGtt                                                                                                                           <  3:151827/148‑1 (MQ=255)
                              aCATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCg                                                                                                                         >  1:78788/1‑149 (MQ=255)
                                aTCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg                                                                                                                       <  4:346200/149‑1 (MQ=255)
                                 tCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg‑‑‑g                                                                                                                   <  8:316659/149‑1 (MQ=255)
                                 tCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg‑‑‑g                                                                                                                   <  3:174042/149‑1 (MQ=255)
                                 tCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATACAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg‑‑‑g                                                                                                                   >  2:199038/1‑149 (MQ=255)
                                   gCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCt                                                                                                                 <  8:195632/149‑1 (MQ=255)
                                     aGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCg                                                                                                               <  2:78788/149‑1 (MQ=255)
                                           ggCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGc                                                                                                         >  1:240094/1‑149 (MQ=255)
                                               ggTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCAc                                                                                                     >  6:185332/1‑149 (MQ=255)
                                                 tGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACgc                                                                                                   >  7:116212/1‑149 (MQ=255)
                                                     cgAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACgctg‑ct                                                                                              <  5:204516/149‑1 (MQ=255)
                                                     cgAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACgctg‑ct                                                                                              <  5:185332/149‑1 (MQ=255)
                                                      gAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGctg‑ctg                                                                                             >  7:133259/1‑149 (MQ=255)
                                                           gATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GATCGACATGCTCACGCTG‑CTGcccaa                                                                                        >  7:196171/1‑146 (MQ=255)
                                                               tGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTgg                                                                                    <  4:133885/149‑1 (MQ=255)
                                                                 ggTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑Gc                                                                                                                  <  3:324515/118‑1 (MQ=255)
                                                                 ggTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑Gc                                                                                                                  >  4:324515/1‑118 (MQ=255)
                                                                 ggTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCacac                                                                                       <  8:116212/144‑1 (MQ=255)
                                                                       cAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTGCCTGCCCCAGTGTggcagca                                                                             <  8:105504/149‑3 (MQ=255)
                                                                       cAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGggcaggca                                                                            <  6:282960/149‑1 (MQ=255)
                                                                          cAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCcgc                                                                                                                                                 >  6:294958/1‑81 (MQ=255)
                                                                          cAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGAAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTGCCTGCCCCAGTGTggcagcagcg                                                                          <  8:245886/149‑6 (MQ=255)
                                                                           aaCGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGc                                                                        <  8:49838/149‑1 (MQ=255)
                                                                                gTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGc                                                                   >  6:258599/1‑149 (MQ=255)
                                                                                      accGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGa                                                             <  6:206403/149‑1 (MQ=255)
                                                                                          aGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCtt                                                         >  3:270319/1‑149 (MQ=255)
                                                                                            ggCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGc                                                       >  2:46703/1‑149 (MQ=255)
                                                                                               gaagaaGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGc                                                    >  7:33904/1‑149 (MQ=255)
                                                                                                aagaagTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCt                                                   >  4:352490/1‑149 (MQ=255)
                                                                                                aagaagTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGGCGACCTTGCCGCt                                                   <  1:199038/149‑1 (MQ=255)
                                                                                                  gaagTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTgg                                                 >  1:31324/1‑149 (MQ=255)
                                                                                                         tcCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCaa                                          >  6:353677/1‑149 (MQ=255)
                                                                                                           ctgctgGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATg                                        <  4:270319/149‑1 (MQ=255)
                                                                                                            tGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATgg                                       >  4:305278/1‑149 (MQ=255)
                                                                                                             gCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATgg                                       >  8:202900/1‑148 (MQ=255)
                                                                                                                    aaTGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCtgtg                               >  6:86545/1‑149 (MQ=255)
                                                                                                                    aaTGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGAGGATGGTCGACATCCTCACGCTC‑CTGCCACACTCGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGCCATCt                                  >  1:289846/1‑149 (MQ=255)
                                                                                                                        cgcgGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCa                            >  4:310706/1‑148 (MQ=255)
                                                                                                                          cgGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAg                         <  2:31324/149‑1 (MQ=255)
                                                                                                                              aGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGAt                     <  5:176544/149‑1 (MQ=255)
                                                                                                                              aGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGAt                     <  1:46703/149‑1 (MQ=255)
                                                                                                                                 ttCTACACGTAGC‑GCCAGATTGCCGCGACCTGTTGTTGCAGCTCGTTTCGTT‑‑‑GCTCGACATGCTCACGCTG‑CTGACACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAATGGATGa                   <  8:196171/148‑1 (MQ=255)
                                                                                                                                 ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATc                                   >  2:300311/1‑132 (MQ=255)
                                                                                                                                 ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATc                                   <  1:300311/132‑1 (MQ=255)
                                                                                                                                 ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAg                  >  5:354830/1‑149 (MQ=255)
                                                                                                                                 ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAg                  <  3:310706/149‑1 (MQ=255)
                                                                                                                                  acAACAGGTCGCGGCCATCTGG‑CGCGACGTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTggg                                                >  3:152443/2‑118 (MQ=37)
                                                                                                                                  acAACAGGTCGCGGCCATCTGG‑CGCGACGTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTggg                                                <  4:152443/117‑1 (MQ=37)
                                                                                                                                   cAACAGGTCGCGGCCATCTGG‑CGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCt                                  <  7:164505/131‑1 (MQ=255)
                                                                                                                                   cAACAGGTCGCGGCCATCTGG‑CGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCt                                  >  8:164505/1‑131 (MQ=255)
                                                                                                                                    aaCACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGtt                <  4:151666/148‑1 (MQ=255)
                                                                                                                                    aaCACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTg               >  2:294309/1‑149 (MQ=255)
                                                                                                                                       acGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTgcgc            <  2:240094/149‑1 (MQ=255)
                                                                                                                                            gc‑gcCAGATGGCGGCGACCTGTTGGTGCAGCGCGGTGCGTG‑‑‑GCTCGCCTTGCTCTCGGTG‑GTGCCACACTGGGGCTGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCa       <  2:289846/149‑1 (MQ=255)
                                                                                                                                                 cAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCAt   >  1:341388/1‑149 (MQ=255)
                                                                                                                                                  aGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATg  >  7:3573/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATG  >  NC_002947/4826136‑4826429

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: ‑