Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,457,062 C→A 7.6% R76L (CGC→CTC)  hbdH ← 3‑hydroxybutyrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,457,0620CA 7.6% 234.1 / 6.5 74R76L (CGC→CTC) hbdH3‑hydroxybutyrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base A (4/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGC  >  NC_002947/3456915‑3457207
                                                                                                                                                          |                                                                                                                                                 
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                                                                                                                 >  2:236281/1‑148 (MQ=255)
                                  ccGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGAc                                                                                                                 <  1:236281/148‑1 (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                                                                                                 >  6:135792/1‑139 (MQ=255)
                                                           gTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCg                                                                                                 <  5:135792/139‑1 (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                                                                                    >  3:288609/1‑149 (MQ=255)
                                                              ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGca                                                                                    >  1:119402/1‑149 (MQ=255)
                                                                   tttCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACaaa                                                                                                                               <  2:278140/101‑1 (MQ=255)
                                                                     tCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGAc                                                                                                       <  7:263606/123‑1 (MQ=255)
                                                                     tCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGAc                                                                                                       >  8:263606/1‑123 (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                             <  4:93017/149‑1 (MQ=255)
                                                                                                                     ttgACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCaaa                             <  7:87517/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:218387/3‑149 (MQ=255)
                                                                                                                                              tggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGtt           >  4:184843/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                         <  3:233343/133‑1 (MQ=255)
                                                                                                                                                gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCg                         >  4:233343/1‑133 (MQ=255)
                                                                                                                                                 ccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgc        >  7:164817/1‑149 (MQ=255)
                                                                                                                                                      cTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAg   >  5:178362/1‑149 (MQ=255)
                                                                                                                                                       tCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACATCCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGc  >  1:346553/1‑149 (MQ=255)
                                                                                                                                                         gcgTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGc  <  3:183411/147‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                 
CCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGAC‑‑‑‑‑‑‑GCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGC  >  NC_002947/3456915‑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: ‑