Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,572,249 T→A 5.6% L124* (TTG→TAG) ‡ 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,2490TA 5.6% 218.8 / 2.8 71L124* (TTG→TAG) ‡pcaB3‑carboxy‑cis,cis‑muconate cycloisomerase
Reads supporting (aligned to +/- strand):  ref base T (28/39);  new base A (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 = 7.19e-01

GCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTGCTT  >  NC_002947/1572102‑1572395
                                                                                                                                                   |                                                                                                                                                   
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                                                                                              <  1:313328/149‑1 (MQ=255)
                                                     gTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGaa                                                                                              <  4:282002/149‑1 (MQ=255)
                                                     gTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGaa                                                                                              <  3:389396/149‑1 (MQ=255)
                                                            agGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCt                                                                                                                                      <  4:341721/103‑1 (MQ=255)
                                                            agGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCt                                                                                                                                      >  3:341721/1‑103 (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                                                                                                       >  6:81366/1‑111 (MQ=255)
                                                                                  tGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCgcagca                                                                                                       <  5:81366/111‑1 (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)
                                                                                                                                                   |                                                                                                                                                   
GCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGTGCTT  >  NC_002947/1572102‑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: