Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,638,813 G→C 9.7% R270P (CGC→CCC)  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,8130GC 9.7% 215.9 / 12.7 72R270P (CGC→CCC) PP_3205fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base C (4/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGAC  >  NC_002947/3638673‑3638961
                                                                                                                                            |                                                                                                                                                    
cTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCCCacgctgc                                                                                                                                              >  3:225421/1‑143 (MQ=255)
   gAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCt                                                                                                                                           >  2:311544/1‑149 (MQ=255)
      cGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCgg                                                                                                                                        <  8:256294/149‑1 (MQ=255)
        gCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCCCACGCTGCGCTCGGCg                                                                                                                                      <  4:225421/149‑1 (MQ=255)
          gAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGTGCAGCGTGGGCTCGGCGTg                                                                                                                                    >  5:68482/1‑149 (MQ=255)
            gTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCa                                                                                                                                  >  1:6826/1‑149 (MQ=255)
               gTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCg                                                                                                                               >  8:84407/1‑149 (MQ=255)
                      cATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgcg                                                                                                                        >  6:161602/1‑149 (MQ=255)
                       aagtgcATGGCGCCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgcgc                                                                                                                       >  3:321845/6‑149 (MQ=255)
                           ccATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAgcgccgcg                                                                                                                   >  4:130038/1‑149 (MQ=255)
                                gcggcgCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGTGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATc                                                                                                              <  6:68482/149‑1 (MQ=255)
                                         tGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATc                                                                                                     <  3:130038/149‑1 (MQ=255)
                                           cACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCg                                                                                                    <  7:229677/148‑1 (MQ=255)
                                             cTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCTACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCACCGAGCGCCGCGCCATCGAAACCATCGCCt                                                                                                 <  1:30408/149‑1 (MQ=255)
                                              ttCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTg                                                                                                >  6:20727/1‑149 (MQ=255)
                                              ttCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCAGCGAAACCATCGCCTg                                                                                                >  6:339964/1‑149 (MQ=255)
                                                cGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTg                                                                                                <  6:114431/147‑1 (MQ=255)
                                                cGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTg                                                                                                <  7:84407/147‑1 (MQ=255)
                                                       caccaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCgatgc                                                                                                                                  <  1:341904/106‑5 (MQ=255)
                                                       caccaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCgatgc                                                                                                                                  >  2:341904/1‑102 (MQ=255)
                                                       caccaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATc                                                                                                     >  1:22742/1‑135 (MQ=255)
                                                       caccaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATc                                                                                                     <  2:22742/135‑1 (MQ=255)
                                                       caccaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGcc                                                                                       >  6:12740/1‑149 (MQ=255)
                                                       caccaGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGcc                                                                                       >  7:128462/1‑149 (MQ=255)
                                                            gTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCcacgca                                                                                  >  2:138845/1‑149 (MQ=255)
                                                               gATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATTGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACgg                                                                               >  6:249939/1‑149 (MQ=255)
                                                               gATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACgg                                                                               >  4:138742/1‑149 (MQ=255)
                                                                    cccATGCGGCGCTGACGCGACGCTTGGGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTc                                                                          <  6:128375/149‑1 (MQ=255)
                                                                    cccATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTc                                                                          <  5:12740/149‑1 (MQ=255)
                                                                     ccATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCa                                                                         >  8:53770/1‑149 (MQ=255)
                                                                          cggcgCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATgcgc                                                                    >  2:303968/1‑149 (MQ=255)
                                                                            gcgcTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCtt                                                                  <  1:303968/149‑1 (MQ=255)
                                                                              gcTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTg                                                                <  5:20727/149‑1 (MQ=255)
                                                                               cTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCCCACGCTGCGCTCGGCGttg                                                                                                                                   <  4:248078/81‑3 (MQ=25)
                                                                               cTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCCCACGCTGCGCTCGGCGttg                                                                                                                                   >  3:248078/1‑79 (MQ=25)
                                                                                  acgcgacgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGa                                                            >  3:26784/1‑149 (MQ=255)
                                                                                    gcgacgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGTAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGAcc                                                          >  4:267425/1‑149 (MQ=255)
                                                                                    gcgacgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGAcc                                                          >  3:228975/1‑149 (MQ=255)
                                                                                        cgcTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGt                                                      <  8:128462/149‑1 (MQ=255)
                                                                                              aGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCAt                                                <  6:87790/149‑1 (MQ=255)
                                                                                                  cGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCTGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCgg                                                        <  6:50054/137‑1 (MQ=255)
                                                                                                  cGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCTGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCgg                                                        >  5:50054/1‑137 (MQ=255)
                                                                                                  cGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGa                                             <  1:311544/148‑1 (MQ=255)
                                                                                                  cGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGa                                             <  3:138742/148‑1 (MQ=255)
                                                                                                  cGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAg                                            <  7:53770/149‑1 (MQ=255)
                                                                                                       aCGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTcg                                       >  5:256213/1‑149 (MQ=255)
                                                                                                         gCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCa                                                 <  4:38194/137‑1 (MQ=255)
                                                                                                         gCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCa                                                 >  3:38194/1‑137 (MQ=255)
                                                                                                          cTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCgg                                    <  5:161602/149‑1 (MQ=255)
                                                                                                                 gCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAg                             <  4:26784/149‑1 (MQ=255)
                                                                                                                 gCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAg                             <  4:321845/149‑1 (MQ=255)
                                                                                                                   tCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAggg                           >  8:165190/1‑149 (MQ=255)
                                                                                                                       gCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGTAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAg                       <  3:267425/149‑1 (MQ=255)
                                                                                                                          ccGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGGGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTg                    >  2:61672/1‑149 (MQ=255)
                                                                                                                            gtgtCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCt                  <  6:256213/149‑1 (MQ=255)
                                                                                                                              gtCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGt                >  8:182320/1‑149 (MQ=255)
                                                                                                                               tCCAACGCCGAGCGCAGCGTGGGCTCGGCGTTCATCGCCGAGCGCCGCTCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCATTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGATGtt               >  6:284060/1‑149 (MQ=255)
                                                                                                                               tCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATTGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGtt               <  5:249939/149‑1 (MQ=255)
                                                                                                                               tCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGa                        <  2:233280/140‑1 (MQ=255)
                                                                                                                               tCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGa                        >  1:233280/1‑140 (MQ=255)
                                                                                                                                 cAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAg                       <  7:298579/139‑1 (MQ=255)
                                                                                                                                 cAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAg                       >  8:298579/1‑139 (MQ=255)
                                                                                                                                    cGCCGAGCCCACGCTGCGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTgggg          >  6:125659/1‑149 (MQ=255)
                                                                                                                                    cGCCGAGCCCACGCTGCGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTgggg          <  5:125659/149‑1 (MQ=255)
                                                                                                                                     gCCGAGCCCACGCTGCGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGc         >  8:317803/1‑149 (MQ=255)
                                                                                                                                      ccGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGc         >  2:36229/1‑148 (MQ=255)
                                                                                                                                       cGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGc         <  7:131435/147‑1 (MQ=255)
                                                                                                                                       cGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGc         >  8:131435/1‑147 (MQ=255)
                                                                                                                                        gAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCGTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGa                              >  1:109718/1‑125 (MQ=255)
                                                                                                                                        gAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGa                              <  2:109718/125‑1 (MQ=255)
                                                                                                                                        gAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAt      >  8:228745/1‑149 (MQ=255)
                                                                                                                                         aGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAACCCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCAtt     <  4:40907/149‑1 (MQ=255)
                                                                                                                                            gcAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGAc  >  4:168500/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGGTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCAGCGTGGGCTCGGCGTGCATCGCCGAGCGCCGCGCCATCGAAACCATCGCCTGCGGTGCGCCACGCACGGCCTTCATGCGCTTTGGCGACCGGGTGCGCATGGAGGCTCGCGGGCCCGAGGGCGAGGTGCTGTTTGGGGCCATTGAC  >  NC_002947/3638673‑3638961

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