Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,638,816 G→C 6.0% S271T (AGC→ACC) ‡ 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,8160GC 6.0% 218.4 / 3.8 67S271T (AGC→ACC) ‡PP_3205fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (30/33);  new base C (2/2);  total (32/35)
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.80e-01

GTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAG  >  NC_002947/3638668‑3638964
                                                                                                                                                     |                                                                                                                                                     
gTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAg                                                                                                                                                       >  3:77200/1‑149 (MQ=255)
  gCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCg                                                                                                                                                     >  7:293173/1‑149 (MQ=255)
    aCTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCACCGCCGAGCGCAGCGTg                                                                                                                                                   >  5:273756/1‑149 (MQ=255)
      tGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTggg                                                                                                                                                 >  7:90792/1‑149 (MQ=255)
        gAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTggg                                                                                                                                                 <  8:293173/147‑1 (MQ=255)
          aCGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCg                                                                                                                                             >  4:181536/1‑149 (MQ=255)
          aCGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCg                                                                                                                                             >  2:109229/1‑149 (MQ=255)
           cGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCgg                                                                                                                                            >  6:256567/1‑149 (MQ=255)
            ggCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCCCACGCTGCGCTCGGc                                                                                                                                           <  6:153209/149‑1 (MQ=255)
              cGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGt                                                                                                                                        <  4:284481/149‑1 (MQ=255)
              cGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGt                                                                                                                                        <  7:218543/149‑1 (MQ=255)
                 gTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCgatgc                                                                                                                                      <  8:95530/149‑5 (MQ=255)
                    gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGc                                                                                                                                                      <  2:35670/130‑1 (MQ=255)
                    gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTg                                                                                                                                                   <  7:79451/133‑1 (MQ=255)
                    gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTg                                                                                                                                                   >  8:79451/1‑133 (MQ=255)
                       cGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCg                                                                                                                               >  6:93739/1‑149 (MQ=255)
                          ccATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGGCGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgc                                                                                                                            >  7:237658/1‑149 (MQ=255)
                          ccATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGATGCACGCCGAgc                                                                                                                            <  8:114639/149‑1 (MQ=255)
                          ccATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgc                                                                                                                            >  1:53543/1‑149 (MQ=255)
                             tCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgcgcc                                                                                                                         <  1:330677/149‑1 (MQ=255)
                                ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCATCGCCGAGCCCACGCTGCGCTCGGCGtt                                                                                                                                       >  7:162743/1‑131 (MQ=255)
                                ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCATCGCCGAGCCCACGCTGCGCTCGGCGtt                                                                                                                                       <  8:162743/132‑2 (MQ=255)
                                ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgcgccgcg                                                                                                                      >  6:185950/1‑149 (MQ=255)
                                         cgcCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGCTCAGTTGAGTGCTGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGaa                                                                                                              <  1:279463/149‑1 (MQ=255)
                                            cATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCa                                                                                                          >  5:44828/1‑149 (MQ=255)
                                                 aCTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGcc                                                                                                     >  6:87262/1‑149 (MQ=255)
                                                   ttCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTg                                                                                                   >  6:330886/1‑149 (MQ=255)
                                                     cGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTg                                                                                                   <  8:90792/147‑1 (MQ=255)
                                                             accaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGc                                                                                                  <  5:251581/140‑1 (MQ=255)
                                                             accaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGc                                                                                                  >  6:251581/1‑140 (MQ=255)
                                                             accaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                         >  8:44427/1‑149 (MQ=255)
                                                             accaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                         >  7:76511/1‑149 (MQ=255)
                                                             accaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                         >  7:158155/1‑149 (MQ=255)
                                                             accaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                         >  3:26459/1‑149 (MQ=255)
                                                             accaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGACCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCca                                                                                         <  1:109229/149‑1 (MQ=255)
                                                                 gTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCcacgca                                                                                     >  5:271471/1‑149 (MQ=255)
                                                                         cccATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCCCGCACGGCCtt                                                                              >  2:248718/1‑148 (MQ=255)
                                                                         cccATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCtt                                                                              <  6:44828/148‑1 (MQ=255)
                                                                            aTGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGGGGTGCGCCACGCACGGCCTTCATg                                                                          <  1:304428/149‑1 (MQ=255)
                                                                            aTGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGGCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGGGGTGCGCCACGCACGGCCTTCATg                                                                          <  1:308424/149‑1 (MQ=255)
                                                                             tGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATgc                                                                         <  7:44427/149‑1 (MQ=255)
                                                                               cggcgCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATgcgc                                                                       <  5:185950/149‑1 (MQ=255)
                                                                                  cgcTGGCGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCttt                                                                    <  8:237658/149‑1 (MQ=255)
                                                                                       acgcgacgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGc                                                                 <  8:76511/147‑1 (MQ=255)
                                                                                       acgcgacgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGc                                                                 <  8:158155/147‑1 (MQ=255)
                                                                                               gcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTg                                                        >  4:24213/1‑149 (MQ=255)
                                                                                               gcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTg                                                        >  3:125727/1‑149 (MQ=255)
                                                                                                 ttGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACAGGGTgcg                                                      >  1:2815/1‑149 (MQ=255)
                                                                                                   gAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTgcgc                                                     <  4:26459/148‑1 (MQ=255)
                                                                                                        cGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGTCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGACTTCATGCGCTTTGGCGACCGGGTGCGCATGGa                                                <  5:93739/148‑1 (MQ=255)
                                                                                                         gggCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAgg                                              <  5:90691/149‑1 (MQ=255)
                                                                                                               gCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGCGCTCGGCGATGCACGCCGAgccca                                                                                                                         <  3:273120/68‑4 (MQ=25)
                                                                                                               gCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGCGCTCGGCGATGCACGCCGAgccca                                                                                                                         >  4:273120/1‑65 (MQ=25)
                                                                                                                   cTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGc                                     <  5:256567/148‑1 (MQ=255)
                                                                                                                        cTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGTGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCAGGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAgg                               <  2:2815/149‑1 (MQ=255)
                                                                                                                        cTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAgg                               >  5:330509/1‑149 (MQ=255)
                                                                                                                        cTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAgg                               <  6:330509/149‑1 (MQ=255)
                                                                                                                             gCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATCGAGGCTCGCGGGCCCGAGGGCGAg                          >  5:193860/1‑149 (MQ=255)
                                                                                                                                ccGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTg                       <  4:77200/149‑1 (MQ=255)
                                                                                                                                 cGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGc                      >  5:313931/1‑149 (MQ=255)
                                                                                                                                       cAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTg                >  2:331554/1‑149 (MQ=255)
                                                                                                                                        aaCGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTgg               >  8:179464/1‑149 (MQ=255)
                                                                                                                                          cGCCGAGCCCACGCTGCGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTgggg             >  4:310483/1‑149 (MQ=255)
                                                                                                                                           gCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGc            <  5:330886/149‑1 (MQ=255)
                                                                                                                                              gAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAt         <  2:53543/149‑1 (MQ=255)
                                                                                                                                              gAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAt         <  7:179464/149‑1 (MQ=255)
                                                                                                                                               aGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGGGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAtt        <  3:24213/149‑1 (MQ=255)
                                                                                                                                                  gcAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGAc     <  4:125727/149‑1 (MQ=255)
                                                                                                                                                     gCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAg  >  1:284586/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGACCAG  >  NC_002947/3638668‑3638964

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