Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,572,248 T→C 5.6% L124L (TTG→CTG) ‡ pcaB → 3‑carboxy‑cis,cis‑muconate cycloisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,572,2480TC 5.6% 212.9 / 2.6 71L124L (TTG→CTG) ‡pcaB3‑carboxy‑cis,cis‑muconate cycloisomerase
Reads supporting (aligned to +/- strand):  ref base T (28/39);  new base C (2/2);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

GGGCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTGCTT  >  NC_002947/1572100‑1572395
                                                                                                                                                    |                                                                                                                                                    
gggCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGAtt                                                                                                                                                      <  2:282134/149‑1 (MQ=255)
  gCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGAttt                                                                                                                                                     >  2:315171/1‑148 (MQ=255)
  gCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTg                                                                                                                                                    <  4:65794/149‑1 (MQ=255)
  gCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATCTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGAttt                                                                                                                                                     >  5:29051/1‑148 (MQ=255)
    aaCTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCCATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGAt                                                                                                                                                  >  3:116604/1‑149 (MQ=255)
      cTCGGCCATTCCGCTCGTGAAAGCGTTGGCCAAGGTGATCGCCACCGGCGTGCCAGAGTCTGAGCGCTATGTGCAACTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATc                                                                                                                                                 >  4:249485/1‑148 (MQ=255)
       tCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATc                                                                                                                                                 >  1:72962/1‑147 (MQ=255)
         ggCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAg                                                                                                                                              <  1:305648/148‑1 (MQ=255)
         ggCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATCAAATCGAgg                                                                                                                                             <  3:17667/149‑1 (MQ=255)
          gCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGc                                                                                                                                            >  5:82552/1‑149 (MQ=255)
           ccATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGcc                                                                                                                                           >  8:17064/1‑149 (MQ=255)
            ctttCCGCTGGTGAATGCTTTTGTCAAGGTGAGAGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCg                                                                                                                                          <  6:29051/147‑1 (MQ=255)
            cATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGcc                                                                                                                                           >  8:305215/1‑148 (MQ=255)
              ttCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGAc                                                                                                                                        >  5:104731/1‑149 (MQ=255)
              ttCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATCAAATCGAgggcatc                                                                                                                                        <  3:98719/149‑6 (MQ=255)
                  gCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCg                                                                                                                                    <  1:234994/149‑1 (MQ=255)
                        gAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGc                                                                                                                             >  8:295963/1‑149 (MQ=255)
                         aaaGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGc                                                                                                                             <  2:77866/148‑1 (MQ=255)
                          aaGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTg                                                                                                                           <  6:51812/149‑1 (MQ=255)
                             cGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGcc                                                                                                                        <  6:104731/149‑1 (MQ=255)
                                  ggCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAc                                                                                                                   >  1:147533/1‑149 (MQ=255)
                                  ggCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCCATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAc                                                                                                                   <  4:116604/149‑1 (MQ=255)
                                     aaGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCt                                                                                                                <  7:295963/149‑1 (MQ=255)
                                      aGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTg                                                                                                               >  8:55895/1‑149 (MQ=255)
                                       ggTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGt                                                                                                              <  2:72962/149‑1 (MQ=255)
                                        gTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTc                                                                                                             <  7:17064/149‑1 (MQ=255)
                                              gcgaCCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCgcagc                                                                                                        <  8:347313/145‑1 (MQ=255)
                                                cACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAgg                                                                                                     <  7:55895/149‑1 (MQ=255)
                                                cACCCGCGTGCCACCCGCTGACCGCTATGTTGACCTTGGGGGCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAgg                                                                                                     <  4:285931/149‑1 (MQ=255)
                                                  ccGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGcc                                                                                                   >  2:68152/1‑149 (MQ=255)
                                                     gCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCtt                                                                                                 >  5:74898/1‑148 (MQ=255)
                                                       gTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGaa                                                                                              <  3:389396/149‑1 (MQ=255)
                                                       gTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGaa                                                                                              <  1:313328/149‑1 (MQ=255)
                                                       gTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGaa                                                                                              <  4:282002/149‑1 (MQ=255)
                                                              agGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCt                                                                                                                                      >  3:341721/1‑103 (MQ=255)
                                                              agGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCt                                                                                                                                      <  4:341721/103‑1 (MQ=255)
                                                                 cTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGAt                                                                                    <  2:147533/149‑1 (MQ=255)
                                                                    aGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATa                                                                                   <  3:249485/147‑1 (MQ=255)
                                                                            gtgCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTtggt                                                                         >  7:155099/1‑149 (MQ=255)
                                                                             tgCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTtggt                                                                         <  7:305215/148‑1 (MQ=255)
                                                                              gCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTtggtgg                                                                       >  6:212974/1‑149 (MQ=255)
                                                                               caacTTGGGGCCAACATCCATGTTGCGTAGGTAACCGGTCTTGTTCTGGATGTTCGGGATGGCCTTGTTTTTATCGATGCCGTACTCGGGACTGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTtggtg                                                                        <  1:315171/145‑1 (MQ=255)
                                                                               cACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTggg                                                                      <  2:134509/149‑1 (MQ=255)
                                                                                aCCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGt                                                                     >  5:149045/1‑149 (MQ=255)
                                                                                aCCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGt                                                                     >  1:275080/1‑149 (MQ=255)
                                                                                 ccTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTc                                                                    <  2:275080/149‑1 (MQ=255)
                                                                                  cTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAc                                                                                                                   <  4:382909/101‑1 (MQ=255)
                                                                                  cTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAc                                                                                                                   >  3:382909/1‑101 (MQ=255)
                                                                                   ttGGGGCAACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCg                                                                   >  5:290755/1‑148 (MQ=255)
                                                                                    tGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCgcagca                                                                                                       <  5:81366/111‑1 (MQ=255)
                                                                                    tGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCgcagca                                                                                                       >  6:81366/1‑111 (MQ=255)
                                                                                         ccaccaGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTgg                                                            <  3:313785/149‑1 (MQ=255)
                                                                                           accaGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGtt                                                          <  3:8878/149‑1 (MQ=255)
                                                                                                ccagGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCaa                                                      >  6:371912/1‑148 (MQ=255)
                                                                                                   ggATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACg                                                  <  6:74898/149‑1 (MQ=255)
                                                                                                    gATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGc                                                 >  4:171883/1‑149 (MQ=255)
                                                                                                     aTGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGcc                                                <  6:82552/149‑1 (MQ=255)
                                                                                                     aTGCGATGGATACCGGGCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGcc                                                <  5:212974/149‑1 (MQ=255)
                                                                                                                 ccGGTCTTGTTGTGCAGCTGCGCGATGCGCTCGATTTGATCGAGGCCGAACTCGGGAAGCTGGCCGATACCCTGTGGCAGCAGGCCTTGAAGCACGCCGATACGCGCTTGGTGGGTCGTACCTGGTTGCAACACGCGACCCCGGTGAcc                                    <  1:190377/149‑1 (MQ=255)
                                                                                                                 ccGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGAcc                                    <  5:137635/149‑1 (MQ=255)
                                                                                                                  cGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGAccc                                   >  2:409218/1‑149 (MQ=255)
                                                                                                                     tCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTgg                                <  6:149045/149‑1 (MQ=255)
                                                                                                                        ggTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCa                             >  6:174048/1‑149 (MQ=255)
                                                                                                                              gCAGCTGCGCGATGCCCTCGATTTCATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCAGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGCTGGGTCGTACCTGGTTCCAACACGCCACCCCGCTGACCCTGGGCATGaaa                        >  4:61260/1‑148 (MQ=255)
                                                                                                                               cAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGaaa                        <  5:276229/147‑1 (MQ=255)
                                                                                                                               cAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAAc                       <  3:171883/148‑1 (MQ=255)
                                                                                                                                aGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTg                     <  5:371912/149‑1 (MQ=255)
                                                                                                                                      gcgATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGgt               <  1:409218/149‑1 (MQ=255)
                                                                                                                                          cggccTCGATCAAATCGAGGGCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGAGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTAc           >  5:140306/4‑149 (MQ=255)
                                                                                                                                             ccTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTgg        >  7:104272/1‑149 (MQ=255)
                                                                                                                                             ccTCGATCAAATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGACCGGTGTACTgg        >  7:301884/1‑149 (MQ=255)
                                                                                                                                               tCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGt      >  4:87196/1‑149 (MQ=255)
                                                                                                                                               tCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGt      >  3:39345/1‑149 (MQ=255)
                                                                                                                                                 gATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGt      <  6:290755/147‑1 (MQ=255)
                                                                                                                                                    ttGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTGCtt  <  8:104272/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTGCTT  >  NC_002947/1572100‑1572395

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: