Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,238,061 G→A 6.1% A78V (GCG→GTG)  catB ← muconate cycloisomerase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,238,0610GA 6.1% 187.0 / 2.0 65A78V (GCG→GTG) catBmuconate cycloisomerase 1
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base A (2/2);  total (28/37)
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.79e-01

CAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGGGTCTGCTGCTGCATG  >  NC_002947/4237913‑4238210
                                                                                                                                                    |                                                                                                                                                     
cagTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGc                                                                                                                                                      <  2:709022/150‑1 (MQ=255)
  ggtCGCTGGCCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTg                                                                                                                                                    <  2:1272406/148‑1 (MQ=255)
        tGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGa                                                                                                                                              <  1:1054889/150‑1 (MQ=255)
        tGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGa                                                                                                                                              <  1:200817/150‑1 (MQ=255)
         gACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAg                                                                                                                                             <  2:79980/150‑1 (MQ=255)
         gACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCGTGCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCACCTTGTCCAGCTTGACCATGGCCGCCTTCATGTTGTCTACCCTCACCCCAATCAACGCTGTGGCGAg                                                                                                                                             <  2:1165002/150‑1 (MQ=255)
          aCCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTACCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGt                                                                                                                                           >  1:659397/1‑150 (MQ=255)
             ggCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTAc                                                                                                                                         <  2:1890505/150‑1 (MQ=255)
                     cAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGAt                                                                                                                                 >  2:1245437/1‑150 (MQ=255)
                      aGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATg                                                                                                                                <  1:1712889/150‑1 (MQ=255)
                        gCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGtt                                                                                                                              <  1:1631448/150‑1 (MQ=255)
                         cGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTg                                                                                                                             <  1:1073199/150‑1 (MQ=255)
                                ccTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAGCGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGa                                                                                                                      <  1:294864/150‑1 (MQ=255)
                                ccTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGa                                                                                                                      <  1:1666296/150‑1 (MQ=255)
                                ccTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGa                                                                                                                      >  2:1082793/1‑150 (MQ=255)
                                         tgAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCTTTCgg                                                                                                             <  2:479504/149‑1 (MQ=255)
                                            gCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTg                                                                                                                       <  1:1198637/137‑1 (MQ=255)
                                            gCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTg                                                                                                                       >  2:1198637/1‑137 (MQ=255)
                                              aaGGCGCTTTCTATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCtt                                                                                                        <  1:662974/150‑1 (MQ=255)
                                                ggCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTc                                                                                                      >  1:1757784/1‑150 (MQ=255)
                                                   gcTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGATTTCGTa                                                                                                   <  1:211815/150‑1 (MQ=255)
                                                     tttCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGAt                                                                                                                     >  2:447587/1‑130 (MQ=255)
                                                     tttCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGAt                                                                                                                     <  1:447587/130‑1 (MQ=255)
                                                           tGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCACCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACg                                                                                           >  2:1265004/1‑150 (MQ=255)
                                                               ggACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTCCGGGGCTTTCGTAGCCATACGCCAg                                                                                       <  1:1245437/150‑1 (MQ=255)
                                                                    tGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCAc                                                                                  <  1:1393736/150‑1 (MQ=255)
                                                                        gAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAt                                                                              >  1:27072/1‑150 (MQ=255)
                                                                        gAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAt                                                                              >  1:1003716/1‑150 (MQ=255)
                                                                        gAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAt                                                                              >  1:911956/1‑150 (MQ=255)
                                                                              gtTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtggtggt                                                                        <  2:27072/150‑1 (MQ=255)
                                                                                tGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtggtggtgg                                                                      <  1:1082793/150‑1 (MQ=255)
                                                                                   ccTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGa                                                                               <  1:1205511/138‑1 (MQ=255)
                                                                                   ccTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGa                                                                               >  2:1205511/1‑138 (MQ=255)
                                                                                   ccTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCtt                                                                   >  2:1447897/1‑150 (MQ=255)
                                                                                         ccAGCTTGTCGAGCTTGAGCATGGCGGCATTGATGTTGTCTGCGCGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGa                                                             <  2:1297543/150‑1 (MQ=255)
                                                                                         ccAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGa                                                             <  2:1298316/150‑1 (MQ=255)
                                                                                          cAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtggtg                                                                          >  1:477282/1‑136 (MQ=255)
                                                                                          cAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtggtg                                                                          <  2:477282/136‑1 (MQ=255)
                                                                                          cAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGAt                                                            >  1:550601/1‑150 (MQ=255)
                                                                                          cAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGAt                                                            >  2:70110/1‑150 (MQ=255)
                                                                                            gCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTgg                                                                                                                            <  2:765424/84‑1 (MQ=255)
                                                                                            gCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTgg                                                                                                                            >  1:765424/1‑84 (MQ=255)
                                                                                              ttGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCt                                                        >  2:243723/1‑150 (MQ=255)
                                                                                                 tCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGc                                                                                          <  1:1798170/113‑1 (MQ=255)
                                                                                                 tCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGc                                                                                          >  2:1798170/1‑113 (MQ=255)
                                                                                                 tCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACACTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTcc                                                     >  1:274368/1‑150 (MQ=255)
                                                                                                 tCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACACTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTcc                                                     >  1:274302/1‑150 (MQ=255)
                                                                                                   cAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCAc                                                   >  1:1264758/1‑150 (MQ=255)
                                                                                                     gCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCggg                                                                                                            <  2:304880/91‑1 (MQ=255)
                                                                                                     gCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCggg                                                                                                            >  1:304880/1‑91 (MQ=255)
                                                                                                     gCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGAAGCTTCAACGATTCCTTCCACGc                                                 >  2:107092/1‑150 (MQ=255)
                                                                                                      cTTGAGCATGGCGGCATTGATGTTGTCTGCAGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGcc                                                <  2:911956/150‑1 (MQ=255)
                                                                                                       ttGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCa                                               <  1:1675525/150‑1 (MQ=255)
                                                                                                        tGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCAt                                              <  1:1265004/150‑1 (MQ=255)
                                                                                                          aGCAGGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCg                                            <  2:1435072/150‑1 (MQ=255)
                                                                                                              tGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAACCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGc                                        <  2:985233/150‑1 (MQ=255)
                                                                                                                          tgttgtCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCa                            <  2:181440/150‑1 (MQ=255)
                                                                                                                            ttgtCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAAt                          >  1:775323/1‑150 (MQ=255)
                                                                                                                                 tGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATaccac                     >  2:1868805/1‑150 (MQ=255)
                                                                                                                                    cGGCAGGCCAATCAACACTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCGGTCGCAATACCACCAg                  <  2:274302/150‑1 (MQ=255)
                                                                                                                                    cGGCAGGCCAATCAACACTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAg                  <  2:274368/150‑1 (MQ=255)
                                                                                                                                     ggCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAgg                 <  2:1765648/150‑1 (MQ=255)
                                                                                                                                          gCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGCTGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGGGTctg            >  2:853805/1‑150 (MQ=255)
                                                                                                                                              aTCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGGGTctgctgc        >  1:322380/1‑150 (MQ=255)
                                                                                                                                                 aaCGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGGGTCtgctgctgc     >  1:532176/1‑150 (MQ=255)
                                                                                                                                                    gCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGGGTCTGCTGCTGCATg  <  2:1003716/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTACGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGGGTCTGCTGCTGCATG  >  NC_002947/4237913‑4238210

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