Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,638,703 C→A 6.1% P219Q (CCG→CAG)  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,7030CA 6.1% 354.1 / 4.7 111P219Q (CCG→CAG) PP_3205fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (54/49);  new base A (4/3);  total (59/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 112 total aligned reads displayed.

CAACGATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAAC  >  NC_002947/3638554‑3638848
                                                                                                                                                      |                                                                                                                                                 
cAACGATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGcc                                                                                                                                                   >  2:1109268/1‑150 (MQ=255)
 aaCGATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCg                                                                                                                                                  >  1:1030137/1‑150 (MQ=255)
 aaCGATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCg                                                                                                                                                  >  1:1990539/1‑150 (MQ=255)
 aaCGATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCg                                                                                                                                                  >  2:1837667/1‑150 (MQ=255)
   cGATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCt                                                                                                                                                <  2:2176554/150‑1 (MQ=255)
    gATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTg                                                                                                                                               >  2:2168282/1‑150 (MQ=255)
      tGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCAGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCa                                                                                                                                             >  1:2105765/1‑150 (MQ=255)
           gCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGGGCACTTGCCGCTGCACGTGc                                                                                                                                        <  1:1210729/150‑1 (MQ=255)
             cTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCAc                                                                                                                                      <  1:841478/150‑1 (MQ=255)
              tGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACGTGCAGCGGCAAGTGCAcc                                                                                                                                     <  2:2117120/150‑2 (MQ=255)
               gCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACGTGCAGCGGCAAGTGCAccc                                                                                                                                    <  1:1024661/150‑3 (MQ=255)
                cGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACGTGCAGCGGCAAGTGCAcccg                                                                                                                                   <  1:1958372/150‑4 (MQ=255)
                      cTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGACAGGTTTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGc                                                                                                                             <  1:61837/150‑1 (MQ=255)
                        ggCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGa                                                                                                                           <  1:384804/150‑1 (MQ=255)
                           actgcgGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGTCGCTGCACGTGCACTGGAACGGCGAGTg                                                                                                                        >  1:2427601/4‑150 (MQ=255)
                             cgcgGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGt                                                                                                                      <  1:1144210/150‑1 (MQ=255)
                               cgGGAAATGAAGACCGGCTTCGGCTTGCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCCCCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTc                                                                                                                    <  2:1343699/150‑1 (MQ=255)
                                gggAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCg                                                                                                                   >  2:1458500/1‑150 (MQ=255)
                                   aaaTGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGcc                                                                                                                >  1:446394/1‑150 (MQ=255)
                                   aaaTGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGcc                                                                                                                >  1:446942/1‑150 (MQ=255)
                                   aaaTGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGcc                                                                                                                >  1:453518/1‑150 (MQ=255)
                                       gAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATcc                                                                                                            >  2:1562232/1‑150 (MQ=255)
                                           aCCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCAt                                                                                                        <  2:1441738/150‑1 (MQ=255)
                                           aCCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCAt                                                                                                        <  2:1441367/150‑1 (MQ=255)
                                                   cggcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGcgccgc                                                                                                <  2:1030137/150‑1 (MQ=255)
                                                     gcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCa                                                                                              >  1:915118/1‑150 (MQ=255)
                                                       tttCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                            <  2:446942/150‑1 (MQ=255)
                                                       tttCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                            <  2:453518/150‑1 (MQ=255)
                                                       tttCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                            <  2:446394/150‑1 (MQ=255)
                                                         tctcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGTCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGCCGCCATggg                                                                                          <  2:2020832/150‑3 (MQ=255)
                                                           tcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                            <  1:1724821/146‑1 (MQ=255)
                                                           tcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                            >  2:1724821/1‑146 (MQ=255)
                                                           tcCAAGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGCCCGAACTCGGCGAGGCCTGGCGCGACGGTCGGGGGCACGTGACGCTGGACCCTCCCTGGAACCGCGAGTGGTTCGGCCATCCCCATGGCTGCGCCGTGCAcc                                                                                        >  2:1146213/1‑149 (MQ=255)
                                                             cAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTc                                                                                      >  1:2130434/1‑150 (MQ=255)
                                                             cAGGCCAAGCCCGCCAGCAGTTTCGCCCCCGTTGACGGTGGCACC‑ACCAAGCCCTGCGATGCCTGGCGCGACAAGTGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTc                                                                                      <  1:2378011/150‑1 (MQ=255)
                                                              aGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACtt                                                                                       >  1:837027/1‑148 (MQ=255)
                                                              aGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACtt                                                                                       <  2:837027/148‑1 (MQ=255)
                                                                  cAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGttt                                                                                 >  1:1624875/1‑150 (MQ=255)
                                                                     gCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCAt                                                                                                        >  2:1271516/1‑124 (MQ=255)
                                                                     gCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCAt                                                                                                        <  1:1271516/124‑1 (MQ=255)
                                                                        cgacagcagTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTccacca                                                                           <  2:43593/147‑1 (MQ=255)
                                                                        cGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTccacca                                                                           <  1:2445799/150‑1 (MQ=255)
                                                                        cGCAAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCACCACCGGCGCCATGCACTTCGGTTTccacca                                                                           <  1:1146213/150‑1 (MQ=255)
                                                                          ccAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCAt                                                                                             <  1:2263623/130‑1 (MQ=255)
                                                                          ccAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCAt                                                                                             >  2:2263623/1‑130 (MQ=255)
                                                                           cagcagTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGtt                                                                        <  1:675643/150‑1 (MQ=255)
                                                                            agcagTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGAAAGACCGGTGCACTTGTTTCGTTTCGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTg                                                                       >  1:2123418/1‑150 (MQ=255)
                                                                                gTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTg                                                                   <  1:2338761/150‑1 (MQ=255)
                                                                                 tttCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGc                                                                  <  2:2130434/150‑1 (MQ=255)
                                                                                      cgcCCGTTG‑CGGTCACCCAGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATg                                                             <  2:2105765/150‑1 (MQ=255)
                                                                                       gcCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATgc                                                            >  1:2239615/1‑150 (MQ=255)
                                                                                         ccGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATgcgg                                                          >  2:1595439/1‑150 (MQ=255)
                                                                                           gTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACgg                                                                                                                              >  1:1549513/1‑80 (MQ=255)
                                                                                           gTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACgg                                                                                                                              <  2:1549513/80‑1 (MQ=255)
                                                                                            ttG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGGACTGGAACGGCGAGTgg                                                                                                                       >  1:1355570/1‑86 (MQ=255)
                                                                                            ttG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTgg                                                                                                                       <  2:1355570/86‑1 (MQ=255)
                                                                                             tG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCt                                                      <  1:2168282/150‑1 (MQ=255)
                                                                                                 ggTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgac                                                   <  2:2400798/150‑1 (MQ=255)
                                                                                                 ggTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgac                                                   <  2:1368962/150‑1 (MQ=255)
                                                                                                  gTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacg                                                  >  1:304764/1‑150 (MQ=255)
                                                                                                   tCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCg                                                                                     <  2:114008/114‑1 (MQ=255)
                                                                                                   tCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCg                                                                                     >  1:114008/1‑114 (MQ=255)
                                                                                                   tCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacgc                                                 >  1:72101/1‑150 (MQ=255)
                                                                                                       cccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacgcgacg                                             >  1:1941569/1‑150 (MQ=255)
                                                                                                       cccGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTgacgcgacg                                             >  1:1915541/1‑150 (MQ=255)
                                                                                                          ggACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCtt                                          <  2:2239615/150‑1 (MQ=255)
                                                                                                                aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCa                                                                                                         >  1:201885/1‑81 (MQ=255)
                                                                                                                aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCa                                                                                                         <  2:201885/81‑1 (MQ=255)
                                                                                                                aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCAACCCGCTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGc                                    >  2:1145444/1‑150 (MQ=255)
                                                                                                                 cTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCg                                   >  2:2441085/1‑150 (MQ=255)
                                                                                                                  tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATgcggc                                                         <  1:1417858/127‑1 (MQ=255)
                                                                                                                  tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATgcggc                                                         >  2:1417858/1‑127 (MQ=255)
                                                                                                                  tCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                  >  2:1091252/1‑150 (MQ=255)
                                                                                                                   cGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCggg                                 >  2:2414437/1‑150 (MQ=255)
                                                                                                                    ggCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGc                                >  2:579938/1‑150 (MQ=255)
                                                                                                                    ggCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGc                                >  1:1806523/1‑150 (MQ=255)
                                                                                                                       gAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACg                             >  2:1889355/1‑150 (MQ=255)
                                                                                                                          gCCTGGCGAGACGGCCGGGTGCAATTGGAGCTGCACGTGCGCTGGAAGGGCGAGTGGTTCGGGCATCCCCATGGCGGCGCCATGCACTTAGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTa                          >  2:1048183/1‑150 (MQ=255)
                                                                                                                            cTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACt                        >  2:249988/1‑150 (MQ=255)
                                                                                                                            cTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACt                        >  1:1312844/1‑150 (MQ=255)
                                                                                                                            cTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACt                        <  2:2445275/150‑1 (MQ=255)
                                                                                                                                 gcgACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCAc                                                                                         <  1:1756734/80‑1 (MQ=255)
                                                                                                                                 gcgACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCAc                                                                                         >  2:1756734/1‑80 (MQ=255)
                                                                                                                                   gACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTCGGCTcc                 >  2:44788/1‑150 (MQ=255)
                                                                                                                                    aCGGCCGGGTGCACTTGTCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCg                <  2:2427601/150‑1 (MQ=255)
                                                                                                                                     cGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGAGGTTCGGCCATCCCCATGGCAGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTAGGCTCCgg               <  1:511882/150‑1 (MQ=255)
                                                                                                                                      gttccagTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGc              >  1:2269952/7‑150 (MQ=255)
                                                                                                                                       gCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCtt                                          >  2:671040/1‑121 (MQ=255)
                                                                                                                                           agTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACACGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgt         >  1:2380783/2‑150 (MQ=255)
                                                                                                                                            gTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                  >  1:1603400/1‑124 (MQ=255)
                                                                                                                                            gTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                  <  2:1603400/124‑1 (MQ=255)
                                                                                                                                            gTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgtg        >  2:1369961/1‑150 (MQ=255)
                                                                                                                                            gTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACAGGGCTCCGGCACCgtg        >  2:1222230/1‑150 (MQ=255)
                                                                                                                                            gTGCACGTGCAGCGGCAAGTGCCATGTGATGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTCGGCTCCGGCACCgtg        >  1:918563/1‑150 (MQ=255)
                                                                                                                                             ggCACGTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCCTGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgtgt       <  1:991767/149‑1 (MQ=255)
                                                                                                                                                aCTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCCCCAGTTGATTGCCCa                                                               >  2:2412945/1‑91 (MQ=255)
                                                                                                                                                aCTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCa                                                               <  1:2412945/91‑1 (MQ=255)
                                                                                                                                                 cTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCt                                                      <  2:1119154/99‑1 (MQ=255)
                                                                                                                                                 cTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCt                                                      >  1:1119154/1‑99 (MQ=255)
                                                                                                                                                  ttGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAAc  >  2:2075959/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CAACGATGTCAGCCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTG‑CGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAAC  >  NC_002947/3638554‑3638848

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