Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,457,059 T→G 6.6% E77A (GAG→GCG)  hbdH ← 3‑hydroxybutyrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,457,0590TG 6.6% 237.8 / 7.5 76E77A (GAG→GCG) hbdH3‑hydroxybutyrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base G (3/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

ATGCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGC  >  NC_002947/3456912‑3457207
                                                                                                                                                          |                                                                                                                                                    
aTGCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTc                                                                                                                                                     <  4:213329/149‑1 (MQ=255)
   cccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTcgc                                                                                                                                                   <  7:276838/148‑1 (MQ=255)
   cccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTcgc                                                                                                                                                   <  2:147894/148‑1 (MQ=255)
    ccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGt                                                                                                                                                 <  6:110431/149‑1 (MQ=255)
         aGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCGAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACGCGAGTTCGGc                                                                                                                                            <  5:1633/149‑1 (MQ=255)
                 gcgcGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTg                                                                                                                                     >  5:263766/1‑148 (MQ=255)
                  cgcgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATTCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGaa                                                                                                                                   <  6:79162/149‑1 (MQ=255)
                  cgcgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGaa                                                                                                                                   <  2:205629/149‑1 (MQ=255)
                  cgcgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCCGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGaa                                                                                                                                   <  1:199395/149‑1 (MQ=255)
                   gcgTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGaa                                                                                                                                   >  1:55581/1‑148 (MQ=255)
                              aaTACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGAt                                                                                                                       <  3:141901/149‑1 (MQ=255)
                               aTACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATc                                                                                                                      <  6:263766/149‑1 (MQ=255)
                                  cGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTgg                                                                                                                   <  2:138384/149‑1 (MQ=255)
                                   ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                                  <  3:252246/149‑1 (MQ=255)
                                     ccGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAcc                                                                                                                <  3:4212/149‑1 (MQ=255)
                                     ccGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAc                                                                                                                 <  1:236281/148‑1 (MQ=255)
                                     ccGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAc                                                                                                                 >  2:236281/1‑148 (MQ=255)
                                        aCAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGAt                                                                                                                       <  4:357613/139‑1 (MQ=255)
                                        aCAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGAt                                                                                                                       >  3:357613/1‑139 (MQ=255)
                                            gTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCg                                                                                                                         >  4:236176/1‑133 (MQ=255)
                                            gTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATcc                                                                                                          >  7:182755/1‑148 (MQ=255)
                                                 gggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACGCGAGTTCGGcg                                                                                                                                           <  2:87874/110‑2 (MQ=255)
                                                  ggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGt                                                                                                   <  2:55581/149‑1 (MQ=255)
                                                          tCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCggg                                                                                           >  3:269935/1‑149 (MQ=255)
                                                             tGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGGGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggtg                                                                                         <  7:224434/148‑1 (MQ=255)
                                                              gTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCg                                                                                                 <  5:135792/139‑1 (MQ=255)
                                                              gTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCg                                                                                                 >  6:135792/1‑139 (MQ=255)
                                                               tCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGAAGGGGGGCGCCGAAAAACACTACGGCCAGGATGAACAAACACTCGATCATGACCACATCCTACAGGTCGGCCg                                                                                             >  1:297526/1‑149 (MQ=255)
                                                                 ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                                  >  3:119290/1‑119 (MQ=255)
                                                                 ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGa                                                                                                                  <  4:119290/119‑1 (MQ=255)
                                                                 ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGca                                                                                    >  1:119402/1‑149 (MQ=255)
                                                                 ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGca                                                                                    >  3:288609/1‑149 (MQ=255)
                                                                      tttCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACaaa                                                                                                                               <  2:278140/101‑1 (MQ=255)
                                                                        tCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGAc                                                                                                       >  8:263606/1‑123 (MQ=255)
                                                                        tCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGAc                                                                                                       <  7:263606/123‑1 (MQ=255)
                                                                              gggAACTGCTCAACTGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCAcccc                                                                       >  3:222359/1‑149 (MQ=255)
                                                                              gggAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCAcccc                                                                       >  6:151011/1‑149 (MQ=255)
                                                                                    tGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAc                                                                                                                 <  6:301407/101‑1 (MQ=255)
                                                                                    tGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAc                                                                                                                 >  5:301407/1‑101 (MQ=255)
                                                                                    tGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGAAGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGATCGGCCAGGTAGAACAAAGCCTCGATCTGGACCACATCCGACAGGACGGCCGGGTGGTGGCCAACCTTCACCCCGTGCCg                                                                 >  4:210324/1‑149 (MQ=255)
                                                                                           ccGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCt                                                                                                                           <  5:360390/84‑1 (MQ=255)
                                                                                           ccGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCt                                                                                                                           >  6:360390/1‑84 (MQ=255)
                                                                                              gTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTcc                                                       >  7:321724/1‑149 (MQ=255)
                                                                                               tGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCg                                                      <  4:269935/149‑1 (MQ=255)
                                                                                                  cacaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCa                                                   >  8:87517/1‑149 (MQ=255)
                                                                                                  cacaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCATTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGcc                                                                                              >  5:124336/1‑106 (MQ=255)
                                                                                                  cacaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCATTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGcc                                                                                              <  6:124336/106‑1 (MQ=255)
                                                                                                   agaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGCC‑‑‑‑‑‑‑GCCGCCGAACTCGCGGTCGGGCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGGCAGGTCGGCCGGGTGGTGCACAAGCTTCCCCACGAGCCGGGGAAGTTCCGGCAt                                                  >  6:69711/3‑149 (MQ=255)
                                                                                                         tGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGc                                                                                     >  1:317846/1‑108 (MQ=255)
                                                                                                         tGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGc                                                                                     <  2:317846/108‑1 (MQ=255)
                                                                                                         tGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCggg                                            >  5:68442/1‑149 (MQ=255)
                                                                                                            aTACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGcc                                         <  6:68442/149‑1 (MQ=255)
                                                                                                             tACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGcc                                         <  8:321724/148‑1 (MQ=255)
                                                                                                               ccGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                      <  8:182755/149‑1 (MQ=255)
                                                                                                               ccGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                      <  2:119402/149‑1 (MQ=255)
                                                                                                                    gttgttGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCAc                                 <  4:288609/149‑1 (MQ=255)
                                                                                                                     ttgttgACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCAc                                 >  3:93017/1‑148 (MQ=255)
                                                                                                                        ttgACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCaaa                             <  7:87517/149‑1 (MQ=255)
                                                                                                                        ttgACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCaaa                             <  4:93017/149‑1 (MQ=255)
                                                                                                                         tgACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCaaa                             <  2:293730/148‑1 (MQ=255)
                                                                                                                          gACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGc                                                                                     <  3:332970/91‑1 (MQ=255)
                                                                                                                          gACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGc                                                                                     >  4:332970/1‑91 (MQ=255)
                                                                                                                            ccAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCCGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCg                         >  7:229933/1‑149 (MQ=255)
                                                                                                                              aGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                       <  5:151011/149‑1 (MQ=255)
                                                                                                                               ggATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                       >  3:300328/1‑148 (MQ=255)
                                                                                                                                      gAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACg                <  1:326192/148‑1 (MQ=255)
                                                                                                                                       aC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGAt              <  4:222359/149‑1 (MQ=255)
                                                                                                                                                 tggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGtt           >  4:184843/3‑149 (MQ=255)
                                                                                                                                                 tggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGtt           >  4:218387/3‑149 (MQ=255)
                                                                                                                                                 tggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCAGGTCACCAAAGCCGTTGAGCACGATGtt           >  6:56305/3‑149 (MQ=255)
                                                                                                                                                  cgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggtg                                                                                         <  8:196304/70‑1 (MQ=255)
                                                                                                                                                  cgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTg          >  7:11137/1‑149 (MQ=255)
                                                                                                                                                   gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCg                         >  4:233343/1‑133 (MQ=255)
                                                                                                                                                   gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCg                         <  3:233343/133‑1 (MQ=255)
                                                                                                                                                    ccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgc        >  7:164817/1‑149 (MQ=255)
                                                                                                                                                         cTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAg   >  5:178362/1‑149 (MQ=255)
                                                                                                                                                          tCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACATCCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGc  >  1:346553/1‑149 (MQ=255)
                                                                                                                                                          |                                                                                                                                                    
ATGCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGC  >  NC_002947/3456912‑3457207

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