Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,002,525 G→C 6.4% G189A (GGT→GCT) ‡ PP_2626 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,002,5250GC 6.4% 198.0 / 3.3 62G189A (GGT→GCT) ‡PP_2626hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base C (2/2);  total (28/34)
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.65e-01

ATCGTTTTACCTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATG  >  NC_002947/3002381‑3002661
                                                                                                                                                 |                                                                                                                                        
aTCGTTTTACCTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGtt                                                                                                                                      >  6:347295/1‑149 (MQ=255)
  cGTTTTACCTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTgg                                                                                                                                    >  1:15960/1‑149 (MQ=255)
   gTTTTACCTACACCCCGACCCACGCAGCGGCGAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGACAACGGTGTTgg                                                                                                                                    >  1:39596/1‑149 (MQ=255)
      ttACCTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGa                                                                                                                                 <  6:375102/148‑1 (MQ=255)
      ttACCTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGa                                                                                                                                 <  3:215358/148‑1 (MQ=255)
         ccTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGgct                                                                                                                             <  7:44440/149‑1 (MQ=255)
          cTACACCCCGACCCACGCGGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTCGCCAACACCGTTGGCAAgttcc                                                                                                                            >  3:89058/1‑145 (MQ=255)
          cTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGgctg                                                                                                                            >  1:27914/1‑149 (MQ=255)
             cacCCCGACCCACGCGGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTCGCCAACACCGTTGGCAAgttccga                                                                                                                          <  4:89058/148‑7 (MQ=255)
               ccccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTg                                                                                                                                     >  4:193391/1‑135 (MQ=255)
               ccccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTg                                                                                                                                     <  3:193391/135‑1 (MQ=255)
               ccccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCtgtg                                                                                                                       <  4:50620/149‑1 (MQ=255)
                cccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGCCACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGc                                                                                                                      >  4:221556/1‑149 (MQ=255)
                cccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGc                                                                                                                      <  3:178105/149‑1 (MQ=255)
                cccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTCGCCAACACCGTTGGCAAgttccga                                                                                                                          <  8:90709/145‑7 (MQ=255)
                cccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTCGCCAACACCGTTGGCAAgttccga                                                                                                                          >  7:90709/1‑139 (MQ=255)
                 ccGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGcc                                                                                                                     >  6:12339/1‑149 (MQ=255)
                       cACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAgc                                                                                                               >  8:13013/1‑149 (MQ=255)
                             gcggc‑agcgCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGc                                                                                                         <  1:297428/149‑1 (MQ=255)
                               ggc‑agcgCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTg                                                                                                       >  8:328746/1‑149 (MQ=255)
                                   agcgCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCaa                                                                                                    <  7:13013/149‑1 (MQ=255)
                                      gcAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAAcc                                                                                                  >  3:316119/1‑148 (MQ=255)
                                       cAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCga                                                                                                <  2:27914/149‑1 (MQ=255)
                                       cAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCga                                                                                                <  4:338996/149‑1 (MQ=255)
                                          gCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCTGCTGCAACCgag                                                                                               >  4:53612/1‑147 (MQ=255)
                                              ccTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCt                                                                                         <  5:12339/149‑1 (MQ=255)
                                                ttGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTg                                                                                        >  5:230927/1‑148 (MQ=255)
                                                   aaaCGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCAcc                                                                                    >  3:29835/1‑149 (MQ=255)
                                                     aCGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCAGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCgg                                                                                                          <  2:50411/125‑1 (MQ=255)
                                                     aCGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCAGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCgg                                                                                                          >  1:50411/1‑125 (MQ=255)
                                                       gCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTg                                                                                <  4:29835/149‑1 (MQ=255)
                                                       gCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTg                                                                                <  2:15960/149‑1 (MQ=255)
                                                       ccgatctacaacaGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTg                                                                                <  5:247821/142‑1 (MQ=255)
                                                         cTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGc                                                                              <  7:328746/149‑1 (MQ=255)
                                                             aacaacAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGAt                                                                          <  2:341119/149‑1 (MQ=255)
                                                               caacaGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTg                                                                                >  6:247821/1‑141 (MQ=255)
                                                                  caGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCAt                                                                     <  5:347295/149‑1 (MQ=255)
                                                                  caGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACCTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGTCGATGTCAt                                                                     >  8:108518/1‑149 (MQ=255)
                                                                     gAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGc                                                                  >  3:215621/1‑149 (MQ=255)
                                                                          aGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCacac                                                             >  5:1068/1‑149 (MQ=255)
                                                                              cAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCgt                                                         <  2:362623/149‑1 (MQ=255)
                                                                               aGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGCAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCACCCGAGAAGCCTGCACCGTTCGCCGATGTCATCGCCCCACCCgag                                                        >  8:339346/1‑147 (MQ=255)
                                                                                 cTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACCTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGTCGATGTCATCGCCACACCCGTGTc                                                      <  7:108518/149‑1 (MQ=255)
                                                                                    gACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCAc                                                    <  8:70799/148‑1 (MQ=255)
                                                                                       acGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGtt                                                                                 >  2:264627/1‑116 (MQ=255)
                                                                                       acGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGtt                                                                                 <  1:264627/116‑1 (MQ=255)
                                                                                         ggCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATgc                                              <  3:221556/149‑2 (MQ=255)
                                                                                         ggCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGt                                              >  5:318556/1‑149 (MQ=255)
                                                                                             gcaCTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGAt                                          <  6:318556/149‑1 (MQ=255)
                                                                                                cTGGGCGAGGTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATggg                                       <  7:339346/149‑1 (MQ=255)
                                                                                                          ttGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGt                             >  7:310405/1‑149 (MQ=255)
                                                                                                            gAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTg                            <  6:1068/148‑1 (MQ=255)
                                                                                                             aGCGCCTGCTAGGCAAGTACTCGGAACTTGGCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTga                           <  4:198453/148‑1 (MQ=255)
                                                                                                                  cTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACg                     >  1:269881/1‑149 (MQ=255)
                                                                                                                        ggCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGACGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCgg                 <  3:53612/147‑1 (MQ=255)
                                                                                                                            aGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCTATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTg           <  4:215621/149‑1 (MQ=255)
                                                                                                                            aGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTg           <  5:322894/149‑1 (MQ=255)
                                                                                                                                 tCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGt       >  7:309192/1‑148 (MQ=255)
                                                                                                                                   ggAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGt       <  8:309192/146‑1 (MQ=255)
                                                                                                                                   ggAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCa    >  6:273627/1‑149 (MQ=255)
                                                                                                                                     aaCTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATg  >  7:321917/1‑149 (MQ=255)
                                                                                                                                     aaCTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATg  <  4:316119/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
ATCGTTTTACCTACACCCCGACCCACGCAGCGGC‑AGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCTGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATG  >  NC_002947/3002381‑3002661

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