Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,237,903 T→A 6.5% Q116L (CAG→CTG)  catB ← muconate cycloisomerase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,237,9030TA 6.5% 183.3 / 2.8 61Q116L (CAG→CTG) catBmuconate cycloisomerase 1
Reads supporting (aligned to +/- strand):  ref base T (32/25);  new base A (2/2);  total (34/27)
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.93e-01

ACGCGGTGCCGGCGAATGTCCAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGA  >  NC_002947/4237756‑4238050
                                                                                                                                                   |                                                                                                                                                   
acGCGGTGCCGGCGAATGTCCAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTg                                                                                                                                                    >  8:86628/1‑149 (MQ=255)
 cgcgGTGCCGGCGAATGTCCAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTgg                                                                                                                                                   <  1:164104/149‑1 (MQ=255)
 cgcgGTGCCGGCGAATGTCCAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTgg                                                                                                                                                   <  3:85315/149‑1 (MQ=255)
      tGCCGGCGAATGTCCAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTc                                                                                                                                              <  4:166899/149‑1 (MQ=255)
           gCGAATGTCCAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGc                                                                                                                                          >  4:315752/1‑148 (MQ=255)
                    cAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCtt                                                                                                                                 >  2:289722/1‑148 (MQ=255)
                       cATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCtt                                                                                                                                 >  6:221976/1‑145 (MQ=255)
                       cATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCtt                                                                                                                                 <  5:221976/145‑1 (MQ=255)
                           tgctgtgctTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGc                                                                                                                          <  1:289722/148‑1 (MQ=255)
                                tgctTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCttt                                                                                                                                >  7:132598/1‑137 (MQ=255)
                                tgctTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCttt                                                                                                                                <  8:132598/137‑1 (MQ=255)
                                 gctTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACtt                                                                                                                   >  2:267153/1‑149 (MQ=255)
                                  ctTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTg                                                                                                                  >  3:100992/1‑149 (MQ=255)
                                  ctTCGGCGATGTCGCGGGCGGTGGCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTg                                                                                                                  <  4:178064/149‑1 (MQ=255)
                                   ttCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTgg                                                                                                                 >  8:164958/1‑149 (MQ=255)
                                   ttCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTgg                                                                                                                 >  6:75501/1‑149 (MQ=255)
                                     cGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGc                                                                                                                <  4:100992/148‑1 (MQ=255)
                                      ggCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTgg                                                                                                                                                   >  2:180726/1‑112 (MQ=255)
                                      ggCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTgg                                                                                                                                                   <  1:180726/112‑1 (MQ=255)
                                               gcgggcggTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTg                                                                                                      <  5:75501/148‑1 (MQ=255)
                                                 gggcggTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGccc                                                                                                   >  7:3989/1‑149 (MQ=255)
                                                    cggTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTg                                                                                                <  3:312222/149‑1 (MQ=255)
                                                           gccgcTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGc                                                                                          <  7:30838/148‑1 (MQ=255)
                                                            ccgcTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCtt                                                                                        >  1:222986/1‑149 (MQ=255)
                                                               cTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTc                                                                                     >  7:156179/1‑149 (MQ=255)
                                                                    cAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCt                                                                                >  6:144876/1‑149 (MQ=255)
                                                                       ggTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCtt                                                                               >  2:277204/1‑147 (MQ=255)
                                                                       ggTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGa                                                                             <  7:86628/149‑1 (MQ=255)
                                                                       ggTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGa                                                                             <  1:277204/149‑1 (MQ=255)
                                                                       ggTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGa                                                                             <  2:222986/149‑1 (MQ=255)
                                                                       ggTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGa                                                                             <  6:117528/149‑1 (MQ=255)
                                                                             ggCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATgg                                                                       >  1:315459/1‑149 (MQ=255)
                                                                             ggCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATgg                                                                       >  1:310623/1‑149 (MQ=255)
                                                                              gCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGc                                                                                          >  5:325636/1‑129 (MQ=255)
                                                                              gCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGc                                                                                          <  6:325636/129‑1 (MQ=255)
                                                                                cACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATggcgg                                                                    <  7:266216/149‑1 (MQ=255)
                                                                                   ttCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCAt                                                                 >  5:129377/1‑149 (MQ=255)
                                                                                       aGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGAt                                                             >  3:23851/1‑149 (MQ=255)
                                                                                       aGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGAt                                                             <  4:200356/149‑1 (MQ=255)
                                                                                       aGGCTGTCACGCACGCGGCAACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTc                                                                                     >  8:67909/1‑125 (MQ=255)
                                                                                       aGGCTGTCACGCACGCGGCAACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTc                                                                                     <  7:67909/125‑1 (MQ=255)
                                                                                                                   gTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCt                                 >  6:198032/1‑149 (MQ=255)
                                                                                                                    ttCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGc                                  <  8:78152/147‑1 (MQ=255)
                                                                                                                    ttCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGc                                  >  7:78152/1‑147 (MQ=255)
                                                                                                                    ttCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTg                                <  8:156179/149‑1 (MQ=255)
                                                                                                                    ttCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTg                                >  7:9983/1‑149 (MQ=255)
                                                                                                                       gCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTgggg                             >  5:7921/1‑149 (MQ=255)
                                                                                                                       gCTGACCGGCAGGCCCAGGCGTTTGCCCAGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTg                                                     <  7:18643/125‑1 (MQ=255)
                                                                                                                       gCTGACCGGCAGGCCCAGGCGTTTGCCCAGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTg                                                     <  7:13811/125‑1 (MQ=255)
                                                                                                                       gCTGACCGGCAGGCCCAGGCGTTTGCCCAGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTg                                                     >  8:13811/1‑125 (MQ=255)
                                                                                                                       gCTGACCGGCAGGCCCAGGCGTTTGCCCAGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTg                                                     >  8:18643/1‑125 (MQ=255)
                                                                                                                           ggtggCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAg                         <  6:152782/146‑1 (MQ=255)
                                                                                                                           aCCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGTGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAg                         <  3:315752/149‑1 (MQ=255)
                                                                                                                             cGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAgg                        >  3:62058/1‑148 (MQ=255)
                                                                                                                             cGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGAAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGAGGCATTGATGTTGTCTGCCGGCAGGCCAAGCAAAGGTGTGGCGAgg                        >  8:210980/1‑148 (MQ=255)
                                                                                                                                   gCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGt                 >  4:215113/1‑149 (MQ=255)
                                                                                                                                    cccAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTc                >  1:138742/1‑149 (MQ=255)
                                                                                                                                        ggCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATg            >  7:204624/1‑149 (MQ=255)
                                                                                                                                                cccTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCtt    >  2:44324/1‑149 (MQ=255)
                                                                                                                                                 ccTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCtt    >  6:6645/1‑148 (MQ=255)
                                                                                                                                                  cTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGa  <  6:129377/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGCGGTGCCGGCGAATGTCCAGCATGTGCTGTGCTTCGGCGATGTCGCGGGCGGTGTCGCCGCTGGCCAGGGTCCAGGCCACTTCCAGGCTGTCACGCACGCGGCCACCCAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGA  >  NC_002947/4237756‑4238050

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