Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,457,228 T→C 6.9% I21V (ATC→GTC) ‡ hbdH ← 3‑hydroxybutyrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,457,2280TC 6.9% 205.2 / 4.8 72I21V (ATC→GTC) ‡hbdH3‑hydroxybutyrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/42);  new base C (2/3);  total (27/45)
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.03e-01

AAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAGGGTCTTGATCAGGGTAATACCAA  >  NC_002947/3457081‑3457372
                                                                                                                                                   |                                                                                                                                                
aaGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATg                                                                                                                                                 >  1:14732/1‑149 (MQ=255)
    cTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCa                                                                                                                                             <  4:3260/149‑1 (MQ=255)
    cTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCa                                                                                                                                             <  2:285457/149‑1 (MQ=255)
     tagATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAAACTTCACCCCGTGCCGGGCAATTTCCGCCATGGAGGGGCCCGGCTCACCAAAGCCGTTGCGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAg                                                                                                                                            >  2:348080/3‑149 (MQ=255)
       gATCTGGACCACATCCGATAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCATCGCCCAgg                                                                                                                                           <  6:173533/148‑1 (MQ=255)
       gATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAg                                                                                                                                            <  3:39723/147‑1 (MQ=255)
           tGGACCACATCCGACGGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGACAAt                                                                                                                                      >  2:30241/1‑149 (MQ=255)
             gACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAAt                                                                                                                                      <  8:176004/147‑1 (MQ=255)
             gACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGc                                                                                                                                    <  5:219012/149‑1 (MQ=255)
                cacaTCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGAt                                                                                                                                                  <  6:359020/132‑1 (MQ=255)
                cacaTCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGAt                                                                                                                                                  >  5:359020/1‑132 (MQ=255)
                cacaTCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGc                                                                                                                                 <  3:93884/149‑1 (MQ=255)
                   aTCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCtgg                                                                                                                              <  5:105250/149‑1 (MQ=255)
                     ccGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCtggtg                                                                                                                            <  4:241605/149‑1 (MQ=255)
                      cGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCtggtg                                                                                                                            <  8:212802/148‑1 (MQ=255)
                            gTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCgg                                                                                                                     >  3:53154/1‑149 (MQ=255)
                            gTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCgg                                                                                                                     <  7:268410/149‑1 (MQ=255)
                                  cGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGa                                                                                                               <  4:60645/149‑1 (MQ=255)
                                    ggtggtGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGt                                                                                                             <  8:164883/149‑1 (MQ=255)
                                        gtgCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGATCCGGTGACGAGTGCAg                                                                                                         >  5:200810/1‑149 (MQ=255)
                                        gtgCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCa                                                                                                          <  2:160686/148‑1 (MQ=255)
                                            acaACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCGGCTGGTGGAACCGGTGACGAGTGCAGtttt                                                                                                     <  1:32411/149‑1 (MQ=255)
                                            acaACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGtttt                                                                                                     >  2:273922/1‑149 (MQ=255)
                                            acaACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGtttt                                                                                                     <  3:207012/149‑1 (MQ=255)
                                            acaACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGtttt                                                                                                     <  1:273922/149‑1 (MQ=255)
                                              aaCCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGc                                                                                                   <  6:371422/149‑1 (MQ=255)
                                              aaCCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGc                                                                                                   <  1:27/149‑1 (MQ=255)
                                                     aCCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCaa                                                                                             <  5:101945/148‑1 (MQ=255)
                                                     aCCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAg                                                                                            >  6:332941/1‑149 (MQ=255)
                                                     aCCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGACAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAg                                                                                            <  3:169775/149‑1 (MQ=255)
                                                      ccccGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAgg                                                                                           <  3:327360/149‑1 (MQ=255)
                                                       cccGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAggg                                                                                          <  4:9794/149‑1 (MQ=255)
                                                        ccGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGtttt                                                                                                     <  1:146306/137‑1 (MQ=255)
                                                        ccGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGtttt                                                                                                     >  2:146306/1‑137 (MQ=255)
                                                             ccGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATg                                                                                     <  4:295804/148‑1 (MQ=255)
                                                                ggCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCgg                                                                                                                     <  3:310088/113‑1 (MQ=255)
                                                                ggCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCgg                                                                                                                     >  4:310088/1‑113 (MQ=255)
                                                                     tttCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTcc                                                                            <  1:176249/149‑1 (MQ=255)
                                                                     tttCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTcc                                                                            >  7:334886/1‑149 (MQ=255)
                                                                        ccGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCtt                                                                          <  3:3436/148‑1 (MQ=255)
                                                                         cGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTa                                                                         <  6:190113/148‑1 (MQ=255)
                                                                          gCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCa                                                                                                                                             <  2:78012/79‑1 (MQ=255)
                                                                          gCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCa                                                                                                                                             >  1:78012/1‑79 (MQ=255)
                                                                               ggCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATg                                                                  >  2:101693/1‑149 (MQ=255)
                                                                                      cccGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTgg                                                           <  3:101246/149‑1 (MQ=255)
                                                                                              aCCAACGCCGTTGAGCACGCTGTTGGCGCCCGCCCGGGCCAATACCTGGGCGATCCCCAGGCCAATGCCGCTGGGGGAACCGGTCACGAGGGCAGTGTTGCCTTCCAGGGTCATGCAATCCTCCTTAGACGATCCCAGTGGCGGAGAAc                                                   <  2:14732/149‑1 (MQ=255)
                                                                                               ccAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACg                                                  >  6:257684/1‑149 (MQ=255)
                                                                                                cAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGt                                                                                                                    <  7:293241/82‑1 (MQ=255)
                                                                                                cAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGt                                                                                                                    >  8:293241/1‑82 (MQ=255)
                                                                                                cAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGt                                                 >  4:48010/1‑149 (MQ=255)
                                                                                                 aaaGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTg                                                <  1:348080/149‑1 (MQ=255)
                                                                                                  aaGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTgg                                               >  3:306295/1‑149 (MQ=255)
                                                                                                  aaGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTgg                                               >  3:346668/1‑149 (MQ=255)
                                                                                                     ccGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCa                                                                                              >  2:201640/1‑99 (MQ=255)
                                                                                                     ccGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCa                                                                                              <  1:201640/99‑1 (MQ=255)
                                                                                                     ccGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGa                                            >  2:119560/1‑149 (MQ=255)
                                                                                                           aGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGAcc                                       <  5:257684/148‑1 (MQ=255)
                                                                                                              aCGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACg                                    >  7:33929/1‑148 (MQ=255)
                                                                                                                gATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGaa                                  >  7:212814/1‑148 (MQ=255)
                                                                                                                 aTGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAaca                                <  8:212814/149‑1 (MQ=255)
                                                                                                                          ccAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAggg                       >  5:309869/1‑149 (MQ=255)
                                                                                                                          ccAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGGGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAggg                       <  4:12536/149‑1 (MQ=255)
                                                                                                                                  ggCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGa                                         <  5:376177/123‑1 (MQ=255)
                                                                                                                                  ggCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGa                                         >  6:376177/1‑123 (MQ=255)
                                                                                                                                      aaTACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGNCCACGAACACCGCCAGGGTCTTGATCAggg           <  1:30241/149‑1 (MQ=255)
                                                                                                                                      aaTACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAGGGTCTTGATCAggg           <  8:334886/149‑1 (MQ=255)
                                                                                                                                         gccTGGGCATCGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGa                                            >  6:350461/2‑113 (MQ=255)
                                                                                                                                         gccTGGGCATCGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGa                                            <  5:350461/112‑1 (MQ=255)
                                                                                                                                          ccTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAGGGTCTTGATCAGGGTaa        <  1:101693/148‑1 (MQ=255)
                                                                                                                                          ccTGGGCATCGCCCAGGTCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGt                                                             >  3:142796/1‑95 (MQ=255)
                                                                                                                                          ccTGGGCATCGCCCAGGTCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGt                                                             <  4:142796/95‑1 (MQ=255)
                                                                                                                                           cTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAGGGTCTTGa                 <  8:266744/138‑1 (MQ=255)
                                                                                                                                           cTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAGGGTCTTGa                 >  7:266744/1‑138 (MQ=255)
                                                                                                                                                cGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAGGGTCTTGATCAGGGTAATACCaa  >  2:183867/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCCCGGGCCAATACCTGGGCGATGCCCAGGCCAATGCCGCTGGTGGAACCGGTGACGAGTGCAGTTTTGCCTTCAAGGGTCATGCAATGCTCCTTAGACGATGCCAGTGGCGTAGAACGTGGCGATGACCACGAACACCGCCAGGGTCTTGATCAGGGTAATACCAA  >  NC_002947/3457081‑3457372

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