Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,638,710 C→A 7.1% H221Q (CAC→CAA)  PP_3205 → fumarylacetoacetate hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,638,7100CA 7.1% 314.8 / 8.9 98H221Q (CAC→CAA) PP_3205fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (51/40);  new base A (5/2);  total (56/42)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAG  >  NC_002947/3638561‑3638859
                                                                                                                                                     |                                                                                                                                                     
gTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTTGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCAGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCAc                                                                                                                                                       <  2:716304/150‑1 (MQ=255)
   aGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTg                                                                                                                                                    >  2:2312656/1‑150 (MQ=255)
     ccTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCa                                                                                                                                                  <  1:2312656/150‑1 (MQ=255)
      cTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCAc                                                                                                                                                 <  1:1566661/150‑1 (MQ=255)
       tGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACt                                                                                                                                                <  2:1072715/150‑1 (MQ=255)
        gCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTg                                                                                                                                               <  1:1512234/150‑1 (MQ=255)
           ggCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGaa                                                                                                                                            >  2:464579/1‑150 (MQ=255)
            gCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAAc                                                                                                                                           <  1:393569/150‑1 (MQ=255)
            gCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAAc                                                                                                                                           >  2:2389041/1‑150 (MQ=255)
             cccTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACg                                                                                                                                          >  1:1563710/1‑150 (MQ=255)
             cccTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACg                                                                                                                                          <  1:279938/150‑1 (MQ=255)
                 ggCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGa                                                                                                                                      <  1:906466/150‑1 (MQ=255)
                 ggCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGa                                                                                                                                      >  1:2183124/1‑150 (MQ=255)
                  gCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAg                                                                                                                                     <  2:2183124/150‑1 (MQ=255)
                  gCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAg                                                                                                                                     >  2:2074931/1‑150 (MQ=255)
                    aCCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTg                                                                                                                                   <  1:247398/150‑1 (MQ=255)
                               tGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCCGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGACGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATc                                                                                                                        >  1:19330/1‑150 (MQ=255)
                                   gACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCa                                                                                                                    >  1:1826988/1‑150 (MQ=255)
                                   gACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGCGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCa                                                                                                                    >  2:1347283/1‑150 (MQ=255)
                                      cggcttcggcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGAGAGTGGTTCGGCCATCCCCATgg                                                                                                                 >  1:1284353/1‑150 (MQ=255)
                                              gcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCa                                                                                                         >  2:886567/1‑150 (MQ=255)
                                                tttCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCCTGGCGgcgcaaag                                                                                                       >  1:646970/1‑146 (MQ=255)
                                                tttCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCACGGAAGACATCGGCGAGGCCTGTTTCGACGTCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                                       >  1:1498832/1‑150 (MQ=255)
                                                   ctcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGCCGCCATgggg                                                                                                    <  1:1915813/150‑4 (MQ=255)
                                                    tcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACt                                                                                                   >  2:238816/1‑150 (MQ=255)
                                                      cAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTACACTGGAACGGCGAg                                                                                                                                     <  2:1946344/114‑1 (MQ=255)
                                                      cAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGGCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTACACTGGAACGGCGAg                                                                                                                                     >  1:1946344/1‑114 (MQ=255)
                                                       aGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCTCTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCg                                                                                                <  2:866872/150‑1 (MQ=255)
                                                       aGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCg                                                                                                <  1:2328219/150‑1 (MQ=255)
                                                         gCCAAGCCCGCCGGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACCCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGt                                                                                              >  1:2065811/1‑150 (MQ=255)
                                                         gCCAAGCCCGCCAGCAGTTTCGCGCACGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGt                                                                                              <  2:265520/150‑1 (MQ=255)
                                                           cAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGttt                                                                                            >  2:2298112/1‑150 (MQ=255)
                                                           cAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGttt                                                                                            >  1:2320046/1‑150 (MQ=255)
                                                           cAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGttt                                                                                            <  1:464579/150‑1 (MQ=255)
                                                           cAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATAGCGGCGCCATGCACTTCGGttt                                                                                            >  2:986915/1‑150 (MQ=255)
                                                           cAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGCACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGttt                                                                                            >  1:2293840/1‑150 (MQ=255)
                                                              gCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTcca                                                                                         <  1:2074931/150‑1 (MQ=255)
                                                               cccGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTg                                                                                                                                               >  1:992090/1‑95 (MQ=255)
                                                               cccGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTg                                                                                                                                               <  2:992090/95‑1 (MQ=255)
                                                                ccGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTccacc                                                                                       >  2:1437261/1‑150 (MQ=255)
                                                                ccGCCAGCAGTTTCGCGCCCGTAGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTccacc                                                                                       >  1:933346/1‑150 (MQ=255)
                                                                 cGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTccacca                                                                                      >  1:1575765/1‑150 (MQ=255)
                                                                    cagcagTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGtt                                                                                   <  2:1445605/150‑1 (MQ=255)
                                                                     agcagTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTg                                                                                  >  2:1251425/1‑150 (MQ=255)
                                                                      gcagTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGa                                                                                 <  1:748442/150‑1 (MQ=255)
                                                                          tttCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGc                                                                             >  2:1677517/1‑150 (MQ=255)
                                                                                gcCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                                       <  1:2097258/118‑1 (MQ=255)
                                                                                gcCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                                       >  2:2097258/1‑118 (MQ=255)
                                                                                gcCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATgc                                                                       <  2:19330/150‑1 (MQ=255)
                                                                                 cccGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCCCCAGTTGATTGCCCATgcg                                                                      >  2:1089098/1‑150 (MQ=255)
                                                                                     ttGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGgcgc                                                                  <  1:2298112/150‑1 (MQ=255)
                                                                                               cccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacgcgacg                                                        >  2:1447800/1‑150 (MQ=255)
                                                                                                ccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCAt                                                                                                        >  2:2363144/1‑101 (MQ=255)
                                                                                                ccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCAt                                                                                                        <  1:2363144/101‑1 (MQ=255)
                                                                                                ccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacgcgacgc                                                       >  2:613450/1‑150 (MQ=255)
                                                                                                ccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacgcgacgc                                                       <  1:1914748/150‑1 (MQ=255)
                                                                                                 cGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgac                                                              >  2:758702/1‑142 (MQ=255)
                                                                                                 cGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgac                                                              <  1:758702/142‑1 (MQ=255)
                                                                                                    aCGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGa                                                   <  1:1447800/150‑1 (MQ=255)
                                                                                                        aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTg                                                    <  1:1597618/145‑1 (MQ=255)
                                                                                                        aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTg                                                    >  2:1597618/1‑145 (MQ=255)
                                                                                                         cTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCg                                              >  1:2166525/1‑150 (MQ=255)
                                                                                                         cTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCg                                              <  2:1575765/150‑1 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGCCGGCCGCGTGCACGTGCCGCCGCACGTGCACCGGAACGGCGAGTGGTTCGGCCAGCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                             <  1:413396/150‑1 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCTCTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                             <  2:1563710/150‑1 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCg                                                                                                <  2:497100/99‑1 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCg                                                                                                >  1:497100/1‑99 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                             >  2:474973/1‑150 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                             >  2:2186847/1‑150 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                             <  1:238816/150‑1 (MQ=255)
                                                                                                          tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGAACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                             >  2:651562/1‑150 (MQ=255)
                                                                                                              ggAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCAc                                         <  2:285026/149‑1 (MQ=255)
                                                                                                               gAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGc                                                                                                      >  1:902711/1‑88 (MQ=255)
                                                                                                               gAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGc                                                                                                      <  2:902711/88‑1 (MQ=255)
                                                                                                                  gCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCTACGCTTGAGTGCGGGCACGCTa                                     >  2:737984/1‑150 (MQ=255)
                                                                                                                    cTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgac                                                              <  2:571571/123‑1 (MQ=255)
                                                                                                                    cTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgac                                                              >  1:571571/1‑123 (MQ=255)
                                                                                                                       gcgcGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGTGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTggg                                >  1:515980/1‑150 (MQ=255)
                                                                                                                          cgACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTc                             <  1:1677517/150‑1 (MQ=255)
                                                                                                                           gACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTc                                                                                                                               <  1:378576/51‑1 (MQ=255)
                                                                                                                           gACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTc                                                                                                                               >  2:378576/1‑51 (MQ=255)
                                                                                                                                ccGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGGGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTg                                  <  1:1949709/139‑1 (MQ=255)
                                                                                                                                ccGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTg                                  >  2:1949709/1‑139 (MQ=255)
                                                                                                                                 cGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATAGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCAcc                      <  2:1826988/150‑1 (MQ=255)
                                                                                                                                   agTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgt                    >  1:1178490/2‑150 (MQ=255)
                                                                                                                                    gTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGtt                                                                                             <  1:832741/76‑1 (MQ=37)
                                                                                                                                    gTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGtt                                                                                             >  2:832741/1‑76 (MQ=37)
                                                                                                                                    gTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgtg                   >  1:700268/1‑150 (MQ=255)
                                                                                                                                    gTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgtg                   >  1:1431531/1‑150 (MQ=255)
                                                                                                                                    gTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCACCCCCATGGCgccgccg                                                                                                         <  2:663609/64‑2 (MQ=16)
                                                                                                                                    gTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCACCCCCATGGCgccgccg                                                                                                         >  1:663609/1‑63 (MQ=18)
                                                                                                                                      gCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTc                 >  1:1758875/1‑150 (MQ=255)
                                                                                                                                        aCTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCa               >  2:1711070/1‑150 (MQ=255)
                                                                                                                                         cTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCAt                                                                                                                   >  1:415322/1‑49 (MQ=255)
                                                                                                                                         cTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCAt                                                                                                                   <  2:415322/49‑1 (MQ=255)
                                                                                                                                          ttGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAAc             >  1:710664/1‑150 (MQ=255)
                                                                                                                                            gCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGc           >  2:1564081/1‑150 (MQ=255)
                                                                                                                                               gCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGa        >  1:510905/1‑150 (MQ=255)
                                                                                                                                                   cACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAgcgc    >  2:1551975/1‑150 (MQ=255)
                                                                                                                                                     cGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAg  <  2:1498832/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAG  >  NC_002947/3638561‑3638859

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