Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,718,979 C→G 7.2% G141A (GGC→GCC)  paaF ← enoyl‑CoA hydratase/isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,718,9790CG 7.2% 211.2 / 6.8 69G141A (GGC→GCC) paaFenoyl‑CoA hydratase/isomerase
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (3/2);  total (36/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 = 8.65e-01

TGGCGCGTTCTACGGTGAGTTCGGGCTGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGTTTGGCGAAT  >  NC_002947/3718833‑3719125
                                                                                                                                                  |                                                                                                                                                  
tGGCGCGTTCTACGGTGAGTTCGGGCTGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCACCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCa                                                                                                                                                  >  2:308196/1‑149 (MQ=255)
          tACGGTGAGTTCGGGCTGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGga                                                                                                                                        >  3:270342/1‑149 (MQ=255)
            cGGTGAGTTCGGGCTGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGgat                                                                                                                                       >  6:43257/1‑148 (MQ=255)
             ggTGAGTTCGGGCTGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGgatga                                                                                                                                     >  4:187172/1‑149 (MQ=255)
                       ggCTGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCACCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGTGATGATGCCAAGGtt                                                                                                                           <  1:308196/149‑1 (MQ=255)
                          tGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGAt                                                                                                                        >  8:300975/1‑149 (MQ=255)
                           gggTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGAtt                                                                                                                       >  7:203108/1‑149 (MQ=255)
                            ggTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGAttt                                                                                                                      >  6:27744/1‑149 (MQ=255)
                            ggTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTGGAttt                                                                                                                      >  5:389663/1‑149 (MQ=255)
                             gtcaattCGCTGACCAGGCCCGCCAGCTGCGCACCGCGGCCAGTGAAGGCTTCCCCCGTCAGCACCATAAGCATGGCCAAACGCTTACCGACGGCACGAAACAGCCGATGGCTCCCGCCAGCAACCGGCATGATGACAAGGATGATTTc                                                                                                                     >  4:202531/6‑149 (MQ=255)
                              tGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTcc                                                                                                                    <  5:387817/149‑1 (MQ=255)
                              tGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTcc                                                                                                                    <  3:66776/149‑1 (MQ=255)
                              tGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTcc                                                                                                                    <  3:72060/149‑1 (MQ=255)
                                           cAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAAc                                                                                                       <  3:180864/149‑1 (MQ=255)
                                           cAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAAc                                                                                                       <  7:294512/149‑1 (MQ=255)
                                           cAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAAc                                                                                                       >  4:76728/1‑149 (MQ=255)
                                                   ccTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCg                                                                                               <  1:211387/149‑1 (MQ=255)
                                                          ggCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCGGGGGGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTgg                                                                                        <  2:381602/149‑1 (MQ=255)
                                                           gcgcGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGc                                                                                       >  2:162915/1‑149 (MQ=255)
                                                             gcgGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGa                                                                                     <  5:43257/149‑1 (MQ=255)
                                                                   ggCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGTTGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGa                                                                               >  5:247555/1‑149 (MQ=255)
                                                                   ggCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGCCGa                                                                               >  7:75359/1‑149 (MQ=255)
                                                                   ggCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGa                                                                               <  4:66167/149‑1 (MQ=255)
                                                                    gCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGGCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGAt                                                                              <  8:157876/149‑1 (MQ=255)
                                                                      aGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGt                                                                            >  3:376654/1‑149 (MQ=255)
                                                                       gTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTc                                                                           >  6:41363/1‑149 (MQ=255)
                                                                        tGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTa                                                                                            >  2:87105/1‑131 (MQ=255)
                                                                        tGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTa                                                                                            <  1:87105/131‑1 (MQ=255)
                                                                          aTGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCgg                                                                         <  8:203108/148‑1 (MQ=255)
                                                                                  ccccGTCAGCACCATCTGCATGGCCAACGGTTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTg                                                                >  2:39607/1‑149 (MQ=255)
                                                                                   cccGTCAGCACCATCTGCATGGCCTACGGCTTAGCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGTCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGc                                                               <  3:202531/149‑1 (MQ=255)
                                                                                    ccGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACgg                                                                                                     >  1:105988/1‑110 (MQ=255)
                                                                                    ccGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACgg                                                                                                     <  2:105988/110‑1 (MQ=255)
                                                                                    ccGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGc                                                               >  2:372177/1‑148 (MQ=255)
                                                                                    ccGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGc                                                               >  5:302257/1‑148 (MQ=255)
                                                                                     cGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGCACTGGCGATGACGATGTCGGCGCACATTGCCa                                                             <  6:234045/149‑1 (MQ=255)
                                                                                         aGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTc                                                         <  6:343391/149‑1 (MQ=255)
                                                                                           cACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCg                                                        >  2:344020/1‑148 (MQ=255)
                                                                                              cATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAAc                                                    >  4:81123/1‑149 (MQ=255)
                                                                                                   gCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACAtt                                                                 >  4:328461/1‑131 (MQ=255)
                                                                                                   gCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACAtt                                                                 <  3:328461/131‑1 (MQ=255)
                                                                                                   gCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCa                                               >  6:46756/1‑149 (MQ=255)
                                                                                                    cATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCAc                                              >  3:324912/1‑149 (MQ=255)
                                                                                                         ccAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAgg                                         >  5:18704/1‑149 (MQ=255)
                                                                                                           aaCGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGATGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAggg                                        >  7:199953/1‑148 (MQ=255)
                                                                                                             cGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCAt                                     <  3:330767/149‑1 (MQ=255)
                                                                                                              ggCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATa                                    <  6:18704/149‑1 (MQ=255)
                                                                                                              cgcgtACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCCAGGTTGATTTCCGGCTGGCCGAACCGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATa                                    <  7:374214/145‑1 (MQ=255)
                                                                                                                cTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATa                                    <  8:201363/147‑1 (MQ=255)
                                                                                                                cTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATa                                    >  7:201363/1‑147 (MQ=255)
                                                                                                                 ttACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGcc                                 >  7:365138/1‑149 (MQ=255)
                                                                                                                   aCCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGcc                                 <  6:389663/147‑1 (MQ=255)
                                                                                                                   aCCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCg                                <  7:300975/148‑1 (MQ=255)
                                                                                                                   aCCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGt                               >  4:219195/1‑149 (MQ=255)
                                                                                                                    ccGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGtt                              <  3:234433/149‑1 (MQ=255)
                                                                                                                      gACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGa                            <  5:27744/149‑1 (MQ=255)
                                                                                                                         ggCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTg                         <  1:372177/149‑1 (MQ=255)
                                                                                                                          gCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTg                         >  5:109809/1‑148 (MQ=255)
                                                                                                                           cACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTgca                       <  1:162915/149‑1 (MQ=255)
                                                                                                                           cACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTgca                       <  5:46756/149‑1 (MQ=255)
                                                                                                                              gTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAACCGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCCCCCAGGGCATAGCCGTTGACTgcagc                     >  8:329087/1‑148 (MQ=255)
                                                                                                                                aaCAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAAt                  >  6:15992/1‑149 (MQ=255)
                                                                                                                                        cccggTGCTGGCGGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGAt                                                                              >  4:320791/4‑81 (MQ=255)
                                                                                                                                        cccggTGCTGGCGGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGAt                                                                              <  3:320791/78‑1 (MQ=255)
                                                                                                                                         ccggTGCTGGCGGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGa                            >  6:131259/3‑130 (MQ=255)
                                                                                                                                         ccggTGCTGGCGGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGa                            <  5:131259/128‑1 (MQ=255)
                                                                                                                                         ccggTGCTGGCGGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGACAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGt           >  3:124658/3‑147 (MQ=255)
                                                                                                                                           ggTGCCGCCAGCAGCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGttt         <  4:270342/147‑1 (MQ=255)
                                                                                                                                           ggTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGTTTg        <  1:216607/148‑1 (MQ=255)
                                                                                                                                               ccgccAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGTTTGGCGaa   <  5:41363/149‑1 (MQ=255)
                                                                                                                                                cgccAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGTTTGGCGAAt  >  5:104869/1‑149 (MQ=255)
                                                                                                                                                cgccAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGTTTGGCGAAt  <  3:76728/149‑1 (MQ=255)
                                                                                                                                                  ccAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTgcag                      <  2:347336/127‑1 (MQ=255)
                                                                                                                                                  ccAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTgcag                      >  1:347336/1‑127 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGGCGCGTTCTACGGTGAGTTCGGGCTGGGTGATTTCGCTGACCAGGCCGGCCTGCTGGGCGCGGAGGGCAGTGATGGCTTCCCCCGTCAGCACCATCTGCATGGCCAACGGCTTACCGACGGCACGTAACAGGCGCTGGGTGCCGCCAGCACCGGGGATGATGCCAAGGTTGATTTCCGGCTGGCCGAAACGGGCGTCGGTACTGGCGATGACGATGTCGGCGCACATTGCCAGTTCGCAACCGCCACCCAGGGCATAGCCGTTGACTGCAGCAATCAGCGGTTTGGCGAAT  >  NC_002947/3718833‑3719125

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