Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,457,102 T→G 5.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,1020TG 5.6% 200.3 / 4.7 71E77A (GAG→GCG) hbdH3‑hydroxybutyrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base G (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

GCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCC  >  NC_002947/3456957‑3457251
                                                                                                                                                 |                                                                                                                                                     
gCCCGGTAGCGCAAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGTTCTTGTCCCAGCTTTCTGGCGGGAACTGCTTAACCGGTGCCACTTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTcgcg                                                                                                                                                   <  2:132652/150‑1 (MQ=255)
  ccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGtt                                                                                                                                                 >  1:33334/1‑150 (MQ=255)
  ccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGtt                                                                                                                                                 >  1:1649069/1‑150 (MQ=255)
  ccGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGtt                                                                                                                                                 >  1:1621585/1‑150 (MQ=255)
    ggCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACGAGGATGTCGACGCCGCCGAACTCGCGTTCg                                                                                                                                               >  1:1376349/1‑150 (MQ=255)
       aGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATACTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACGCCAGTTCGGcg                                                                                                                                            <  2:181692/150‑2 (MQ=255)
         cgcCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAg                                                                                                                                          <  1:1587950/150‑1 (MQ=255)
                   gtcccATGCAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACaaa                                                                                                                                <  2:665932/147‑1 (MQ=255)
                    ggCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAg                                                                                                                               <  1:421908/149‑1 (MQ=255)
                       cATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAg                                                                                                                               <  2:250132/147‑1 (MQ=255)
                       cATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAg                                                                                                                               >  1:250132/1‑147 (MQ=255)
                       cATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACGCGAGTTCgg                                                                                                                                              >  1:1258912/1‑132 (MQ=255)
                       cATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACGCGAGTTCgg                                                                                                                                              <  2:1258912/132‑1 (MQ=255)
                        aTGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTc                                                                                                                           >  1:684965/1‑150 (MQ=255)
                                 ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTcgc                                                                                                                                                    >  2:259925/1‑116 (MQ=255)
                                 ggCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTcgc                                                                                                                                                    <  1:259925/116‑1 (MQ=255)
                                     gACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGATAGGACAGCTTGAAAACAGCCTCGGTCAGGAcctta                                                                                                              >  2:416993/1‑147 (MQ=255)
                                       cAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATc                                                                                                            <  1:254649/150‑1 (MQ=255)
                                             cAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCAAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAg                                                                                                      <  2:178679/150‑1 (MQ=255)
                                              aGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAgg                                                                                                     <  1:1638928/150‑1 (MQ=255)
                                                ggCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTc                                                                                                   >  2:1219993/1‑150 (MQ=255)
                                                   gatgatCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATcc                                                                                                           <  1:774705/139‑1 (MQ=255)
                                                   gatgatCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATcc                                                                                                           >  2:774705/1‑139 (MQ=255)
                                                   gatgatCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATACTGGAGACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGc                                                                                                <  2:1717650/150‑1 (MQ=255)
                                                       atCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCggg                                                                                            >  1:281731/1‑150 (MQ=255)
                                                        tCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTCAACAACGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggt                                                                                           <  1:1340069/150‑1 (MQ=255)
                                                          ttGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggtgg                                                                                         <  2:281731/150‑1 (MQ=255)
                                                           tGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGggtggt                                                                                        <  2:1436081/150‑1 (MQ=255)
                                                             tCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGc                                                                                      >  1:788715/1‑150 (MQ=255)
                                                              cccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGca                                                                                     >  1:1167953/1‑150 (MQ=255)
                                                               ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCg                                                                                                                                               >  2:590603/1‑91 (MQ=255)
                                                               ccAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCg                                                                                                                                               <  1:590603/91‑1 (MQ=255)
                                                                  gCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAAc                                                                                 >  2:1743533/1‑150 (MQ=255)
                                                                              gAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTg                                                                     >  2:1515165/1‑150 (MQ=255)
                                                                                  tGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACaaa                                                                                                                                <  2:1717879/87‑1 (MQ=255)
                                                                                  tGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACaaa                                                                                                                                >  1:1717879/1‑87 (MQ=255)
                                                                                    cTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGc                                                               >  1:1062711/1‑150 (MQ=255)
                                                                                    cTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGc                                                               >  2:1873044/1‑150 (MQ=255)
                                                                                       aaCCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGACCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAAt                                                            <  2:488939/150‑1 (MQ=255)
                                                                                           ggTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAATAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTACCGGGCAATTTcc                                                        <  1:1219993/150‑1 (MQ=255)
                                                                                              gCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGcc                                                     >  2:1429500/1‑150 (MQ=255)
                                                                                               ccACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCa                                                    >  2:1821517/1‑150 (MQ=255)
                                                                                               ccACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCa                                                    >  2:446484/1‑150 (MQ=255)
                                                                                                cacaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCAt                                                   >  2:1587122/1‑150 (MQ=255)
                                                                                                 acaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCAt                                                   >  1:868261/1‑149 (MQ=255)
                                                                                                 acaTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCAt                                                   <  2:868261/149‑1 (MQ=255)
                                                                                                   aTGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGc                                                <  1:1515165/150‑1 (MQ=255)
                                                                                                            gccGGCGTTGTTGACCGGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCgg                                       <  2:536058/149‑1 (MQ=255)
                                                                                                              cGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGc                                                                                                                              <  2:1172494/61‑1 (MQ=255)
                                                                                                              cGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGc                                                                                                                              >  1:1172494/1‑61 (MQ=255)
                                                                                                              cGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCgg                                                                 <  2:1085531/122‑1 (MQ=255)
                                                                                                              cGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCgg                                                                 >  1:1085531/1‑122 (MQ=255)
                                                                                                                     gttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCaaa                              <  2:788715/150‑1 (MQ=255)
                                                                                                                      ttgACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAg                             >  1:798301/1‑150 (MQ=255)
                                                                                                                          ccAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGt                         >  2:439534/1‑150 (MQ=255)
                                                                                                                           cAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGtt                        >  1:1747879/1‑150 (MQ=255)
                                                                                                                             ggATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGa                      <  2:1649069/150‑1 (MQ=255)
                                                                                                                             ggATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGGCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCATGTTTCACCAAAGCCGTTGa                      >  1:1069316/1‑150 (MQ=255)
                                                                                                                                tGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCAGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCa                   >  2:644106/1‑150 (MQ=255)
                                                                                                                                   caacctggccgAACGCGAGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCCCCCCCGGAACATCAATTTCCCCCATGTCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGa                >  1:1400638/8‑150 (MQ=255)
                                                                                                                                     aCGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATg              <  1:1743533/150‑1 (MQ=255)
                                                                                                                                        tggccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTg           >  1:1184685/3‑150 (MQ=255)
                                                                                                                                         cgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTgg          >  1:1435678/1‑150 (MQ=255)
                                                                                                                                         cgccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTgg          >  1:1333677/1‑150 (MQ=255)
                                                                                                                                          gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCggg                                                                <  2:1136569/95‑1 (MQ=255)
                                                                                                                                          gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCggg                                                                >  1:1136569/1‑95 (MQ=255)
                                                                                                                                          gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgc         <  1:1821517/150‑1 (MQ=255)
                                                                                                                                          gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgc         <  2:1184685/150‑1 (MQ=255)
                                                                                                                                          gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGAACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCgg                                                                 <  1:1360642/94‑1 (MQ=255)
                                                                                                                                          gccgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGAACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCgg                                                                 >  2:1360642/1‑94 (MQ=255)
                                                                                                                                            cgAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGgcgc       <  1:416993/150‑1 (MQ=255)
                                                                                                                                                 tCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGcc  <  2:1062711/150‑1 (MQ=255)
                                                                                                                                                 tCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGcc  >  1:1739189/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCCCGGCAGCGCCAGGCGCGTGCCATGGAATACGGCCGACAGGTTCAGGGCGATGATCTTGTCCCAGCTTTCTGGCGGGAACTGCTCAACCGGTGCCACATGCTGGATACCGGCGTTGTTGACCAGGATGTCGACGCCGCCGAACTCGCGTTCGGCCAGGTTGAACAAAGCCTCGATCTGGACCACATCCGACAGGTCGGCCGGGTGGTGCACAACCTTCACCCCGTGCCGGGCAATTTCCGCCATGGCGGGGCCCGGGTCACCAAAGCCGTTGAGCACGATGTTGGCGCCAGCC  >  NC_002947/3456957‑3457251

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